Source author record

Emmanuel Lorin

Emmanuel Lorin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Computation of the Time-Dependent Dirac Equation with Physics-Informed Neural Networks

We propose to compute the time-dependent Dirac equation using physics-informed neural networks (PINNs), a new powerful tool in scientific machine learning avoiding the use of approximate derivatives of differential operators. PINNs search solutions in the form of parameterized (deep) neural networks, whose derivatives (in time and space) are performed by automatic differentiation. The computational cost comes from the need to solve high-dimensional optimization problems using stochastic gradient methods and train the network with a large number of points. Specifically, we derive PINNs-based algorithms and present some key fundamental properties of these algorithms when applied to the Dirac equations in different physical frameworks.

preprint2019arXiv

Pseudospectral computational methods for the time-dependent Dirac equation in static curved spaces

Pseudospectral numerical schemes for solving the Dirac equation in general static curved space are derived using a pseudodifferential representation of the Dirac equation along with a simple Fourier-basis technique. Owing to the presence of non-constant coefficients in the curved space Dirac equation, convolution products usually appear when the Fourier transform is performed. The strategy based on pseudodifferential operators allows for efficient computations of these convolution products by using an ordinary fast Fourier transform algorithm. The resulting numerical methods are efficient and have spectral convergence. Simultaneously, wave absorption at the boundary can be easily derived using absorbing layers to cope with some potential negative effects of periodic conditions inherent to spectral methods. The numerical schemes are first tested on simple systems to verify the convergence and are then applied to the dynamics of charge carriers in strained graphene.

preprint2013arXiv

A Split-Step Numerical Method for the Time-Dependent Dirac Equation in 3-D Axisymmetric Geometry

A numerical method is developed to solve the time-dependent Dirac equation in cylindrical coordinates for 3-D axisymmetric systems. The time evolution is treated by a splitting scheme in coordinate space using alternate direction iteration, while the wave function is discretized spatially on a uniform grid. The longitudinal coordinate evolution is performed exactly by the method of characteristics while the radial coordinates evolution uses Poisson's integral solution, which allows to implement the radial symmetry of the wave function. The latter is evaluated on a time staggered mesh by using Hermite polynomial interpolation and by performing the integration analytically. The cylindrical coordinate singularity problem at $r=0$ is circumvented by this method as the integral is well-defined at the origin. The resulting scheme is reminiscent of non-standard finite differences. In the last step of the splitting, the remaining equation has a solution in terms of a time-ordered exponential, which is approximated to a higher order than the time evolution scheme. We study the time evolution of Gaussian wave packets, and we evaluate the eigenstates of hydrogen-like systems by using a spectral method. We compare the numerical results to analytical solutions to validate the method. In addition, we present three-dimensional simulations of relativistic laser-matter interactions, using the Dirac equation.

preprint2013arXiv

Atoms and Molecules in Intense Laser Fields: Gauge Invariance of Theory and Models

Gauge invariance was discovered in the development of classical electromagnetism and was required when the latter was formulated in terms of the scalar and vector potentials. It is now considered to be a fundamental principle of nature, stating that different forms of these potentials yield the same physical description: they describe the same electromagnetic field as long as they are related to each other by gauge transformations. Gauge invariance can also be included into the quantum description of matter interacting with an electromagnetic field by assuming that the wave function transforms under a given local unitary transformation. The result of this procedure is a quantum theory describing the coupling of electrons, nuclei and photons. Therefore, it is a very important concept: it is used in almost every fields of physics and it has been generalized to describe electroweak and strong interactions in the standard model of particles. A review of quantum mechanical gauge invariance and general unitary transformations is presented for atoms and molecules in interaction with intense short laser pulses, spanning the perturbative to highly nonlinear nonperturbative interaction regimes. Various unitary transformations for single spinless particle Time Dependent Schrödinger Equations, TDSE, are shown to correspond to different time-dependent Hamiltonians and wave functions. Accuracy of approximation methods involved in solutions of TDSE's such as perturbation theory and popular numerical methods depend on gauge or representation choices which can be more convenient due to faster convergence criteria. We focus on three main representations: length and velocity gauges, in addition to the acceleration form which is not a gauge, to describe perturbative and nonperturbative radiative interactions. Numerical schemes for solving TDSE's in different representations are also discussed.

preprint2013arXiv

Enhanced Schwinger's pair production in many-center systems

Electron-positron pair production is considered for many-center systems with multiple bare nuclei immersed in a constant electric field. It is shown that there are two distinct regimes where the pair production rate is enhanced. At small interatomic distance, the effective charge of the nuclei approaches the critical charge where the ground state dives into the negative continuum. This facilitates the transition from the negative to the positive energy states, which in turns, increases the pair production rate. At larger atomic distance, the enhancement is due to the crossing of resonances and the pair production proceeds by the Resonantly Enhanced Pair Production (REPP) mechanism. These processes are studied within a simple one-dimensional model. A numerical method is developed to evaluate the transmission coefficient in relativistic quantum mechanics, which is required in the calculation of the pair production rate. The latter is evaluated for systems with many (up to five) nuclei. It is shown that the production rate for many-center systems can reach a few orders of magnitude above Schwinger's tunnelling result in a static field.

preprint2013arXiv

On the Saturation of the Nonlinear Refractive Index in Atomic Gases

Motivated by the ongoing controversy on the origin of the nonlinear index saturation and subsequent intensity clamping in femtosecond filaments, we study the atomic nonlinear polarization induced by a high-intensity and ultrashort laser pulse in hydrogen by numerically solving the time dependent Schrödinger equation. Special emphasis is given to the efficient modeling of the nonlinear polarization at central laser frequency corresponding to 800 nm wavelength. Here, the recently proposed model of the Higher-Order Kerr Effect (HOKE) and two versions of the Standard model for femtosecond filamentation, including either a multi-photon or tunnel ionization rate, are compared. We find that around the clamping intensity the instantaneous HOKE model does not reproduce the temporal structure of the nonlinear response obtained from the quantum mechanical results. In contrast, the non-instantaneous charge contributions included in the Standard models ensure a reasonable quantitative agreement. Therefore, the physical origin for the observed saturation of the overall electron response is confirmed to mainly result from contributions of free or nearly free electrons.

preprint2012arXiv

Landau-Zener-Stückelberg interferometry in pair production from counterpropagating lasers

The rate of electron-positron pair production in linearly polarized counter-propagating lasers is evaluated from a recently discovered solution of the time-dependent Dirac equation. The latter is solved in momentum space where it is formally equivalent to the Schrödinger equation describing a strongly driven two-level system. The solution is found from a simple transformation of the Dirac equation and is given in compact form in terms of the doubly-confluent Heun's function. By using the analogy with the two-level system, it is shown that for high-intensity lasers, pair production occurs through periodic non-adiabatic transitions when the adiabatic energy gap is minimal. These transitions give rise to an intricate interference pattern in the pair spectrum, reminiscent of the Landau-Zener-Stückelberg phenomenon in molecular physics: the accumulated phase result in constructive or destructive interference. The adiabatic-impulse model is used to study this phenomenon and shows an excellent agreement with the exact result.

preprint2012arXiv

Numerical Solution of the Time-Dependent Dirac Equation in Coordinate Space without Fermion-Doubling

The validation and parallel implementation of a numerical method for the solution of the time-dependent Dirac equation is presented. This numerical method is based on a split operator scheme where the space-time dependence is computed in coordinate space using the method of characteristics. Thus, most of the steps in the splitting are calculated exactly, making for a very efficient and unconditionally stable method. We show that it is free from spurious solutions related to the fermion-doubling problem and that it can be parallelized very efficiently. We consider a few simple physical systems such as the time evolution of Gaussian wave packets and the Klein paradox. The numerical results obtained are compared to analytical formulas for the validation of the method.

preprint2012arXiv

Relativistic Stark resonances in a simple exactly soluble model for a diatomic molecule

A simple 1-D relativistic model for a diatomic molecule with a double point interaction potential is solved exactly in a constant electric field. The Weyl-Titchmarsh-Kodaira method is used to evaluate the spectral density function, allowing the correct normalization of continuum states. The boundary conditions at the potential wells are evaluated using Colombeau's generalized function theory along with charge conjugation invariance and general properties of self-adjoint extensions for point-like interactions. The resulting spectral density function exhibits resonances for quasibound states which move in the complex energy plane as the model parameters are varied. It is observed that for a monotonically increasing interatomic distance, the ground state resonance can either go deeper into the negative continuum or can give rise to a sequence of avoided crossings, depending on the strength of the potential wells. For sufficiently low electric field strength or small interatomic distance, the behavior of resonances is qualitatively similar to non-relativistic results.

preprint2012arXiv

Resonantly enhanced pair production in a simple diatomic model

A new mechanism for the production of electron-positron pairs from the interaction of a laser field and a fully stripped diatomic molecule in the tunneling regime is presented. When the laser field is turned off, the Dirac operator has resonances in both the positive and the negative energy continua while bound states are in the mass gap. When this system is immersed in a strong laser field, the resonances move in the complex energy plane: the negative energy resonances are pushed to higher energies while the bound states are Stark shifted. It is argued here that there is a pair production enhancement at the crossing of resonances by looking at a simple 1-D model: the nuclei are modeled simply by Dirac delta potential wells while the laser field is assumed to be static and of finite spatial extent. The average rate for the number of electron-positron pairs produced is evaluated and the results are compared to the single nucleus and to the free cases. It is shown that positrons are produced by the Resonantly Enhanced Pair Production (REPP) mechanism, which is analogous to the resonantly enhanced ionization of molecular physics. This phenomenon could be used to increase the number of pairs produced at low field strength, allowing the study of the Dirac vacuum.

preprint2011arXiv

Efficient determination of critical parameters of nonlinear Schrödinger equation with point-like potential using generalized polynomial chaos methods

We consider the nonlinear Schrödinger equation with a point-like source term. The soliton interaction with such a singular potential yields a critical solution behavior. That is, for the given value of the potential strength and the soliton amplitude, there exists a critical velocity of the initial soliton solution, around which the solution is either trapped by or transmitted through the potential. In this paper, we propose an efficient method for finding such a critical velocity by using the generalized polynomial chaos method. For the proposed method, we assume that the soliton velocity is a random variable and expand the solution in the random space using the orthogonal polynomials. The proposed method finds the critical velocity accurately with spectral convergence. Thus the computational complexity is much reduced. Numerical results for the smaller and higher values of the potential strength confirm the spectral convergence of the proposed method.