-
Transitions in Xenes between excitonic, topological and trivial insulator phases: influence of screening, band dispersion and external electric field
Authors:
Olivia Pulci,
Paola Gori,
Davide Grassano,
Marco D'Alessandro,
Friedhelm Bechstedt
Abstract:
Using a variational approach, the binding energies $E_b$ of the lowest bound excitons in Xenes under varying electric field are investigated. The internal exciton motion is described both by Dirac electron dispersion and in effective-mass approximation, while the screened electron-hole attraction is modeled by a Rytova-Keldysh potential with a 2D electronic polarizability $α_{2{\rm D}}$. The most…
▽ More
Using a variational approach, the binding energies $E_b$ of the lowest bound excitons in Xenes under varying electric field are investigated. The internal exciton motion is described both by Dirac electron dispersion and in effective-mass approximation, while the screened electron-hole attraction is modeled by a Rytova-Keldysh potential with a 2D electronic polarizability $α_{2{\rm D}}$. The most important parameters as spin-orbit-induced gap $E_g$, Fermi velocity $v_F$ and $α_{2{\rm D}}$ are taken from ab initio density functional theory calculations. In addition, $α_{2{\rm D}}$ is approximated in two different ways. The relation of $E_b$ and $E_g$ is ruled by the screening. The existence of an excitonic insulator phase with $E_b>E_g$ sensitively depends on the chosen $α_{2{\rm D}}$. The values of $E_g$ and $α_{2{\rm D}}$ are strongly modified by a vertical external electric bias $U$, which defines a transition from the topological into a trivial insulator at $U=E_g/2$, with the exception of plumbene. Within the Dirac approximation, but also within the effective mass description of the kinetic energy, the treatment of screening dominates the appearance or non-appearance of an excitonic insulator phase. Gating does not change the results: the prediction done at zero electric field is confirmed when a vertical electric field is applied. Finally, Many-Body perturbation theory approaches based on the Green's function method, applied to stanene, confirm the absence of an excitonic insulator phase, thus validating our results obtained by ab initio modeling of $α_{2{\rm D}}$.
△ Less
Submitted 27 April, 2023; v1 submitted 20 January, 2023;
originally announced January 2023.
-
Exciton-Exciton transitions involving strongly bound excitons: an ab initio approach
Authors:
D. Sangalli,
M. D'Alessandro,
C. Attaccalite
Abstract:
In pump-probe spectroscopy, two laser pulses are employed to garner dynamical information from the sample of interest. The pump initiates the optical process by exciting a portion of the sample from the electronic ground state to an accessible electronic excited state, an exciton. Thereafter, the probe interacts with the already excited sample. The change in the absorbance after pump provides info…
▽ More
In pump-probe spectroscopy, two laser pulses are employed to garner dynamical information from the sample of interest. The pump initiates the optical process by exciting a portion of the sample from the electronic ground state to an accessible electronic excited state, an exciton. Thereafter, the probe interacts with the already excited sample. The change in the absorbance after pump provides information on transitions between the excited states and their dynamics. In this work we study these exciton-exciton transitions by means of an ab initio real time propagation scheme based on dynamical Berry phase formulation. The results are then analyzed taking advantage of a Fermi-golden rule approach formulated in the excitonic basis-set and in terms of the symmetries of the excitonic states. Using bulk LiF and 2D hBN as two prototype materials, we discuss the selection rules for transitions involving strongly bound excitons, for which the hydrogen model cannot be used.
△ Less
Submitted 15 April, 2023; v1 submitted 22 November, 2022;
originally announced November 2022.
-
Real-time ab initio description of the photon-echo mechanisms in extended systems: the case study of bulk GaAs
Authors:
M. D'Alessandro,
D. Sangalli
Abstract:
In this paper we present an ab initio real-time analysis of free polarization decay and photon echo in extended systems. As a prototype material, we study bulk GaAs driven by ultra-short laser pulses of 10 fs (energy spread of 0.4 eV), with frequency tuned in the continuum of the optical spectrum. We compute the electronic polarization P(t), and define a computational procedure to extract the echo…
▽ More
In this paper we present an ab initio real-time analysis of free polarization decay and photon echo in extended systems. As a prototype material, we study bulk GaAs driven by ultra-short laser pulses of 10 fs (energy spread of 0.4 eV), with frequency tuned in the continuum of the optical spectrum. We compute the electronic polarization P(t), and define a computational procedure to extract the echo signal in the dipole approximation. Results are obtained in both the low and high field regime, and compared with a two-levels system (TLS) model, with parameters extracted from the ab initio simulations. ab initio results are in optimal agreement with the TLS in the low-field case, whereas some differences are observed in the high-field regime where the multi-band nature of GaAs becomes relevant. In the high field regime we compute the pulse area, and look for fluences with pulse area close to π. We highlight that such fluences are well below the damage threshold of GaAs. However a unique value of the area cannot be defined, due to the strong dependence of the transition dipoles in the energy window excited by the laser pulse.
△ Less
Submitted 22 April, 2022; v1 submitted 27 December, 2021;
originally announced December 2021.
-
Real-time modelling of Optical orientation in GaAs: generation and decay of the degree of spin polarization
Authors:
Marco D'Alessandro,
Davide Sangalli
Abstract:
We present a real-time abinitio description of optical orientation in bulk GaAs due to the coupling with an ultrashort circular polarized laser source. The injection of spin polarized electrons in the conduction band is correctly reproduced, and a non vanishing spin polarization ($\mathbf{P}$) parallel to the direction of propagation of the laser ($z$) emerges. A detailed analysis of the generatio…
▽ More
We present a real-time abinitio description of optical orientation in bulk GaAs due to the coupling with an ultrashort circular polarized laser source. The injection of spin polarized electrons in the conduction band is correctly reproduced, and a non vanishing spin polarization ($\mathbf{P}$) parallel to the direction of propagation of the laser ($z$) emerges. A detailed analysis of the generation and the evolution of $\mathbf{P}(t)$ is discussed. The single $\mathbf{k}$-point dynamics is a motion of precession around a fixed axis with constant $|\mathbf{P}|$ and fixed frequency. Instead, the $\mathbf{k}$-integrated signal shows only a time dependent $P_z(t)$ and decays few pico seconds after the end of the laser pump due to decoherence. Decoherence emerges since the individual contributions activated by the pump give rise to destructive interference. We interpret the results in terms of the \emph{free induction decay} mechanism proposed some years ago. For the first time we are able to reproduce such effect in a full abinitio fashion, giving a quantitative estimate of the associated decay time. Our result also shows a possible explanation to the time decay of spin magnetization observed in many real-time abinitio simulations.
△ Less
Submitted 12 June, 2020;
originally announced June 2020.
-
Work function, deformation potential, and collapse of Landau levels in strained graphene and silicene
Authors:
D. Grassano,
M. D'Alessandro,
O. Pulci,
S. G. Sharapov,
V. P. Gusynin,
A. A. Varlamov
Abstract:
We perform a systematic {\it ab initio} study of the work function and its uniform strain dependence for graphene and silicene for both tensile and compressive strains. The Poisson ratios associated with armchair and zigzag strains are also computed. Based on these results, we obtain the deformation potential, crucial for straintronics, as a function of the applied strain. Further, we propose a pa…
▽ More
We perform a systematic {\it ab initio} study of the work function and its uniform strain dependence for graphene and silicene for both tensile and compressive strains. The Poisson ratios associated with armchair and zigzag strains are also computed. Based on these results, we obtain the deformation potential, crucial for straintronics, as a function of the applied strain. Further, we propose a particular experimental setup with a special strain configuration that generates only the electric field, while the pseudomagnetic field is absent. Then, applying a real magnetic field, one should be able to realize experimentally the spectacular phenomenon of the collapse of Landau levels in graphene or related two-dimensional materials.
△ Less
Submitted 4 June, 2020; v1 submitted 24 January, 2020;
originally announced January 2020.
-
Locality and Computational Reliability of Linear Response Calculations for Molecular Systems
Authors:
Marco D'Alessandro,
Luigi Genovese
Abstract:
By performing a critical analysis of the fundamental equations of linear-response (LR) formalism in molecules, we explore the interplay between locality of the response density operator and numerical convergence of LR-related quantities. We show that for frequencies below the first ionization potential (IP) of the system, it is possible to express the response density by employing localized states…
▽ More
By performing a critical analysis of the fundamental equations of linear-response (LR) formalism in molecules, we explore the interplay between locality of the response density operator and numerical convergence of LR-related quantities. We show that for frequencies below the first ionization potential (IP) of the system, it is possible to express the response density by employing localized states only. Above this threshold energy, such a locality property cannot be achieved. Such considerations may be transposed in terms of the molecule's excited states. We show that not all the system's excitations can be considered on equal footing. There is a discrete sector of excitations -- which may also extend above IP -- that can be parametrized by observable, localized states, which can be computationally expressed with high precision, provided an adequate level of completeness. We present indicators that can help to quantify such potential observable properties of an excitation, that can be evaluated in any discretization scheme. The remaining excitation modes belong to a continuum spectrum that, on the contrary, is not directly associated to observable properties and can only be effectively represented in a given computational setup. Such considerations are important not only for reproducibility of the results among different computer codes employing diverse formalisms, but also in view of providing a deeper understanding on the impact of models' approximations on the scientific outcomes of the simulation.
△ Less
Submitted 25 February, 2019; v1 submitted 12 July, 2018;
originally announced July 2018.
-
Multiparticle correlation expansion of relative entropy in lattice systems
Authors:
Marco D'Alessandro
Abstract:
This paper deals with the construction of the multiparticle correlation expansion of relative entropy for lattice systems. Thanks to this analysis we are able to express the statistical distance between two systems as a series built over clusters of increasing dimension. Each addend is written in terms of correlation functions and expresses the contribution to the relative entropy due to structura…
▽ More
This paper deals with the construction of the multiparticle correlation expansion of relative entropy for lattice systems. Thanks to this analysis we are able to express the statistical distance between two systems as a series built over clusters of increasing dimension. Each addend is written in terms of correlation functions and expresses the contribution to the relative entropy due to structural information inside the selected cluster. We present a general procedure for the explicit construction of all the terms of the series. As a first application of this result, we show that the coefficients of the multiparticle correlation expansion of the excess entropy can be computed from our formula, as a particular case.
△ Less
Submitted 7 July, 2016;
originally announced July 2016.
-
Maximum-entropy Monte Carlo method for the inversion of the structure factor in simple classical systems
Authors:
Marco D'Alessandro
Abstract:
We present a method for the evaluation of the interaction potential of an equilibrium classical system starting from the (partial) knowledge of its structure factor. The procedure is divided into two phases both of which are based on the maximum entropy principle of information theory. First we determine the maximum entropy estimate of the radial distribution function constrained by the informatio…
▽ More
We present a method for the evaluation of the interaction potential of an equilibrium classical system starting from the (partial) knowledge of its structure factor. The procedure is divided into two phases both of which are based on the maximum entropy principle of information theory. First we determine the maximum entropy estimate of the radial distribution function constrained by the information contained in the structure factor. Next we invert the pair function and extract the interaction potential. The method is tested on a Lennard-Jones fluid at high density and the reliability of its results with respect to the missing information in the structure factor data are discussed. Finally, it is applied to the experimental data of liquid sodium at 100$^{\circ}$C.
△ Less
Submitted 21 June, 2012;
originally announced June 2012.
-
Information-theory-based solution of the inverse problem in classical statistical mechanics
Authors:
Marco D'Alessandro,
Francesco Cilloco
Abstract:
We present a procedure for the determination of the interaction potential from the knowledge of the radial pair distribution function. The method, realized inside an inverse Monte Carlo simulation scheme, is based on the application of the Maximum Entropy Principle of information theory and the interaction potential emerges as the asymptotic expression of the transition probability. Results obtain…
▽ More
We present a procedure for the determination of the interaction potential from the knowledge of the radial pair distribution function. The method, realized inside an inverse Monte Carlo simulation scheme, is based on the application of the Maximum Entropy Principle of information theory and the interaction potential emerges as the asymptotic expression of the transition probability. Results obtained for high density monoatomic fluids are very satisfactory and provide an accurate extraction of the potential, despite a modest computational effort.
△ Less
Submitted 3 September, 2010; v1 submitted 18 June, 2010;
originally announced June 2010.