Spacetime geometries and light trapping in travelling refractive index perturbations
Identifieur interne : 001362 ( Istex/Corpus ); précédent : 001361; suivant : 001363Spacetime geometries and light trapping in travelling refractive index perturbations
Auteurs : S L Cacciatori ; F. Belgiorno ; V. Gorini ; G. Ortenzi ; L. Rizzi ; V G Sala ; D. FaccioSource :
- New Journal of Physics [ 1367-2630 ] ; 2010.
Abstract
In the framework of transformation optics, we show that the propagation of a locally superluminal refractive index perturbation (RIP) in a Kerr medium can be described, in the eikonal approximation, by means of a stationary metric, which we prove to be of Gordon type. Under suitable hypotheses on the RIP, we obtain a stationary but not static metric, which is characterized by an ergosphere and by a peculiar behaviour of the geodesics, which are studied numerically, also accounting for material dispersion. Finally, the equation to be satisfied by an event horizon is also displayed and briefly discussed.
Url:
DOI: 10.1088/1367-2630/12/9/095021
Links to Exploration step
ISTEX:DB3DCA224C4A8D0CF41243DB5493DA9F659E1357Le document en format XML
<record><TEI wicri:istexFullTextTei="biblStruct"><teiHeader><fileDesc><titleStmt><title xml:lang="en">Spacetime geometries and light trapping in travelling refractive index perturbations</title><author><name sortKey="Cacciatori, S L" sort="Cacciatori, S L" uniqKey="Cacciatori S" first="S L" last="Cacciatori">S L Cacciatori</name><affiliation><mods:affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</mods:affiliation></affiliation><affiliation><mods:affiliation>INFN sezione di Milano, via Celoria 16, IT-20133 Milano, Italy</mods:affiliation></affiliation><affiliation><mods:affiliation>Author to whom any correspondence should be addressed.</mods:affiliation></affiliation><affiliation><mods:affiliation>E-mail: sergio.cacciatori@uninsubria.it</mods:affiliation></affiliation></author><author><name sortKey="Belgiorno, F" sort="Belgiorno, F" uniqKey="Belgiorno F" first="F" last="Belgiorno">F. Belgiorno</name><affiliation><mods:affiliation>Dipartimento di Fisica, Universit di Milano, Via Celoria 16, IT-20133 Milano, Italy</mods:affiliation></affiliation></author><author><name sortKey="Gorini, V" sort="Gorini, V" uniqKey="Gorini V" first="V" last="Gorini">V. Gorini</name><affiliation><mods:affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</mods:affiliation></affiliation><affiliation><mods:affiliation>INFN sezione di Milano, via Celoria 16, IT-20133 Milano, Italy</mods:affiliation></affiliation></author><author><name sortKey="Ortenzi, G" sort="Ortenzi, G" uniqKey="Ortenzi G" first="G" last="Ortenzi">G. Ortenzi</name><affiliation><mods:affiliation>Dipartimento di Matematica e Applicazioni, Universitdegli Studi di Milano, Bicocca, IT-20125 Milano, Italy</mods:affiliation></affiliation></author><author><name sortKey="Rizzi, L" sort="Rizzi, L" uniqKey="Rizzi L" first="L" last="Rizzi">L. Rizzi</name><affiliation><mods:affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</mods:affiliation></affiliation></author><author><name sortKey="Sala, V G" sort="Sala, V G" uniqKey="Sala V" first="V G" last="Sala">V G Sala</name><affiliation><mods:affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</mods:affiliation></affiliation></author><author><name sortKey="Faccio, D" sort="Faccio, D" uniqKey="Faccio D" first="D" last="Faccio">D. Faccio</name><affiliation><mods:affiliation>CNISM and Department of Physics and Mathematics, Universitdell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</mods:affiliation></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">ISTEX</idno><idno type="RBID">ISTEX:DB3DCA224C4A8D0CF41243DB5493DA9F659E1357</idno><date when="2010" year="2010">2010</date><idno type="doi">10.1088/1367-2630/12/9/095021</idno><idno type="url">https://api.istex.fr/document/DB3DCA224C4A8D0CF41243DB5493DA9F659E1357/fulltext/pdf</idno><idno type="wicri:Area/Istex/Corpus">001362</idno></publicationStmt><sourceDesc><biblStruct><analytic><title level="a" type="main" xml:lang="en">Spacetime geometries and light trapping in travelling refractive index perturbations</title><author><name sortKey="Cacciatori, S L" sort="Cacciatori, S L" uniqKey="Cacciatori S" first="S L" last="Cacciatori">S L Cacciatori</name><affiliation><mods:affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</mods:affiliation></affiliation><affiliation><mods:affiliation>INFN sezione di Milano, via Celoria 16, IT-20133 Milano, Italy</mods:affiliation></affiliation><affiliation><mods:affiliation>Author to whom any correspondence should be addressed.</mods:affiliation></affiliation><affiliation><mods:affiliation>E-mail: sergio.cacciatori@uninsubria.it</mods:affiliation></affiliation></author><author><name sortKey="Belgiorno, F" sort="Belgiorno, F" uniqKey="Belgiorno F" first="F" last="Belgiorno">F. Belgiorno</name><affiliation><mods:affiliation>Dipartimento di Fisica, Universit di Milano, Via Celoria 16, IT-20133 Milano, Italy</mods:affiliation></affiliation></author><author><name sortKey="Gorini, V" sort="Gorini, V" uniqKey="Gorini V" first="V" last="Gorini">V. Gorini</name><affiliation><mods:affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</mods:affiliation></affiliation><affiliation><mods:affiliation>INFN sezione di Milano, via Celoria 16, IT-20133 Milano, Italy</mods:affiliation></affiliation></author><author><name sortKey="Ortenzi, G" sort="Ortenzi, G" uniqKey="Ortenzi G" first="G" last="Ortenzi">G. Ortenzi</name><affiliation><mods:affiliation>Dipartimento di Matematica e Applicazioni, Universitdegli Studi di Milano, Bicocca, IT-20125 Milano, Italy</mods:affiliation></affiliation></author><author><name sortKey="Rizzi, L" sort="Rizzi, L" uniqKey="Rizzi L" first="L" last="Rizzi">L. Rizzi</name><affiliation><mods:affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</mods:affiliation></affiliation></author><author><name sortKey="Sala, V G" sort="Sala, V G" uniqKey="Sala V" first="V G" last="Sala">V G Sala</name><affiliation><mods:affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</mods:affiliation></affiliation></author><author><name sortKey="Faccio, D" sort="Faccio, D" uniqKey="Faccio D" first="D" last="Faccio">D. Faccio</name><affiliation><mods:affiliation>CNISM and Department of Physics and Mathematics, Universitdell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</mods:affiliation></affiliation></author></analytic><monogr></monogr><series><title level="j">New Journal of Physics</title><title level="j" type="abbrev">New J. Phys.</title><idno type="ISSN">1367-2630</idno><idno type="eISSN">1367-2630</idno><imprint><publisher>IOP Publishing</publisher><date type="published" when="2010">2010</date><biblScope unit="volume">12</biblScope><biblScope unit="issue">9</biblScope><biblScope unit="page" from="1">1</biblScope><biblScope unit="page" to="14">14</biblScope></imprint><idno type="ISSN">1367-2630</idno></series><idno type="istex">DB3DCA224C4A8D0CF41243DB5493DA9F659E1357</idno><idno type="DOI">10.1088/1367-2630/12/9/095021</idno><idno type="articleID">355947</idno><idno type="articleNumber">095021</idno></biblStruct></sourceDesc><seriesStmt><idno type="ISSN">1367-2630</idno></seriesStmt></fileDesc><profileDesc><textClass></textClass><langUsage><language ident="en">en</language></langUsage></profileDesc></teiHeader><front><div type="abstract">In the framework of transformation optics, we show that the propagation of a locally superluminal refractive index perturbation (RIP) in a Kerr medium can be described, in the eikonal approximation, by means of a stationary metric, which we prove to be of Gordon type. Under suitable hypotheses on the RIP, we obtain a stationary but not static metric, which is characterized by an ergosphere and by a peculiar behaviour of the geodesics, which are studied numerically, also accounting for material dispersion. Finally, the equation to be satisfied by an event horizon is also displayed and briefly discussed.</div></front></TEI><istex><corpusName>iop</corpusName><author><json:item><name>S L Cacciatori</name><affiliations><json:string>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</json:string><json:string>INFN sezione di Milano, via Celoria 16, IT-20133 Milano, Italy</json:string><json:string>Author to whom any correspondence should be addressed.</json:string><json:string>E-mail: sergio.cacciatori@uninsubria.it</json:string></affiliations></json:item><json:item><name>F Belgiorno</name><affiliations><json:string>Dipartimento di Fisica, Universit di Milano, Via Celoria 16, IT-20133 Milano, Italy</json:string></affiliations></json:item><json:item><name>V Gorini</name><affiliations><json:string>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</json:string><json:string>INFN sezione di Milano, via Celoria 16, IT-20133 Milano, Italy</json:string></affiliations></json:item><json:item><name>G Ortenzi</name><affiliations><json:string>Dipartimento di Matematica e Applicazioni, Universitdegli Studi di Milano, Bicocca, IT-20125 Milano, Italy</json:string></affiliations></json:item><json:item><name>L Rizzi</name><affiliations><json:string>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</json:string></affiliations></json:item><json:item><name>V G Sala</name><affiliations><json:string>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</json:string></affiliations></json:item><json:item><name>D Faccio</name><affiliations><json:string>CNISM and Department of Physics and Mathematics, Universitdell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</json:string></affiliations></json:item></author><language><json:string>eng</json:string></language><originalGenre><json:string>paper</json:string></originalGenre><abstract>In the framework of transformation optics, we show that the propagation of a locally superluminal refractive index perturbation (RIP) in a Kerr medium can be described, in the eikonal approximation, by means of a stationary metric, which we prove to be of Gordon type. Under suitable hypotheses on the RIP, we obtain a stationary but not static metric, which is characterized by an ergosphere and by a peculiar behaviour of the geodesics, which are studied numerically, also accounting for material dispersion. Finally, the equation to be satisfied by an event horizon is also displayed and briefly discussed.</abstract><qualityIndicators><score>6.664</score><pdfVersion>1.4</pdfVersion><pdfPageSize>595.276 x 841.89 pts (A4)</pdfPageSize><refBibsNative>true</refBibsNative><keywordCount>0</keywordCount><abstractCharCount>609</abstractCharCount><pdfWordCount>5947</pdfWordCount><pdfCharCount>30180</pdfCharCount><pdfPageCount>14</pdfPageCount><abstractWordCount>97</abstractWordCount></qualityIndicators><title>Spacetime geometries and light trapping in travelling refractive index perturbations</title><genre><json:string>research-article</json:string></genre><host><volume>12</volume><publisherId><json:string>nj</json:string></publisherId><pages><total>14</total><last>14</last><first>1</first></pages><issn><json:string>1367-2630</json:string></issn><issue>9</issue><genre><json:string>journal</json:string></genre><language><json:string>unknown</json:string></language><eissn><json:string>1367-2630</json:string></eissn><title>New Journal of Physics</title></host><publicationDate>2010</publicationDate><copyrightDate>2010</copyrightDate><doi><json:string>10.1088/1367-2630/12/9/095021</json:string></doi><id>DB3DCA224C4A8D0CF41243DB5493DA9F659E1357</id><score>0.23106353</score><fulltext><json:item><original>true</original><mimetype>application/pdf</mimetype><extension>pdf</extension><uri>https://api.istex.fr/document/DB3DCA224C4A8D0CF41243DB5493DA9F659E1357/fulltext/pdf</uri></json:item><json:item><original>false</original><mimetype>application/zip</mimetype><extension>zip</extension><uri>https://api.istex.fr/document/DB3DCA224C4A8D0CF41243DB5493DA9F659E1357/fulltext/zip</uri></json:item><istex:fulltextTEI uri="https://api.istex.fr/document/DB3DCA224C4A8D0CF41243DB5493DA9F659E1357/fulltext/tei"><teiHeader><fileDesc><titleStmt><title level="a" type="main" xml:lang="en">Spacetime geometries and light trapping in travelling refractive index perturbations</title></titleStmt><publicationStmt><authority>ISTEX</authority><publisher>IOP Publishing</publisher><availability><p>IOP Publishing and Deutsche Physikalische Gesellschaft</p></availability><date>2010</date></publicationStmt><sourceDesc><biblStruct type="inbook"><analytic><title level="a" type="main" xml:lang="en">Spacetime geometries and light trapping in travelling refractive index perturbations</title><author xml:id="author-1"><persName><forename type="first">S L</forename><surname>Cacciatori</surname></persName><email>sergio.cacciatori@uninsubria.it</email><affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</affiliation><affiliation>INFN sezione di Milano, via Celoria 16, IT-20133 Milano, Italy</affiliation><affiliation>Author to whom any correspondence should be addressed.</affiliation></author><author xml:id="author-2"><persName><forename type="first">F</forename><surname>Belgiorno</surname></persName><affiliation>Dipartimento di Fisica, Universit di Milano, Via Celoria 16, IT-20133 Milano, Italy</affiliation></author><author xml:id="author-3"><persName><forename type="first">V</forename><surname>Gorini</surname></persName><affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</affiliation><affiliation>INFN sezione di Milano, via Celoria 16, IT-20133 Milano, Italy</affiliation></author><author xml:id="author-4"><persName><forename type="first">G</forename><surname>Ortenzi</surname></persName><affiliation>Dipartimento di Matematica e Applicazioni, Universitdegli Studi di Milano, Bicocca, IT-20125 Milano, Italy</affiliation></author><author xml:id="author-5"><persName><forename type="first">L</forename><surname>Rizzi</surname></persName><affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</affiliation></author><author xml:id="author-6"><persName><forename type="first">V G</forename><surname>Sala</surname></persName><affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</affiliation></author><author xml:id="author-7"><persName><forename type="first">D</forename><surname>Faccio</surname></persName><affiliation>CNISM and Department of Physics and Mathematics, Universitdell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</affiliation></author></analytic><monogr><title level="j">New Journal of Physics</title><title level="j" type="abbrev">New J. Phys.</title><idno type="pISSN">1367-2630</idno><idno type="eISSN">1367-2630</idno><imprint><publisher>IOP Publishing</publisher><date type="published" when="2010"></date><biblScope unit="volume">12</biblScope><biblScope unit="issue">9</biblScope><biblScope unit="page" from="1">1</biblScope><biblScope unit="page" to="14">14</biblScope></imprint></monogr><idno type="istex">DB3DCA224C4A8D0CF41243DB5493DA9F659E1357</idno><idno type="DOI">10.1088/1367-2630/12/9/095021</idno><idno type="articleID">355947</idno><idno type="articleNumber">095021</idno></biblStruct></sourceDesc></fileDesc><profileDesc><creation><date>2010</date></creation><langUsage><language ident="en">en</language></langUsage><abstract><p>In the framework of transformation optics, we show that the propagation of a locally superluminal refractive index perturbation (RIP) in a Kerr medium can be described, in the eikonal approximation, by means of a stationary metric, which we prove to be of Gordon type. Under suitable hypotheses on the RIP, we obtain a stationary but not static metric, which is characterized by an ergosphere and by a peculiar behaviour of the geodesics, which are studied numerically, also accounting for material dispersion. Finally, the equation to be satisfied by an event horizon is also displayed and briefly discussed.</p></abstract></profileDesc><revisionDesc><change when="2010">Published</change></revisionDesc></teiHeader></istex:fulltextTEI><json:item><original>false</original><mimetype>text/plain</mimetype><extension>txt</extension><uri>https://api.istex.fr/document/DB3DCA224C4A8D0CF41243DB5493DA9F659E1357/fulltext/txt</uri></json:item></fulltext><metadata><istex:metadataXml wicri:clean="corpus iop not found" wicri:toSee="no header"><istex:xmlDeclaration>version="1.0" encoding="ISO-8859-1"</istex:xmlDeclaration><istex:docType SYSTEM="http://ej.iop.org/dtd/iopv1_5_2.dtd" name="istex:docType"></istex:docType><istex:document><article artid="nj355947"><article-metadata><jnl-data jnlid="nj"><jnl-fullname>New Journal of Physics</jnl-fullname><jnl-abbreviation>New J. Phys.</jnl-abbreviation><jnl-shortname>NJP</jnl-shortname><jnl-issn>1367-2630</jnl-issn><jnl-coden>NJOPFM</jnl-coden><jnl-imprint>IOP Publishing</jnl-imprint><jnl-web-address>stacks.iop.org/NJP</jnl-web-address></jnl-data><volume-data><year-publication>2010</year-publication><volume-number>12</volume-number></volume-data><issue-data><issue-number>9</issue-number><coverdate>September 2010</coverdate></issue-data><article-data><article-type type="paper" sort="regular"></article-type><type-number type="paper" numbering="article" artnum="095021">095021</type-number><article-number>355947</article-number><first-page>1</first-page><last-page>14</last-page><length>14</length><pii></pii><doi>10.1088/1367-2630/12/9/095021</doi><copyright>IOP Publishing and Deutsche Physikalische Gesellschaft</copyright><ccc>1367-2630/10/095021+14$30.00</ccc><printed></printed></article-data></article-metadata><header><title-group><title>Spacetime geometries and light trapping in travelling refractive index perturbations</title><short-title>Spacetime geometries and light trapping in travelling refractive index perturbations</short-title><ej-title>Spacetime geometries and light trapping in travelling refractive index perturbations</ej-title></title-group><author-group><author address="nj355947ad1" second-address="nj355947ad2" alt-address="nj355947aad6" email="nj355947ea1"><first-names>S L</first-names><second-name>Cacciatori</second-name></author><author address="nj355947ad3"><first-names>F</first-names><second-name>Belgiorno</second-name></author><author address="nj355947ad1" second-address="nj355947ad2"><first-names>V</first-names><second-name>Gorini</second-name></author><author address="nj355947ad4"><first-names>G</first-names><second-name>Ortenzi</second-name></author><author address="nj355947ad1"><first-names>L</first-names><second-name>Rizzi</second-name></author><author address="nj355947ad1"><first-names>V G</first-names><second-name>Sala</second-name></author><author address="nj355947ad5"><first-names>D</first-names><second-name>Faccio</second-name></author><short-author-list>S L Cacciatori <italic>et al</italic></short-author-list></author-group><address-group><address id="nj355947ad1" showid="yes">Department of Physics and Mathematics, <orgname>Università dell'Insubria</orgname>, Via Valleggio 11, IT-22100 Como, <country>Italy</country></address><address id="nj355947ad2" showid="yes"><orgname>INFN sezione di Milano</orgname>, via Celoria 16, IT-20133 Milano, <country>Italy</country></address><address id="nj355947ad3" showid="yes">Dipartimento di Fisica, <orgname>Università di Milano</orgname>, Via Celoria 16, IT-20133 Milano, <country>Italy</country></address><address id="nj355947ad4" showid="yes">Dipartimento di Matematica e Applicazioni, <orgname>Università degli Studi di Milano</orgname>, Bicocca, IT-20125 Milano, <country>Italy</country></address><address id="nj355947ad5" showid="yes">CNISM and Department of Physics and Mathematics, <orgname>Università dell'Insubria</orgname>, Via Valleggio 11, IT-22100 Como, <country>Italy</country></address><address id="nj355947aad6" alt="yes" showid="yes">Author to whom any correspondence should be addressed.</address><e-address id="nj355947ea1"><email mailto="sergio.cacciatori@uninsubria.it">sergio.cacciatori@uninsubria.it</email></e-address></address-group><history received="30 April 2010" online="30 September 2010"></history><abstract-group><abstract><heading>Abstract</heading><p indent="no">In the framework of transformation optics, we show that the propagation of a locally superluminal refractive index perturbation (RIP) in a Kerr medium can be described, in the eikonal approximation, by means of a stationary metric, which we prove to be of Gordon type. Under suitable hypotheses on the RIP, we obtain a stationary but not static metric, which is characterized by an ergosphere and by a peculiar behaviour of the geodesics, which are studied numerically, also accounting for material dispersion. Finally, the equation to be satisfied by an event horizon is also displayed and briefly discussed.</p></abstract></abstract-group></header><body refstyle="numeric"><sec-level1 id="nj355947s1" label="1"><heading>Introduction</heading><p indent="no">Effective geometries for light, first introduced by Gordon [<cite linkend="nj355947bib1">1</cite>] and then extended to nonlinear electrodynamics, are able to provide analogue black hole metrics and also the possibility of performing experiments involving Hawking analogue radiation [<cite linkend="nj355947bib2">2</cite>]–[<cite linkend="nj355947bib14">14</cite>]. Recently, Philbin <italic>et al</italic> [<cite linkend="nj355947bib15">15</cite>] proposed an optical analogue in which a soliton with intensity <inline-eqn><math-text><italic>I</italic></math-text></inline-eqn>, propagating in an optical fibre, generates through the nonlinear Kerr effect a refractive index perturbation (RIP), <inline-eqn><math-text>δ<italic>n</italic>=<italic>n</italic><sub>2</sub><italic>I</italic></math-text></inline-eqn>, where <inline-eqn><math-text><italic>n</italic><sub>2</sub></math-text></inline-eqn> is the Kerr index. The same mechanism has also been generalized to a full four-dimensional (4D) geometry by Faccio <italic>et al</italic> [<cite linkend="nj355947bib16">16</cite>]. The RIP modifies the spacetime geometry as seen by co-propagating light rays and, similarly to the acoustic analogy, if the RIP is locally superluminal, i.e. if it travels faster than the phase velocity of light in the medium, a trapping horizon is formed and Hawking radiation is to be expected. A full analysis for the case of a static black hole metric has been recently discussed in [<cite linkend="nj355947bib17">17</cite>]. Our aim here is to tackle the following aspects: (i) provide a general derivation of the metric in the eikonal approximation from relativistic nonlinear electrodynamics in a medium and (ii) perform a first analysis for a generic RIP that gives rise to a stationary (but in general non-static) metric. The metric we obtain allows us to determine in a straightforward way the limits under which an ergoregion occurs. Then a study of the behaviour of null geodesics follows, with the aim of identifying conditions such that trapping of light occurs inside the dielectric perturbation. Finally, the equation to be satisfied by an event horizon is provided.</p></sec-level1><sec-level1 id="nj355947s2" label="2"><heading>Effective geometry</heading><p indent="no">In this section, we introduce the specific analogue model that we have analysed. As we need to work both in the laboratory frame and in an inertial reference frame that is moving at relativistic speed with respect to the original frame, it is convenient to adopt a covariant formalism. Then, consider a reference frame in which the dielectric medium is moving with four-velocity <inline-eqn><math-text><italic>u</italic><sup>μ</sup></math-text></inline-eqn>. Moreover, consider a medium for which the permittivity and permeability have, respectively, the form
<display-eqn id="nj355947eqn1" textype="equation" notation="LaTeX" eqnnum="1"></display-eqn><display-eqn id="nj355947eqn2" textype="equation" notation="LaTeX" eqnnum="2"></display-eqn>where <inline-eqn><math-text>η<sub>μν</sub></math-text></inline-eqn> is the Minkowski metric tensor. That is, the permittivity is a function of the invariant electric field amplitude <inline-eqn></inline-eqn> alone. In the laboratory frame, this reduces to the assumption that a polarization is induced in the medium only when it is exposed to an electric field and that the nonlinear effects depend only on the intensity of the latter. The permeability is assumed to be constant and equal to the product <inline-eqn><math-text>μ</math-text></inline-eqn> of the vacuum permeability <inline-eqn><math-text>μ<sub>0</sub></math-text></inline-eqn> and the relative permeability <inline-eqn><math-text>μ<sub><upright>r</upright></sub></math-text></inline-eqn> of the medium. These conditions indeed hold true for diverse nonlinear media. For example, for the case of an isotropic <inline-eqn><math-text>χ<sup>(3)</sup></math-text></inline-eqn> medium, in the laboratory frame
<display-eqn id="nj355947eqn3" textype="equation" notation="LaTeX" eqnnum="3"></display-eqn>where <inline-eqn><math-text>ϵ<sub>0</sub></math-text></inline-eqn> is the vacuum permittivity.</p><p>Nonlinear effects are usually very small and a very strong field is needed to exploit the nonlinearity of the medium itself. Thus, we assume that there exists a strong, background field <inline-eqn><math-text><italic>E</italic></math-text></inline-eqn>, which fixes the properties of the dielectric. Note that these assumptions are rather general. In fact, equations (<eqnref linkend="nj355947eqn1">1</eqnref>) and (<eqnref linkend="nj355947eqn2">2</eqnref>) hold for any nonlinear, homogeneous dielectric, no matter how the strong background pulse is generated [<cite linkend="nj355947bib6">6</cite>]. We are interested in the small and rapidly oscillating fluctuations of the electromagnetic field on top of such a background<fnref linkend="nj355947fn7"><sup>7</sup></fnref>.</p><p>According to the results of [<cite linkend="nj355947bib6">6</cite>, <cite linkend="nj355947bib11">11</cite>, <cite linkend="nj355947bib12">12</cite>], these small fluctuations ‘feel’ an effective curved spacetime. The result is what a relativist would call a ‘bi-metric’ theory. In other words, there exist two classes of rays, namely ordinary and extraordinary rays, which correspond to two different polarizations of the background field. This phenomenon is analogous to birefringence in uni-axial crystals, where the ordinary and extraordinary rays each obey distinct quadratic dispersion relations. In the analogue gravity terminology, each ray ‘feels’ a different geometry, or <italic>analogue metric</italic>: the Gordon's metric
<display-eqn id="nj355947eqn4" textype="equation" notation="LaTeX" eqnnum="4"></display-eqn>for the ordinary rays and the metric
<display-eqn id="nj355947eqn5" textype="equation" notation="LaTeX" eqnnum="5"></display-eqn>for the extraordinary rays. In (<eqnref linkend="nj355947eqn4">4</eqnref>)–(<eqnref linkend="nj355947eqn5">5</eqnref>), we have introduced the ‘nonlinear parameter’ <inline-eqn></inline-eqn>, the background field polarization four-vector <inline-eqn></inline-eqn>, and the refractive index <inline-eqn></inline-eqn>. In the Cartesian coordinates <inline-eqn><math-text>(<italic>t</italic>, <italic>x</italic>, <italic>y</italic>, <italic>z</italic>)</math-text></inline-eqn>, <inline-eqn><math-text>η<sub>μν</sub>=<upright>diag</upright>(<italic>c</italic><sup>2</sup>, −1, −1, −1)</math-text></inline-eqn>.</p><p>In the <secref linkend="nj355947s7">appendix</secref>, we show that, under reasonable approximations, the extraordinary metric is equivalent to the Gordon one, albeit with a different effective Kerr perturbation of the refractive index, i.e. <inline-eqn><math-text><italic>n</italic>=<italic>n</italic><sub>0</sub>+3<italic>dn</italic></math-text></inline-eqn> instead of <inline-eqn><math-text><italic>n</italic>=<italic>n</italic><sub>0</sub>+<italic>dn</italic></math-text></inline-eqn>. Thus, we restrict our attention to the Gordon metric.</p><p>In the following, we study the light-like geodesics of the analogue metrics. It is important to note that these geodesics describe the spacetime trajectories of monochromatic waves. In the present paper, we often refer to these trajectories with the term ‘light rays’. However, these should not be confused with the usual light rays referred to in geometrical optics which, in the Newtonian picture, describe the trajectories of wave packets (or ‘photons’). For this reason, in the forthcoming sections, we analyse the monochromatic components (congruences of light-like geodesics) of the electromagnetic field and the analogue geodesic motion of the related constant phase surfaces, which are precisely the objects involved in the study of the analogue Hawking radiation [<cite linkend="nj355947bib18">18</cite>, <cite linkend="nj355947bib19">19</cite>]. By contrast, we do not consider the notion of the ‘photon’ (or wave packet) here.</p></sec-level1><sec-level1 id="nj355947s3" label="3"><heading>The metric in the refractive index perturbation (RIP) frame</heading><p indent="no">Consider the case of a pulse moving with uniform velocity <inline-eqn><math-text><italic>v</italic></math-text></inline-eqn> along the positive <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn>-axis. According to the above discussion, the pulse excites a localized RIP which also moves accordingly. The latter can then be analytically described as a function of <inline-eqn><math-text><italic>x</italic>−<italic>vt</italic></math-text></inline-eqn> and, assuming cylindrical symmetry, of the transversal coordinate <inline-eqn></inline-eqn>. Thus, in the laboratory frame, <inline-eqn><math-text><italic>n</italic>=<italic>n</italic><sub>0</sub>+δ<italic>n</italic>(<italic>x</italic>−<italic>vt</italic>, ρ)</math-text></inline-eqn>. In this situation, it is convenient to express the Gordon metric in a reference frame which is co-moving with the RIP itself:
<display-eqn id="nj355947eqn6" textype="equation" notation="LaTeX" eqnnum="6" lines="multiline"></display-eqn>where the primed coordinates are relative to the co-moving frame, and
<display-eqn id="nj355947eqn7" textype="equation" notation="LaTeX" eqnnum="7"></display-eqn>Here, <inline-eqn><math-text>δ<italic>n</italic></math-text></inline-eqn> is a <inline-eqn><math-text><italic>C</italic><sup>∞</sup></math-text></inline-eqn> function, rapidly decaying at infinity and with a single maximum of height <inline-eqn><math-text>η</math-text></inline-eqn>, describing the RIP. It is then evident that both <inline-eqn><math-text>∂<sub><italic>t</italic>′</sub></math-text></inline-eqn> and <inline-eqn><math-text>∂<sub>&phis;</sub></math-text></inline-eqn> are Killing vectors for the given class of metrics. Any specific choice for the function <inline-eqn><math-text>δ<italic>n</italic></math-text></inline-eqn> gives rise to a specific metric in the aforementioned class. We point out that an isotropic refractive index in the laboratory frame corresponds to an anisotropic refractive index in the pulse frame due to length contraction associated with a boost in the <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn>-direction. If the refractive index depends explicitly on <inline-eqn><math-text>ρ</math-text></inline-eqn>, the metric is stationary but not static: the integrability conditions stated in the Frobenius theorem [<cite linkend="nj355947bib20">20</cite>] are not satisfied. In particular, if <inline-eqn><math-text>ξ<sub>μ</sub>=(<italic>g</italic><sub>00</sub>, <italic>g</italic><sub>01</sub>, 0, 0)</math-text></inline-eqn> are the covariant components of the time-like Killing vector <inline-eqn><math-text>∂<sub><italic>t</italic>′</sub></math-text></inline-eqn>, the integrability conditions <inline-eqn><math-text>ξ<sub>[λ</sub>∇<sub>μ</sub>ξ<sub>ν]</sub>=0</math-text></inline-eqn>, for <inline-eqn><math-text>λ=0, μ=1</math-text></inline-eqn> give <inline-eqn><math-text>ξ<sub>0</sub>(∇<sub>1</sub>ξ<sub>ν</sub>−∇<sub>ν</sub>ξ<sub>1</sub>)+ξ<sub>1</sub>(∇<sub>ν</sub>ξ<sub>0</sub>−∇<sub>0</sub>ξ<sub>ν</sub>)=0</math-text></inline-eqn>, i.e. for <inline-eqn><math-text>ν=ρ</math-text></inline-eqn>, <inline-eqn><math-text><italic>g</italic><sub>00</sub>(−∂<sub>ρ</sub><italic>g</italic><sub>01</sub>)+<italic>g</italic><sub>01</sub>∂<sub>ρ</sub><italic>g</italic><sub>00</sub>=0</math-text></inline-eqn>, which is easily seen to be not satisfied. The static case is obtained if <inline-eqn><math-text>δ<italic>n</italic></math-text></inline-eqn> does not depend on <inline-eqn><math-text>ρ</math-text></inline-eqn>, and it has been extensively studied in [<cite linkend="nj355947bib17">17</cite>]. Here, we consider the more general stationary, non-static case<fnref linkend="nj355947fn8"><sup>8</sup></fnref>. We first point out that for our class of metrics an ergosurface is allowed, given by the locus of points satisfying the equation
<display-eqn id="nj355947eqn8" textype="equation" notation="LaTeX" eqnnum="8"></display-eqn>This equation admits solutions only for the range of RIP velocities <inline-eqn><math-text><italic>v</italic></math-text></inline-eqn> given by
<display-eqn id="nj355947eqn9" textype="equation" notation="LaTeX" eqnnum="9"></display-eqn>At the ergosurface, which is the boundary of the ergoregion, the Killing vector <inline-eqn><math-text>∂<sub><italic>t</italic>′</sub></math-text></inline-eqn> vanishes, and on entering the ergoregion it changes sign. Here, in contrast with the static case, the ergosurface does not coincide with the event horizon. Indeed, the given hypersurface is not a light-like one, as can be easily verified by looking at the behaviour of its normal vector. However, for <inline-eqn><math-text>ρ=0</math-text></inline-eqn> it actually becomes light-like, i.e. it meets the event horizon in the antipodal points which solve
<display-eqn id="nj355947eqn10" textype="equation" notation="LaTeX" eqnnum="10"></display-eqn>For stationary non-static situations, it is more difficult to identify a non-trivial event horizon in effective geometries (see e.g. Visser's analysis for the case of acoustic black holes [<cite linkend="nj355947bib21">21</cite>]). Before discussing the question of the existence of an event horizon, we investigate whether the path of light rays, that is, the null geodesic curves of metric (<eqnref linkend="nj355947eqn6">6</eqnref>), do suggest the existence of such a horizon.</p></sec-level1><sec-level1 id="nj355947s4" label="4"><heading>Geodesics and trapping of light</heading><p indent="no">In the laboratory frame, the metric (<eqnref linkend="nj355947eqn6">6</eqnref>) is
<display-eqn id="nj355947eqn11" textype="equation" notation="LaTeX" eqnnum="11"></display-eqn>We model the RIP with a Gaussian-shaped function that propagates along the <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn>-axis with velocity <inline-eqn><math-text><italic>v</italic></math-text></inline-eqn>. We also suppose that, during the propagation, it keeps its shape both in space and in time. Assuming cylindrical symmetry, we need to account only for one of the two transverse coordinates, say, <inline-eqn><math-text><italic>y</italic></math-text></inline-eqn>. Then, we write the refractive index in the form
<display-eqn id="nj355947eqn12" textype="equation" notation="LaTeX" eqnnum="12" lines="multiline"></display-eqn>Thus, the geodesic equations for null rays take the explicit form
<display-eqn id="nj355947eqn13" textype="equation" notation="LaTeX" eqnnum="13"></display-eqn><display-eqn id="nj355947eqn14" textype="equation" notation="LaTeX" eqnnum="14"></display-eqn><display-eqn id="nj355947eqn15" textype="equation" notation="LaTeX" eqnnum="15"></display-eqn>where the dot denotes derivation w.r.t. an affine parameter <inline-eqn><math-text>λ</math-text></inline-eqn>. A trivial solution is <inline-eqn></inline-eqn>, <inline-eqn><math-text><italic>y</italic>(<italic>t</italic>)=0</math-text></inline-eqn>. However, in the general case a numerical analysis is required. In the following subsection, we solve the equations above using Matlab's ode45 function. This analysis will provide an insight into the underlying physics.</p><sec-level2 id="nj355947s4.1" label="4.1"><heading>Numerical solutions and trapping of light</heading><p indent="no">Along a geodesic we have <inline-eqn></inline-eqn>, and taking the first derivative (with respect to the affine parameter <inline-eqn><math-text>λ</math-text></inline-eqn>) of the relation <inline-eqn></inline-eqn>, we find an equation for <inline-eqn></inline-eqn> that determines the evolution of the frequency of the ray along a geodesic. At this point a comment is in order: here we have defined <inline-eqn><math-text>ω</math-text></inline-eqn> as the zeroth covariant component of the four-momentum <inline-eqn><math-text><italic>k</italic><sup>μ</sup></math-text></inline-eqn>. It is easy to see that this corresponds to measuring frequencies as the inverse of the lapse of coordinate time between two successive crests (this time is the Minkowskian time measured by the laboratory observer). This is of course different from the analogue proper time <inline-eqn><math-text><upright>d</upright>τ=<upright>d</upright><italic>s</italic>/<italic>c</italic></math-text></inline-eqn> of the analogue metric (<eqnref linkend="nj355947eqn11">11</eqnref>), which has no physical significance: experimental measurements refer to the true Minkowski metric and not to the analogue curved metric.</p><p>Then, for the evolution of the frequency <inline-eqn><math-text>ω</math-text></inline-eqn> along the path of the ray, we obtain
<display-eqn id="nj355947eqn16" textype="equation" notation="LaTeX" eqnnum="16"></display-eqn>For the sake of definiteness assume the propagation to be in fused silica and the rays to have initial wavelength <inline-eqn><math-text>λ<sub><upright>in</upright></sub>=527 <upright>nm</upright></math-text></inline-eqn>. We assume the RIP to have <inline-eqn><math-text>σ=1 <upright>m</upright></math-text></inline-eqn> and amplitude <inline-eqn><math-text>η=0.1</math-text></inline-eqn>, which is quite an unrealistic situation, but simple to deal with. Indeed, the aim of the simulation is merely to study the qualitative behaviour of the geodesics. We have verified that this behaviour does not change, even for much smaller and realistic values of <inline-eqn><math-text>σ∼1</math-text></inline-eqn>–<inline-eqn><math-text>100 μ<upright>m</upright></math-text></inline-eqn>.</p><p>The other parameters of the simulation are the velocity <inline-eqn><math-text><italic>v</italic></math-text></inline-eqn> of the RIP and the initial conditions of the rays, i.e. the array <inline-eqn></inline-eqn>. In the simulation, the rays start away from the RIP (which is right moving) and their initial velocity is directed along the <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn>-axis, as shown in the sketch of figure <figref linkend="nj355947fig1">1</figref>(a). More precisely, the initial conditions are <inline-eqn></inline-eqn>, with the additional constraints <inline-eqn></inline-eqn> and <inline-eqn></inline-eqn>.</p><figure id="nj355947fig1" parts="single" width="column" position="float" pageposition="top" printstyle="normal" orientation="port"><graphic position="indented"><graphic-file version="print" format="EPS" width="26.0pc" printcolour="no" filename="images/nj355947fig1.eps"></graphic-file><graphic-file version="ej" format="JPEG" printcolour="no" filename="images/nj355947fig1.jpg"></graphic-file></graphic><caption type="figure" id="nj355947fc1" label="Figure 1"><p indent="no">(a) Sketch showing the initial conditions of the simulation. (b) Light ray phase velocity (<inline-eqn><math-text><italic>c</italic>/<italic>n</italic></math-text></inline-eqn>). The shaded region highlights the range of velocities of the RIP for which the rays are trapped. The two straight lines indicate the values of <inline-eqn><math-text><italic>v</italic>/<italic>c</italic></math-text></inline-eqn> used in the simulations.</p></caption></figure><p>The results of the simulation show that the behaviour of null geodesics is strongly dependent on the velocity <inline-eqn><math-text><italic>v</italic></math-text></inline-eqn> of the RIP. Namely, there exists a range of velocities for which the rays are trapped<fnref linkend="nj355947fn9"><sup>9</sup></fnref>, and these coincide with those predicted by equation (<eqnref linkend="nj355947eqn9">9</eqnref>) (shaded region in figure <figref linkend="nj355947fig1">1</figref>(b)). On the other hand, for lower values of <inline-eqn><math-text><italic>v</italic></math-text></inline-eqn>, the rays that enter the RIP are bent and eventually escape. In figure <figref linkend="nj355947fig1">1</figref>(b), the black curve represents the velocity of an axial ray, i.e. of a ray propagating along the <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn>-axis, versus <inline-eqn><math-text><italic>x</italic>−<italic>vt</italic></math-text></inline-eqn>. In the same figure, we represent with horizontal lines two choices for the RIP velocity, representing two possible situations: a faster (green line) one and a slower (red line) one. In the first case (faster RIP), the green line intersects the black curve. As the ray approaches the point corresponding to the first intersection, it slows down due to the increasing refractive index, until it acquires the same velocity of the RIP. Asymptotically, the ray is ‘trapped’ and moves together with the RIP. In the other case (slower RIP), the ray velocity is always larger than the RIP's velocity (red line of figure <figref linkend="nj355947fig1">1</figref>(b)). Then, the ray can cross the RIP and escape. In this ‘subluminal’ case, no trapping occurs.</p><sec-level3 id="nj355947s4.1.1" label="4.1.1"><heading>Subluminal RIP</heading><p indent="no">The RIP propagates at <inline-eqn><math-text><italic>v</italic>=0.56<italic>c</italic></math-text></inline-eqn> (red line in figure <figref linkend="nj355947fig1">1</figref>(b)), and it is always slower than the light rays. In figure <figref linkend="nj355947fig2">2</figref>(a) we show the numerical results of the simulation in the co-moving reference frame. The rays start on the left of the RIP. They propagate along straight lines with velocity <inline-eqn><math-text><italic>c</italic>/<italic>n</italic><sub>0</sub></math-text></inline-eqn> until they reach the RIP. Then, inside the RIP, they slow down, bend and finally escape from the RIP, propagating along straight lines in different directions, depending on their impact parameter <inline-eqn><math-text><italic>y</italic><sub><upright>in</upright></sub></math-text></inline-eqn>. The projection of the trajectory on the <inline-eqn><math-text>(<italic>t</italic>′ ,<italic>x</italic>′ )</math-text></inline-eqn> plane, which is also shown figure <figref linkend="nj355947fig2">2</figref>(a), evidently shows a change of slope corresponding to the trajectory being inside the RIP. This means that the rays, while crossing the RIP, slow down. Figure <figref linkend="nj355947fig2">2</figref>(b) represents the evolution of <inline-eqn><math-text>ω</math-text></inline-eqn> along the <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn>-direction in the laboratory reference frame. We observe that, as expected, the frequency varies consistently when the rays reach the RIP. The rays go through a blueshift inside the RIP and then they emerge with a Doppler redshift, due to the fact that the rays are bent and thus emerge at different angles. This behaviour is described in more detail in [<cite linkend="nj355947bib16">16</cite>].</p><figure id="nj355947fig2" parts="single" width="column" position="float" pageposition="top" printstyle="normal" orientation="port"><graphic position="indented"><graphic-file version="print" format="EPS" width="25.6pc" printcolour="no" filename="images/nj355947fig2.eps"></graphic-file><graphic-file version="ej" format="JPEG" printcolour="no" filename="images/nj355947fig2.jpg"></graphic-file></graphic><caption type="figure" id="nj355947fc2" label="Figure 2"><p indent="no">Numerical results in the co-moving reference frame (a) and frequency evolution in the laboratory reference frame (b) for the subluminal case.</p></caption></figure></sec-level3><sec-level3 id="nj355947s4.1.2" label="4.1.2"><heading>Trapping RIP</heading><p indent="no">The RIP propagates at <inline-eqn><math-text><italic>v</italic>=0.66<italic>c</italic></math-text></inline-eqn> (green line in figure <figref linkend="nj355947fig1">1</figref>(b)). In figure <figref linkend="nj355947fig3">3</figref>(a) we show the numerical results in the co-moving reference frame. The rays propagate initially along straight lines with velocity <inline-eqn><math-text><italic>c</italic>/<italic>n</italic><sub>0</sub></math-text></inline-eqn> (for clarity, we recall that rays propagate along geodesics with the <italic>phase</italic> velocity of the corresponding electromagnetic field). When they reach the RIP, they are bent and at the same time slow down. As shown in figure <figref linkend="nj355947fig1">1</figref>(b), at a certain point, away from the centre of the RIP, the rays will have exactly the same velocity as the RIP, and can no longer escape from it. In the co-moving reference frame, they continuously slow down and at the same time they oscillate in the transverse direction with decreasing amplitude oscillations, until, asymptotically, all of them stop at the same point. This point, represented in red in figure <figref linkend="nj355947fig1">1</figref>(a), is the white hole horizon <inline-eqn><math-text><italic>x</italic><sub>−</sub></math-text></inline-eqn> (see the next section). If we look at the slope of the projection <inline-eqn><math-text><italic>t</italic>′ =<italic>t</italic>′ (<italic>x</italic>′ )</math-text></inline-eqn> we see that it diverges, which means that the rays velocity tends to zero corresponding to this point. Clearly, in the laboratory reference frame the rays do not stop, but they move together with the RIP. As shown in figure <figref linkend="nj355947fig3">3</figref>(b), as the rays reach the RIP, their frequency starts growing indefinitely, so we have an infinite blueshift.</p><figure id="nj355947fig3" parts="single" width="column" position="float" pageposition="top" printstyle="normal" orientation="port"><graphic position="indented"><graphic-file version="print" format="EPS" width="25.8pc" printcolour="no" filename="images/nj355947fig3.eps"></graphic-file><graphic-file version="ej" format="JPEG" printcolour="no" filename="images/nj355947fig3.jpg"></graphic-file></graphic><caption type="figure" id="nj355947fc3" label="Figure 3"><p indent="no">Numerical results in the co-moving reference frame (a) and frequency evolution in the laboratory reference frame (b) for the trapping case.</p></caption></figure></sec-level3></sec-level2><sec-level2 id="nj355947s4.2" label="4.2"><heading>Effects of dispersion</heading><p indent="no">In the above considerations, we did not account for the material dispersion, that is, we assumed <inline-eqn><math-text><italic>n</italic>(<bold>r</bold>, <italic>t</italic>)</math-text></inline-eqn> independent of <inline-eqn><math-text>ω</math-text></inline-eqn>. However, the refractive index of a medium is in general a complicated function of the frequency <inline-eqn><math-text>ω</math-text></inline-eqn> of the incident wave, that is, <inline-eqn><math-text><italic>n</italic>=<italic>n</italic>(<bold>r</bold>,<italic>t</italic>, ω)</math-text></inline-eqn>. We restrict our attention to the case of ‘normal dispersion’, that is, we assume that the refraction index grows with frequency. In this case, the full dependence of the refractive index from <inline-eqn><math-text>ω</math-text></inline-eqn> is phenomenologically described by the Sellmeier equation. A standard way of treating dispersion is to expand <inline-eqn><math-text><italic>n</italic>(ω)</math-text></inline-eqn> in power series of <inline-eqn><math-text>ω</math-text></inline-eqn> (about an initial frequency <inline-eqn><math-text>ω<sub>0</sub></math-text></inline-eqn>), or equivalently, by the relation <inline-eqn><math-text><italic>ck</italic>=ω<italic>n</italic></math-text></inline-eqn>, to expand <inline-eqn><math-text><italic>k</italic></math-text></inline-eqn> in power series of <inline-eqn><math-text>ω</math-text></inline-eqn> as follows:
<display-eqn id="nj355947eqn17" textype="equation" notation="LaTeX" eqnnum="17"></display-eqn>The material dispersion also affects the metric, which therefore depends on the frequency:
<display-eqn id="nj355947eqn18" textype="equation" notation="LaTeX" eqnnum="18"></display-eqn>Therefore, different rays propagating in the medium ‘feel’ different metrics, depending on their frequency. Moreover, even a single ray will see an evolving metric while travelling, since its frequency (and hence the effective metric) changes during the propagation. The geodesic equations (<eqnref linkend="nj355947eqn13">13</eqnref>)–(<eqnref linkend="nj355947eqn15">15</eqnref>) keep the same form as in the non-dispersive case, while equation (<eqnref linkend="nj355947eqn16">16</eqnref>) now takes the form
<display-eqn id="nj355947eqn19" textype="equation" notation="LaTeX" eqnnum="19"></display-eqn>We highlight that, contrary to the non-dispersive case, the evolution of the geodesics depends on the change in frequency of the ray, owing to the explicit dependence on <inline-eqn><math-text>ω</math-text></inline-eqn> in the geodesic equations. In the case of a subluminal RIP, dispersion does not modify the qualitative behaviour of the light ray, which remains very similar to the dispersionless case shown in figure <figref linkend="nj355947fig2">2</figref>. On the other hand, for a trapping RIP, the presence of dispersion modifies substantially the geodesics. The results of the numerical simulation are shown in figure <figref linkend="nj355947fig4">4</figref>. The trapping feature of the horizon is still evident. However, taking dispersion into account, the rays are no longer indefinitely trapped inside the RIP and eventually escape. This behaviour is clearly depicted in figure <figref linkend="nj355947fig4">4</figref>(a), where the results are shown in the co-moving reference frame. Note that, at the beginning, the (phase) velocity of the rays is positive, which means that the rays are faster than the RIP, but then it becomes negative, implying that the rays have become slower than the RIP. Let us clarify this aspect. The phase velocity of a ray is defined as <inline-eqn><math-text><italic>v</italic><sub><upright>r</upright></sub>=ω/<italic>k</italic></math-text></inline-eqn>. Under conditions of normal dispersion, i.e. <inline-eqn></inline-eqn>, as in our case, <inline-eqn><math-text><italic>k</italic></math-text></inline-eqn> grows faster than linearly with <inline-eqn><math-text>ω</math-text></inline-eqn> (while without dispersion it would have a linear growth). When the rays approach the RIP, their frequency grows and, at the same time, their phase velocity decreases. Thus, the rays will begin to lag behind with respect to the RIP. When they are out of the RIP they will keep their final frequency. In figure <figref linkend="nj355947fig4">4</figref>(b), we show that the frequency grows only by a finite amount. Therefore, the rays leave the RIP with a finite (albeit very large) blueshift.</p><figure id="nj355947fig4" parts="single" width="column" position="float" pageposition="top" printstyle="normal" orientation="port"><graphic position="indented"><graphic-file version="print" format="EPS" width="26.0pc" printcolour="no" filename="images/nj355947fig4.eps"></graphic-file><graphic-file version="ej" format="JPEG" printcolour="no" filename="images/nj355947fig4.jpg"></graphic-file></graphic><caption type="figure" id="nj355947fc4" label="Figure 4"><p indent="no">Numerical results in the co-moving reference frame (a) and frequency evolution in the laboratory reference frame (b) for the trapping case when we account for dispersion.</p></caption></figure><sec-level3 id="nj355947s4.2.1" label="4.2.1"><heading>Trapping time</heading><p indent="no">By increasing the velocity of the RIP from <inline-eqn><math-text><italic>v</italic>=0.5<italic>c</italic></math-text></inline-eqn> to <inline-eqn><math-text><italic>v</italic>=0.65<italic>c</italic></math-text></inline-eqn> (i.e. approaching the trapping case), the rays spend more and more time inside the RIP. In figure <figref linkend="nj355947fig5">5</figref>(a)–(g), we show this behaviour for different velocities, in the dispersionless case. In figure <figref linkend="nj355947fig5">5</figref>(h), we have plotted the interaction time versus the velocity of the RIP in the dispersion (red curve) and dispersionless (blue curve) case. Neglecting dispersion, the interaction time tends to infinity as the velocity of the RIP approaches the trapping case. We can interpret this divergence as a (meta)stable state where the photon is trapped inside the pulse. An analytical estimate for this trapping time is given, for the axial geodesic, by the formula <inline-eqn></inline-eqn>, where <inline-eqn><math-text><italic>r</italic></math-text></inline-eqn> represents the longitudinal extension of the RIP. In the dispersionless case <inline-eqn><math-text><italic>t</italic><sub><upright>trap</upright></sub></math-text></inline-eqn> diverges, whereas in the dispersive case it is always finite.</p><figure id="nj355947fig5" parts="single" width="column" position="float" pageposition="top" printstyle="normal" orientation="port"><graphic position="left"><graphic-file version="print" format="EPS" width="37.8pc" printcolour="no" filename="images/nj355947fig5.eps"></graphic-file><graphic-file version="ej" format="JPEG" printcolour="no" filename="images/nj355947fig5.jpg"></graphic-file></graphic><caption type="figure" id="nj355947fc5" label="Figure 5"><p indent="no">(a)–(g) The trapping dynamics in the absence of dispersion for different RIP velocities: (a) <inline-eqn><math-text><italic>v</italic>=0.5</math-text></inline-eqn>, (b) <inline-eqn><math-text><italic>v</italic>=0.56</math-text></inline-eqn>, (c) <inline-eqn><math-text><italic>v</italic>=0.6</math-text></inline-eqn>, (d) <inline-eqn><math-text><italic>v</italic>=0.62</math-text></inline-eqn>, (e) <inline-eqn><math-text><italic>v</italic>=0.63</math-text></inline-eqn>, (f) <inline-eqn><math-text><italic>v</italic>=0.64</math-text></inline-eqn> and (g) <inline-eqn><math-text><italic>v</italic>=0.65</math-text></inline-eqn>. (h) The interaction time versus the velocity of the RIP in the case of no dispersion (blue curve) and in the case of dispersion (red curve).</p></caption></figure></sec-level3></sec-level2></sec-level1><sec-level1 id="nj355947s5" label="5"><heading>Horizon in optics</heading><p indent="no">In this section, we provide a definition for the horizons from the optical point of view, in the dispersionless case. For the sake of definiteness, we focus on the ‘white hole’ horizon, but similar considerations hold true for the ‘black hole’ horizon. We work in the laboratory frame, where the pulse is moving in the <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn>-direction, with velocity <inline-eqn><math-text><italic>v</italic></math-text></inline-eqn>. The situation is depicted in figure <figref linkend="nj355947fig6">6</figref>(a).</p><figure id="nj355947fig6" parts="single" width="column" position="float" pageposition="top" printstyle="normal" orientation="port"><graphic position="indented"><graphic-file version="print" format="EPS" width="25.9pc" printcolour="no" filename="images/nj355947fig6.eps"></graphic-file><graphic-file version="ej" format="JPEG" printcolour="no" filename="images/nj355947fig6.jpg"></graphic-file></graphic><caption type="figure" id="nj355947fc6" label="Figure 6"><p indent="no">(a) Picture of the ergosphere and of the hypothetical horizon at a fixed point of time in the laboratory frame. The small segments at <inline-eqn><math-text><italic>P</italic></math-text></inline-eqn> and <inline-eqn><math-text><italic>Q</italic></math-text></inline-eqn> represent infinitesimal portions of null geodesics (light rays) passing, respectively, through <inline-eqn><math-text><italic>P</italic></math-text></inline-eqn> and <inline-eqn><math-text><italic>Q</italic></math-text></inline-eqn> at that given point of time. (b) Normal velocity <italic>u</italic> compared to the pulse velocity, <italic>v</italic>.</p></caption></figure><p>Recall that the ergosphere is defined as the set of points where the velocity of light is the same as the velocity of the pulse, which means that <inline-eqn><math-text><italic>n</italic>(<italic>X</italic>)=<italic>c</italic>/<italic>v</italic></math-text></inline-eqn> for all points <inline-eqn><math-text><italic>X</italic></math-text></inline-eqn> on the ergosphere. In particular, it follows that a light ray reaching the ergosphere, for example in <inline-eqn><math-text><italic>P</italic></math-text></inline-eqn>, will not enter the ball if it is moving in the <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn>-direction. However, if it is directed along the NW–SE direction (for example toward the centre), then it will cross over the ergosphere, thus entering the ball. This is because it ‘sees’ the surface move more slowly in such a direction: indeed, note that a plane forming an angle <inline-eqn><math-text>&thetas;</math-text></inline-eqn> w.r.t. <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn> and moving with velocity <inline-eqn><math-text><italic>v</italic></math-text></inline-eqn> in the <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn>-direction appears to move with velocity <inline-eqn><math-text><italic>u</italic>=<italic>v</italic> <upright>sin</upright> &thetas;<<italic>c</italic>/<italic>n</italic>(<italic>P</italic>)</math-text></inline-eqn> along its normal direction, see figure <figref linkend="nj355947fig6">6</figref>(b).</p><p>Therefore, only a single point on the ergosphere, namely the axial point <inline-eqn><math-text><italic>x</italic><sub>−</sub></math-text></inline-eqn>, behaves as a true white hole horizon. In order to determine whether a full event horizon extends beyond <inline-eqn><math-text><italic>x</italic><sub>−</sub></math-text></inline-eqn>, we impose the condition <inline-eqn><math-text><italic>u</italic>=<italic>c</italic>/<italic>n</italic>(<italic>Q</italic>)</math-text></inline-eqn> on some point internal to the ergosphere (where <inline-eqn><math-text><italic>n</italic>(<italic>Q</italic>)><italic>n</italic>(<italic>P</italic>)=<italic>c</italic>/<italic>v</italic></math-text></inline-eqn>), which then gives
<display-eqn id="nj355947eqn20" textype="equation" notation="LaTeX" eqnnum="20"></display-eqn>Now, suppose that the horizon is described by a surface <inline-eqn><math-text>ρ=<italic>h</italic>(<italic>x</italic>−<italic>vt</italic>)</math-text></inline-eqn>, where <inline-eqn><math-text>ρ</math-text></inline-eqn> is the radial coordinate. Then
<display-eqn id="nj355947eqn21" textype="equation" notation="LaTeX" eqnnum="21"></display-eqn>Note that equations (<eqnref linkend="nj355947eqn20">20</eqnref>) and (<eqnref linkend="nj355947eqn21">21</eqnref>) are indeed satisfied by the axial point <inline-eqn><math-text><italic>x</italic><sub>−</sub></math-text></inline-eqn> (for which equation (<eqnref linkend="nj355947eqn21">21</eqnref>) gives <inline-eqn><math-text>&thetas;=90°</math-text></inline-eqn>). A simple manipulation gives
<display-eqn id="nj355947eqn22" textype="equation" notation="LaTeX" eqnnum="22"></display-eqn>Finally, in the comoving frame, using the coordinates <inline-eqn><math-text>(<italic>x</italic>′ , ρ)</math-text></inline-eqn>, the horizon equation is
<display-eqn id="nj355947eqn23" textype="equation" notation="LaTeX" eqnnum="23"></display-eqn>The same equation can be obtained in the following way. It is a common strategy to parameterize the horizon by means of a level set function <inline-eqn><math-text><italic>F</italic></math-text></inline-eqn> such that <inline-eqn><math-text><italic>F</italic>=0</math-text></inline-eqn> on the horizon; the light-like condition for the event horizon is <inline-eqn><math-text><italic>g</italic><sup>μν</sup>(∂<sub>μ</sub><italic>F</italic>)(∂<sub>ν</sub><italic>F</italic>)=0</math-text></inline-eqn>. In our case, there are two further requirements: axial symmetry and that the horizon is stationary (i.e. <inline-eqn><math-text>∂<sub><italic>t</italic>′</sub><italic>F</italic></math-text></inline-eqn>=0). Then, in the RIP frame the condition that imposes the normal to the event horizon to be light-like is
<display-eqn id="nj355947eqn24" textype="equation" notation="LaTeX" eqnnum="24"></display-eqn>In the static case, <inline-eqn><math-text>∂<sub>ρ</sub><italic>F</italic>=0</math-text></inline-eqn>, and then from (<eqnref linkend="nj355947eqn24">24</eqnref>) we deduce that the horizon consists of two infinite extended planes, orthogonal to the <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn>-axis, located at the solutions <inline-eqn><math-text><italic>x</italic>′ <sub>±</sub></math-text></inline-eqn> of <inline-eqn></inline-eqn>, as described in [<cite linkend="nj355947bib17">17</cite>]. In the non-static case, it is evident that the horizon has to be contained within the ergosurface, because it is necessary that <inline-eqn><math-text>1−<italic>n</italic><sup>2</sup>(<italic>x</italic>, ρ)v<sup>2</sup>/<italic>c</italic><sup>2</sup>⩽0</math-text></inline-eqn>. Moreover (at least locally, due to the implicit function theorem), we assume that <inline-eqn></inline-eqn>, i.e. <inline-eqn><math-text><italic>F</italic>(<italic>x</italic>′ , ρ)=ρ−<italic>h</italic>(<italic>x</italic>′ )</math-text></inline-eqn>. Then, equation (<eqnref linkend="nj355947eqn24">24</eqnref>) becomes <inline-eqn></inline-eqn>. As a consequence, we obtain <inline-eqn></inline-eqn>, which explicitly corresponds to (<eqnref linkend="nj355947eqn23">23</eqnref>).</p><p>According to equations (<eqnref linkend="nj355947eqn20">20</eqnref>) and (<eqnref linkend="nj355947eqn21">21</eqnref>), as described above, the points <inline-eqn><math-text><italic>x</italic><sub>−</sub></math-text></inline-eqn> and <inline-eqn><math-text><italic>x</italic><sub>+</sub></math-text></inline-eqn> for <inline-eqn><math-text>ρ=0</math-text></inline-eqn> belong to the event horizons, which then represent initial points for the ordinary differential equation (<eqnref linkend="nj355947eqn23">23</eqnref>). Allowed solutions should start from the aforementioned points and remain inside the ergosurface. According to numerical evaluations, carried out for <inline-eqn><math-text>η=1−10<sup>−4</sup></math-text></inline-eqn>, in the range of interest for our assumption that the RIP represents a small perturbation with respect to the refractive index value in the absence of the RIP, the RIP exhibits a somewhat odd behaviour in the sense that the solution does not extend beyond <inline-eqn><math-text><italic>x</italic><sub>−</sub></math-text></inline-eqn> and <inline-eqn><math-text><italic>x</italic><sub>+</sub></math-text></inline-eqn>.</p><p>Physically, a light ray passing through a point lying outside the ergoregion can travel in either direction (toward the RIP, or away from it). Such a ray is not ‘trapped’, and therefore such a point cannot belong to an event horizon.</p><p>Figure <figref linkend="nj355947fig7">7</figref> shows the solutions to equation (<eqnref linkend="nj355947eqn23">23</eqnref>) for the case <inline-eqn><math-text>η=0.1</math-text></inline-eqn>. The dashed line shows the ergosurface predicted by equation (<eqnref linkend="nj355947eqn9">9</eqnref>) and the solid lines are solutions to equation (<eqnref linkend="nj355947eqn23">23</eqnref>) for different initial conditions. These curves represent light-like surfaces passing through points inside the ergosphere (in the RIP frame). The different null hypersurfaces that pass near <inline-eqn><math-text><italic>x</italic><sub>+</sub></math-text></inline-eqn> (or <inline-eqn><math-text><italic>x</italic><sub>−</sub></math-text></inline-eqn>) cross the ergosphere at a point that is progressively closer to <inline-eqn><math-text><italic>x</italic><sub>+</sub></math-text></inline-eqn> (or <inline-eqn><math-text><italic>x</italic><sub>−</sub></math-text></inline-eqn>). This is numerical evidence of the fact that a hypothetical event horizon cannot extend beyond the points <inline-eqn><math-text><italic>x</italic><sub>+</sub></math-text></inline-eqn> and <inline-eqn><math-text><italic>x</italic><sub>−</sub></math-text></inline-eqn>. As can be seen, the solutions that pass through <inline-eqn><math-text><italic>x</italic><sub>+</sub></math-text></inline-eqn> and <inline-eqn><math-text><italic>x</italic><sub>−</sub></math-text></inline-eqn> are tangent to the ergosurface and always remain <italic>outside</italic> the latter so that the horizon effectively reduces to a point. We note that strictly speaking even though these points behave like points of a trapping horizon, as shown for example in figure <figref linkend="nj355947fig3">3</figref>, it is of course not possible to define them as event horizons in the sense of a light-like hypersurface. One may wonder whether a ‘trapping point’ is expected to still emit Hawking-like radiation and, if so, with which features. In particular, can we rely on some arguments in order to support the idea that there is a pair-creation process that, moreover, still shares common features with the Hawking effect? We show in a forthcoming work that the answer is affirmative [<cite linkend="nj355947bib18">18</cite>]. Points where the ergosurface becomes light-like, beyond their properties to be trapping points, which behave as standard horizons, share the further property of involving production of particles with a thermal spectrum, again like the case of standard (non-degenerate) horizons.</p><figure id="nj355947fig7" parts="single" width="column" position="float" pageposition="top" printstyle="normal" orientation="port"><graphic position="indented"><graphic-file version="print" format="EPS" width="18.0pc" printcolour="no" filename="images/nj355947fig7.eps"></graphic-file><graphic-file version="ej" format="JPEG" printcolour="no" filename="images/nj355947fig7.jpg"></graphic-file></graphic><caption type="figure" id="nj355947fc7" label="Figure 7"><p indent="no">Solid lines: real part of the solutions to equation (<eqnref linkend="nj355947eqn23">23</eqnref>). Dashed line: the ergosphere predicted by equation (<eqnref linkend="nj355947eqn9">9</eqnref>). The parameters are <inline-eqn><math-text><italic>n</italic><sub>0</sub>=1.45</math-text></inline-eqn>, <inline-eqn><math-text>η=0.1</math-text></inline-eqn> and <inline-eqn><math-text><italic>v</italic>=<italic>c</italic>/1.52</math-text></inline-eqn>.</p></caption></figure></sec-level1><sec-level1 id="nj355947s6" label="6"><heading>Conclusions</heading><p indent="no">We have derived the effective metric associated with an RIP in a nonlinear Kerr medium. Such a metric simulates the effect on light rays of the presence of a dielectric perturbation in a sample of an isotropic <inline-eqn><math-text>χ<sup>(3)</sup></math-text></inline-eqn> dielectric material. Then, we have investigated the possibility for an ergoregion to arise, and the behaviour of geodesics by means of numerical simulations, also taking into account the presence of dispersion. Finally, we have provided an equation for the event horizon that highlights a rather peculiar geometry, as it appears that the solutions may reduce to a single point. Numerical evaluation in the range of interest confirms the former conclusion. Still, the two special points <inline-eqn><math-text><italic>x</italic><sub>−</sub></math-text></inline-eqn> and <inline-eqn><math-text><italic>x</italic><sub>+</sub></math-text></inline-eqn> of the ergosurface exhibit trapping properties like those of a standard trapping horizon from a classical point of view.</p></sec-level1><appendix id="nj355947a1"><sec-level1 id="nj355947s7" type="num" appendix="yes" style="ALPHA" label="Appendix"><heading>A remark on the extraordinary metric</heading><p indent="no">We show that, under reasonable approximations, the extraordinary metric can be recast in the Gordon form. We work in a reference frame at rest w.r.t. the medium (the laboratory frame). Consider a background field <inline-eqn><math-text><italic>E</italic></math-text></inline-eqn> that, in the laboratory frame, has polarization vector orthogonal to its direction of propagation <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn>. We can choose Cartesian coordinates <inline-eqn><math-text>(<italic>t</italic>, <italic>x</italic>, <italic>y</italic>, <italic>z</italic>)</math-text></inline-eqn> so that <inline-eqn><math-text><italic>l</italic><sup>μ</sup>=(0, <italic>l</italic><sub><italic>x</italic></sub>=0, <italic>l</italic><sub><italic>y</italic></sub>=0, <italic>l</italic><sub><italic>z</italic></sub>=1)</math-text></inline-eqn>. The line element for the extraordinary rays takes the form
<display-eqn id="nj355947eqnA.1" textype="equation" notation="LaTeX" eqnnum="A.1"></display-eqn>Note that the factor in front of <inline-eqn><math-text><upright>d</upright><italic>z</italic><sup>2</sup></math-text></inline-eqn> cannot vanish. Moreover, <inline-eqn><math-text>ξ≪1</math-text></inline-eqn> (see also [<cite linkend="nj355947bib17">17</cite>]). Then, we can safely keep only the leading order term in the <inline-eqn><math-text><italic>z</italic></math-text></inline-eqn> part of the metric so that, in so doing, the extraordinary metric reduces to
<display-eqn id="nj355947eqnA.2" textype="equation" notation="LaTeX" eqnnum="A.2"></display-eqn>Under this approximation, the effective metric for extraordinary rays has again the same form as the Gordon one, with an effective refractive index
<display-eqn id="nj355947eqnA.3" textype="equation" notation="LaTeX" eqnnum="A.3"></display-eqn>to be compared to the refractive index for the ordinary rays, which is
<display-eqn id="nj355947eqnA.4" textype="equation" notation="LaTeX" eqnnum="A.4"></display-eqn>Here, <inline-eqn><math-text><italic>n</italic><sub>0</sub></math-text></inline-eqn> is the refractive index of the medium when we disregard nonlinear effects, i.e.
<display-eqn id="nj355947eqnA.5" textype="equation" notation="LaTeX" eqnnum="A.5"></display-eqn>Then, the extraordinary rays ‘feel’ a disturbance of the refractive index three times stronger than the ordinary ones. The same conclusion can be also achieved in the context of Hawking radiation production by an RIP [<cite linkend="nj355947bib18">18</cite>]. Indeed, it can be shown that the differential equation that defines the extraordinary field modes is well approximated by the same equation as for the ordinary ones, apart from an effective RIP three times larger.</p></sec-level1></appendix></body><back><references><heading>References</heading><reference-list type="numeric"><journal-ref id="nj355947bib1" num="1"><authors><au><second-name>Gordon</second-name><first-names>W</first-names></au></authors><year>1923</year><art-title>Zur Lichtfortpflanzung nach der relativitätstheorie</art-title><jnl-title>Ann. Phys.</jnl-title><volume>377</volume><pages>421–56</pages><crossref><cr_doi>http://dx.doi.org/10.1002/andp.19233772202</cr_doi><cr_issn type="print">00033804</cr_issn><cr_issn type="electronic">15213889</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib2" num="2"><authors><au><second-name>Leonhardt</second-name><first-names>U</first-names></au><au><second-name>Piwnicki</second-name><first-names>P</first-names></au></authors><year>1999</year><art-title>Optics of nonuniformly moving media</art-title><jnl-title>Phys. Rev.</jnl-title><part>A</part><volume>60</volume><pages>4301–12</pages><crossref><cr_doi>http://dx.doi.org/10.1103/PhysRevA.60.4301</cr_doi><cr_issn type="print">10502947</cr_issn><cr_issn type="electronic">10941622</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib3" num="3"><authors><au><second-name>Leonhardt</second-name><first-names>U</first-names></au><au><second-name>Piwnicki</second-name><first-names>P</first-names></au></authors><year>2000</year><art-title>Relativistic effects of light in moving media with extremely low group velocity</art-title><jnl-title>Phys. Rev. Lett.</jnl-title><volume>84</volume><pages>822–5</pages><crossref><cr_doi>http://dx.doi.org/10.1103/PhysRevLett.84.822</cr_doi><cr_issn type="print">00319007</cr_issn><cr_issn type="electronic">10797114</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib4" num="4"><authors><au><second-name>Schützhold</second-name><first-names>R</first-names></au><au><second-name>Plunien</second-name><first-names>G</first-names></au><au><second-name>Soff</second-name><first-names>G</first-names></au></authors><year>2002</year><art-title>Dielectric black hole analogs</art-title><jnl-title>Phys. Rev. Lett.</jnl-title><volume>88</volume><pages>061101</pages><crossref><cr_doi>http://dx.doi.org/10.1103/PhysRevLett.88.061101</cr_doi><cr_issn type="print">00319007</cr_issn><cr_issn type="electronic">10797114</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib5" num="5"><authors><au><second-name>Brevik</second-name><first-names>I H</first-names></au><au><second-name>Halnes</second-name><first-names>G</first-names></au></authors><year>2001</year><art-title>Light rays at optical black holes in moving media</art-title><jnl-title>Phys. Rev.</jnl-title><part>D</part><volume>65</volume><pages>024005</pages><crossref><cr_doi>http://dx.doi.org/10.1103/PhysRevD.65.024005</cr_doi><cr_issn type="print">15507998</cr_issn><cr_issn type="electronic">10894918</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib6" num="6"><authors><au><second-name>De Lorenci</second-name><first-names>V A</first-names></au><au><second-name>Klippert</second-name><first-names>R</first-names></au></authors><year>2002</year><art-title>Analogue gravity from electrodynamics in nonlinear media</art-title><jnl-title>Phys. Rev.</jnl-title><part>D</part><volume>65</volume><pages>064027</pages><crossref><cr_doi>http://dx.doi.org/10.1103/PhysRevD.65.064027</cr_doi><cr_issn type="print">15507998</cr_issn><cr_issn type="electronic">10894918</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib7" num="7"><authors><au><second-name>De Lorenci</second-name><first-names>V A</first-names></au><au><second-name>Souza</second-name><first-names>M A</first-names></au></authors><year>2001</year><art-title>Electromagnetic wave propagation inside a material medium: an effective geometry interpretation</art-title><jnl-title>Phys. Lett.</jnl-title><part>B</part><volume>512</volume><pages>417–22</pages><crossref><cr_doi>http://dx.doi.org/10.1016/S0370-2693(01)00588-3</cr_doi><cr_issn type="print">03702693</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib8" num="8"><authors><au><second-name>Novello</second-name><first-names>M</first-names></au><au><second-name>Salim</second-name><first-names>J M</first-names></au></authors><year>2001</year><art-title>Effective electromagnetic geometry</art-title><jnl-title>Phys. Rev.</jnl-title><part>D</part><volume>63</volume><pages>083511</pages><crossref><cr_doi>http://dx.doi.org/10.1103/PhysRevD.63.083511</cr_doi><cr_issn type="print">15507998</cr_issn><cr_issn type="electronic">10894918</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib9" num="9"><authors><au><second-name>Marklund</second-name><first-names>M</first-names></au><au><second-name>Anderson</second-name><first-names>D</first-names></au><au><second-name>Cattani</second-name><first-names>F</first-names></au><au><second-name>Lisak</second-name><first-names>M</first-names></au><au><second-name>Lundgren</second-name><first-names>L</first-names></au></authors><year>2002</year><art-title>Fermat's principle and the variational analysis of an optical model for light propagation exhibiting a critical radius</art-title><jnl-title>Am. J. Phys.</jnl-title><volume>70</volume><pages>680–3</pages><crossref><cr_doi>http://dx.doi.org/10.1119/1.1477429</cr_doi><cr_issn type="print">00029505</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib10" num="10"><authors><au><second-name>De Lorenci</second-name><first-names>V A</first-names></au><au><second-name>Klippert</second-name><first-names>R</first-names></au><au><second-name>Obukhov Yu</second-name><first-names>N</first-names></au></authors><year>2003</year><art-title>On optical black holes in moving dielectrics</art-title><jnl-title>Phys. Rev.</jnl-title><part>D</part><volume>68</volume><pages>061502</pages><crossref><cr_doi>http://dx.doi.org/10.1103/PhysRevD.68.061502</cr_doi><cr_issn type="print">15507998</cr_issn><cr_issn type="electronic">10894918</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib11" num="11"><authors><au><second-name>Novello</second-name><first-names>M</first-names></au><au><second-name>Perez</second-name><first-names>B</first-names></au><au><second-name>Santiago</second-name><first-names>E</first-names></au><au><second-name>Salim</second-name><first-names>J</first-names></au><au><second-name>De Lorenci</second-name><first-names>V</first-names></au><au><second-name>Klippert</second-name><first-names>R</first-names></au></authors><year>2003</year><art-title>Analog black holes in flowing dielectrics</art-title><jnl-title>Class. Quantum Gravity</jnl-title><volume>20</volume><pages>859–72</pages><crossref><cr_doi>http://dx.doi.org/10.1088/0264-9381/20/5/306</cr_doi><cr_issn type="print">02649381</cr_issn><cr_issn type="electronic">13616382</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib12" num="12"><authors><au><second-name>Novello</second-name><first-names>M</first-names></au><au><second-name>Perez</second-name><first-names>B</first-names></au><au><second-name>Santiago</second-name><first-names>E</first-names></au></authors><year>2003</year><art-title>Effective geometry</art-title><jnl-title>AIP Conf. Proc.</jnl-title><volume>668</volume><pages>288–300</pages><crossref><cr_doi>http://dx.doi.org/10.1063/1.1587103</cr_doi></crossref></journal-ref><journal-ref id="nj355947bib13" num="13"><authors><au><second-name>Unruh</second-name><first-names>W G</first-names></au><au><second-name>Schützhold</second-name><first-names>R</first-names></au></authors><year>2003</year><art-title>On slow light as a black hole analogue</art-title><jnl-title>Phys. Rev.</jnl-title><part>D</part><volume>68</volume><pages>024008</pages><crossref><cr_doi>http://dx.doi.org/10.1103/PhysRevD.68.024008</cr_doi><cr_issn type="print">15507998</cr_issn><cr_issn type="electronic">10894918</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib14" num="14"><authors><au><second-name>Schutzhold</second-name><first-names>R</first-names></au><au><second-name>Unruh</second-name><first-names>W G</first-names></au></authors><year>2005</year><art-title>Hawking radiation in an electro-magnetic waveguide?</art-title><jnl-title>Phys. Rev. Lett.</jnl-title><volume>95</volume><pages>031301</pages><crossref><cr_doi>http://dx.doi.org/10.1103/PhysRevLett.95.031301</cr_doi><cr_issn type="print">00319007</cr_issn><cr_issn type="electronic">10797114</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib15" num="15"><authors><au><second-name>Philbin</second-name><first-names>T G</first-names></au><au><second-name>Kuklewicz</second-name><first-names>C</first-names></au><au><second-name>Robertson</second-name><first-names>S</first-names></au><au><second-name>Hill</second-name><first-names>S</first-names></au><au><second-name>Konig</second-name><first-names>F</first-names></au><au><second-name>Leonhardt</second-name><first-names>U</first-names></au></authors><year>2008</year><art-title>Fiber-optical analog of the event horizon</art-title><jnl-title>Science</jnl-title><volume>319</volume><pages>1367</pages><crossref><cr_doi>http://dx.doi.org/10.1126/science.1153625</cr_doi><cr_issn type="print">00368075</cr_issn><cr_issn type="electronic">10959203</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib16" num="16"><authors><au><second-name>Faccio</second-name><first-names>D</first-names></au><au><second-name>Cacciatori</second-name><first-names>S</first-names></au><au><second-name>Gorini</second-name><first-names>V</first-names></au><au><second-name>Sala</second-name><first-names>V G</first-names></au><au><second-name>Averchi</second-name><first-names>A</first-names></au><au><second-name>Lotti</second-name><first-names>A</first-names></au><au><second-name>Kolesik</second-name><first-names>M</first-names></au><au><second-name>Moloney</second-name><first-names>J V</first-names></au></authors><year>2010</year><art-title>Analogue gravity and ultrashort laser pulse filamentation</art-title><jnl-title>Europhys. Lett.</jnl-title><volume>89</volume><pages>34004</pages><crossref><cr_doi>http://dx.doi.org/10.1209/0295-5075/89/34004</cr_doi><cr_issn type="print">02955075</cr_issn><cr_issn type="electronic">12864854</cr_issn></crossref></journal-ref><journal-ref id="nj355947bib17" num="17"><authors><au><second-name>Belgiorno</second-name><first-names>F</first-names></au><au><second-name>Cacciatori</second-name><first-names>S L</first-names></au><au><second-name>Gorini</second-name><first-names>V</first-names></au><au><second-name>Ortenzi</second-name><first-names>G</first-names></au><au><second-name>Rizzi</second-name><first-names>L</first-names></au><au><second-name>Faccio</second-name><first-names>D</first-names></au></authors><year>2010</year><art-title>Dielectric black holes induced by a refractive index perturbation and the Hawking effect</art-title><jnl-title>Phys. Rev.</jnl-title><part>D</part><misc-text>submitted (arXiv:<arxiv url="1003.4150">1003.4150</arxiv>)</misc-text></journal-ref><misc-ref id="nj355947bib18" num="18"><authors><au><second-name>Belgiorno</second-name><first-names>F</first-names></au><au><second-name>Cacciatori</second-name><first-names>S L</first-names></au><au><second-name>Gorini</second-name><first-names>V</first-names></au><au><second-name>Ortenzi</second-name><first-names>G</first-names></au><au><second-name>Rizzi</second-name><first-names>L</first-names></au><au><second-name>Faccio</second-name><first-names>D.</first-names></au></authors><year>2010</year><misc-text>in preparation</misc-text></misc-ref><journal-ref id="nj355947bib19" num="19"><authors><au><second-name>Robertson</second-name><first-names>S</first-names></au><au><second-name>Leonhardt</second-name><first-names>U</first-names></au></authors><year>2010</year><art-title>Frequency shifting at fiber-optical event horizons: the effect of Raman deceleration</art-title><jnl-title>Phys. Rev.</jnl-title><part>A</part><volume>81</volume><pages>063835</pages><crossref><cr_doi>http://dx.doi.org/10.1103/PhysRevA.81.063835</cr_doi><cr_issn type="print">10502947</cr_issn><cr_issn type="electronic">10941622</cr_issn></crossref></journal-ref><book-ref id="nj355947bib20" num="20"><authors><au><second-name>Wald</second-name><first-names>R M</first-names></au></authors><year>1984</year><book-title>General Relativity</book-title><publication><place>Chicago, IL</place><publisher>University of Chicago Press</publisher></publication></book-ref><journal-ref id="nj355947bib21" num="21"><authors><au><second-name>Visser</second-name><first-names>M</first-names></au></authors><year>1998</year><jnl-title>Class. Quantum Gravity</jnl-title><volume>15</volume><pages>1767–91</pages><crossref><cr_doi>http://dx.doi.org/10.1088/0264-9381/15/6/024</cr_doi><cr_issn type="print">02649381</cr_issn><cr_issn type="electronic">13616382</cr_issn></crossref></journal-ref></reference-list></references><footnotes><footnote id="nj355947fn7" source="text" type="number" marker="7"><p>For reasons connected with the actual experimental verification of this theoretical framework, we may think of this strong background field as a spatially localized (albeit rather extended) pulse, which is travelling through the medium.</p></footnote><footnote id="nj355947fn8" source="text" type="number" marker="8"><p>A simple, non-technical definition of static and stationary spacetimes can be given as follows. A spacetime is stationary if it is possible to find a reference frame in which the metric coefficients do not depend on time. This means that time translation is a symmetry of the spacetime. However, such a spacetime might not be invariant under time reversal. In the case of time-reversal invariance, the spacetime is not only stationary but also static. An example of static spacetime is the Schwarzschild solution, describing the static empty spacetime outside a non-rotating black hole. If the black hole has a non-vanishing angular momentum the spacetime is still invariant under time translations; however, time reversal reverses the sign of the angular momentum, turning a clockwise spinning black hole to an anticlockwise one. Then, the spacetime around a spinning black hole is stationary, but not static.</p></footnote><footnote id="nj355947fn9" source="text" type="number" marker="9"><p>This is true also when taking into account the material dispersion, which is neglected in this subsection.</p></footnote></footnotes></back></article></istex:document></istex:metadataXml><mods version="3.6"><titleInfo lang="eng"><title>Spacetime geometries and light trapping in travelling refractive index perturbations</title></titleInfo><titleInfo type="abbreviated"><title>Spacetime geometries and light trapping in travelling refractive index perturbations</title></titleInfo><titleInfo type="alternative" lang="eng"><title>Spacetime geometries and light trapping in travelling refractive index perturbations</title></titleInfo><name type="personal"><namePart type="given">S L</namePart><namePart type="family">Cacciatori</namePart><affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</affiliation><affiliation>INFN sezione di Milano, via Celoria 16, IT-20133 Milano, Italy</affiliation><affiliation>Author to whom any correspondence should be addressed.</affiliation><affiliation>E-mail: sergio.cacciatori@uninsubria.it</affiliation><role><roleTerm type="text">author</roleTerm></role></name><name type="personal"><namePart type="given">F</namePart><namePart type="family">Belgiorno</namePart><affiliation>Dipartimento di Fisica, Universit di Milano, Via Celoria 16, IT-20133 Milano, Italy</affiliation><role><roleTerm type="text">author</roleTerm></role></name><name type="personal"><namePart type="given">V</namePart><namePart type="family">Gorini</namePart><affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</affiliation><affiliation>INFN sezione di Milano, via Celoria 16, IT-20133 Milano, Italy</affiliation><role><roleTerm type="text">author</roleTerm></role></name><name type="personal"><namePart type="given">G</namePart><namePart type="family">Ortenzi</namePart><affiliation>Dipartimento di Matematica e Applicazioni, Universitdegli Studi di Milano, Bicocca, IT-20125 Milano, Italy</affiliation><role><roleTerm type="text">author</roleTerm></role></name><name type="personal"><namePart type="given">L</namePart><namePart type="family">Rizzi</namePart><affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</affiliation><role><roleTerm type="text">author</roleTerm></role></name><name type="personal"><namePart type="given">V G</namePart><namePart type="family">Sala</namePart><affiliation>Department of Physics and Mathematics, Universit dell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</affiliation><role><roleTerm type="text">author</roleTerm></role></name><name type="personal"><namePart type="given">D</namePart><namePart type="family">Faccio</namePart><affiliation>CNISM and Department of Physics and Mathematics, Universitdell'Insubria, Via Valleggio 11, IT-22100 Como, Italy</affiliation><role><roleTerm type="text">author</roleTerm></role></name><typeOfResource>text</typeOfResource><genre type="research-article" displayLabel="paper">paper</genre><originInfo><publisher>IOP Publishing</publisher><dateIssued encoding="w3cdtf">2010</dateIssued><copyrightDate encoding="w3cdtf">2010</copyrightDate></originInfo><language><languageTerm type="code" authority="iso639-2b">eng</languageTerm><languageTerm type="code" authority="rfc3066">en</languageTerm></language><physicalDescription><internetMediaType>text/html</internetMediaType></physicalDescription><abstract>In the framework of transformation optics, we show that the propagation of a locally superluminal refractive index perturbation (RIP) in a Kerr medium can be described, in the eikonal approximation, by means of a stationary metric, which we prove to be of Gordon type. Under suitable hypotheses on the RIP, we obtain a stationary but not static metric, which is characterized by an ergosphere and by a peculiar behaviour of the geodesics, which are studied numerically, also accounting for material dispersion. Finally, the equation to be satisfied by an event horizon is also displayed and briefly discussed.</abstract><relatedItem type="host"><titleInfo><title>New Journal of Physics</title></titleInfo><titleInfo type="abbreviated"><title>New J. Phys.</title></titleInfo><genre type="journal">journal</genre><identifier type="ISSN">1367-2630</identifier><identifier type="eISSN">1367-2630</identifier><identifier type="PublisherID">nj</identifier><identifier type="CODEN">NJOPFM</identifier><identifier type="URL">stacks.iop.org/NJP</identifier><part><date>2010</date><detail type="volume"><caption>vol.</caption><number>12</number></detail><detail type="issue"><caption>no.</caption><number>9</number></detail><extent unit="pages"><start>1</start><end>14</end><total>14</total></extent></part></relatedItem><identifier type="istex">DB3DCA224C4A8D0CF41243DB5493DA9F659E1357</identifier><identifier type="DOI">10.1088/1367-2630/12/9/095021</identifier><identifier type="articleID">355947</identifier><identifier type="articleNumber">095021</identifier><accessCondition type="use and reproduction" contentType="copyright">IOP Publishing and Deutsche Physikalische Gesellschaft</accessCondition><recordInfo><recordContentSource>IOP</recordContentSource><recordOrigin>IOP Publishing and Deutsche Physikalische Gesellschaft</recordOrigin></recordInfo></mods></metadata><enrichments><json:item><type>multicat</type><uri>https://api.istex.fr/document/DB3DCA224C4A8D0CF41243DB5493DA9F659E1357/enrichments/multicat</uri></json:item></enrichments><serie></serie></istex></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Belgique/explor/OpenAccessBelV2/Data/Istex/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 001362 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Istex/Corpus/biblio.hfd -nk 001362 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Belgique
|area= OpenAccessBelV2
|flux= Istex
|étape= Corpus
|type= RBID
|clé= ISTEX:DB3DCA224C4A8D0CF41243DB5493DA9F659E1357
|texte= Spacetime geometries and light trapping in travelling refractive index perturbations
}}
| This area was generated with Dilib version V0.6.25. |  |