Impact of the Electronic Band Structure in High-Harmonic Generation Spectra of Solids

Total Page:16

File Type:pdf, Size:1020Kb

Impact of the Electronic Band Structure in High-Harmonic Generation Spectra of Solids Impact of the Electronic Band Structure in High-Harmonic Generation Spectra of Solids The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation Tancogne-Dejean, Nicolas et al. “Impact of the Electronic Band Structure in High-Harmonic Generation Spectra of Solids.” Physical Review Letters 118.8 (2017): n. pag. © 2017 American Physical Society As Published http://dx.doi.org/10.1103/PhysRevLett.118.087403 Publisher American Physical Society Version Final published version Citable link http://hdl.handle.net/1721.1/107908 Terms of Use Article is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use. week ending PRL 118, 087403 (2017) PHYSICAL REVIEW LETTERS 24 FEBRUARY 2017 Impact of the Electronic Band Structure in High-Harmonic Generation Spectra of Solids † Nicolas Tancogne-Dejean,1,2,* Oliver D. Mücke,3,4 Franz X. Kärtner,3,4,5,6 and Angel Rubio1,2,3,5, 1Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany 2European Theoretical Spectroscopy Facility (ETSF), Luruper Chaussee 149, 22761 Hamburg, Germany 3Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestraße 85, 22607 Hamburg, Germany 4The Hamburg Center for Ultrafast Imaging, Luruper Chaussee 149, 22761 Hamburg, Germany 5Physics Department, University of Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany 6Research Laboratory of Electronics, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, USA (Received 29 September 2016; published 24 February 2017) An accurate analytic model describing the microscopic mechanism of high-harmonic generation (HHG) in solids is derived. Extensive first-principles simulations within a time-dependent density-functional framework corroborate the conclusions of the model. Our results reveal that (i) the emitted HHG spectra are highly anisotropic and laser-polarization dependent even for cubic crystals; (ii) the harmonic emission is enhanced by the inhomogeneity of the electron-nuclei potential; the yield is increased for heavier atoms; and (iii) the cutoff photon energy is driver-wavelength independent. Moreover, we show that it is possible to predict the laser polarization for optimal HHG in bulk crystals solely from the knowledge of their electronic band structure. Our results pave the way to better control and optimize HHG in solids by engineering their band structure. DOI: 10.1103/PhysRevLett.118.087403 Atoms and molecules interacting with strong laser pulses The process of HHG from gases is by now well under- emit high-order harmonics of the fundamental driving laser stood in terms of the three-step model [13–15] in which field. The high-harmonic generation (HHG) in gases is electrons are first promoted from the ground state of the routinely used nowadays to produce isolated attosecond atom (or molecule) to the continuum, then accelerated by pulses [1–4] and coherent radiation ranging from the visible the electric field, and finally recombine with the parent ion. to soft x rays [5]. Because of a higher electronic density, With this simple, intuitive model most of the observed solids are one promising route towards compact, brighter effects are well described, in particular, the dependence of HHG sources. The recent observation of nonperturbative the harmonic cutoff energy on driver wavelength and HHG in solids without damage [6–10], extending even intensity. In the case of solids, electrons are promoted to beyond the atomic limit [10], has opened the door to the discrete conduction bands, where they do not evolve freely. observation and control of attosecond electron dynamics in This leads, for instance, to a different linear field dependence solids [8,9,11], all-optical band-structure reconstruction for the cutoff energy [6], different time-frequency character- [12], and solid-state sources of isolated extreme-ultraviolet istics of the harmonic emission between atoms and solids pulses [9,11]. However, in contrast to HHG from gases, the [8,11,16], and a different ellipticity dependence [6,17]. microscopic mechanism underlying HHG from solids is Historically, HHG in solids was first discussed in terms of still controversially debated in the attoscience community, Bloch oscillations (i.e., pure intraband dynamics) [18–20], in some cases casting doubts on the validity of the proposed and more recently mainly analyzed using simplified models microscopic model and resulting in confusion about the based on numerical solutions of the semiconductor Bloch correct interpretation of experimental data. Various com- equations [21–23] treating the complex, coupled interband peting simplified models have been proposed but they often and intraband dynamics, with the exception of the ab initio are based on strong approximations and a priori assump- simulations of Ref. [24]. Even if these methods have been tions, often stating that there is a strong similarity with the successfully applied to some materials, such as GaSe [7,8,25] processes underlying atomic-gas HHG emission. However, or SiO2 [9], basic questions remain controversial and/or it is clear that many-body effects due to the crystalline unresolved, e.g., which bands are involved in the HHG structure of solids and the fermionic nature of interaction dynamics [26,27]. electrons play a decisive role that fundamentally distin- The first experimental observations of HHG from solids guishes the solid from the gas case. It is the scope of the were explained in terms of Bloch oscillations [6,7].A present work to unravel within an ab initio approach what competing model attributing the HHG mechanism to inter- the impact is of the underlying electronic band structure of band transitions (resembling the three-step model of gas the solids in the observed HHG emission. HHG [13,14]) was introduced [28,29]. For not too strong 0031-9007=17=118(8)=087403(6) 087403-1 © 2017 American Physical Society week ending PRL 118, 087403 (2017) PHYSICAL REVIEW LETTERS 24 FEBRUARY 2017 excitation in ZnO by 3.76 μm pulses, such that the electrons formula for the high-harmonic spectra. Using the current explore only the near-parabolic region of the Brillouin zone expressionR for the HHG spectra, namely, HHGðωÞ¼ 3 2 (BZ), it was found that the magnitude of the interband jFTf Ω d rð∂=∂tÞjðr;tÞgj , and plugging now Eq. (3),we contribution is larger than that of the intraband contribution obtain a general expression for the HHG spectra, [28]. Nevertheless, most theoretical works have used either a Z simplified two-band or five-band model, intrinsically ham- 3 HHGðωÞ ∝ FT d r nðr;tÞ∇v0ðrÞ pering the predictive power of the model and the full Ω microscopic understanding of the HHG process. j r;t B r;t 2 n r;t E r;t ð Þ × ð Þ ; In this Letter, using an ab initio approach based on time- þ ð Þ ð Þþ c ð4Þ dependent density-functional theory (TDDFT) [30,31],we study the microscopic origin of HHG in solids. Effects where the last two terms correspond to the Lorentz force stemming from the full electronic structure (valence and exerted by the external laser on the electronic system [34].If conduction bands) and the real crystal structure are properly we now make the dipole approximation, Eq. (4) further accounted for. We show that the nonperturbative emission of simplifies and we finally get harmonics in solids, arising from the interplay between Z 2 intraband and interband contributions, is enhanced when 3 HHGðωÞ ∝ FT d rnðr;tÞ∇v0ðrÞ þ NeEðωÞ ; ð5Þ the interband contribution is suppressed due to band- Ω structure effects. We identify that the joint density of states (JDOS) along the laser polarization is the key parameter which provides the first important physical result of the governing the weight of the interband contribution. In present work, shedding fundamental insights on the intrinsic addition, we address the still controversial question of the bulk contribution to the HHG spectra. Here, Ne is the number wavelength dependence of the cutoff energy in bulk crystals. of electrons contained in the volume Ω. Note that the HHG We start by presenting some exact analytical results (see spectra depend only on the electronic density. The second Supplemental Material [32] for details). We consider a term does not result in a nonperturbative nonlinearity and thus general interacting many-electron Hamiltonian Hˆ of the form cannot create a plateaulike HHG spectrum. The more ˆ ˆ ˆ ˆ interesting and relevant term for HHG is the first one in HðtÞ¼T þ VðtÞþW; ð1Þ Eq. (5). It shows that higher harmonics are generated by two Tˆ Vˆ t competing terms, the spatial variation of the total electronic where is the kinetic energy, ð Þ is the time-dependent density [nðr;tÞ] and the gradient of the electron-nuclei Wˆ external laser potential, and is the electron-electron potential [∇v0ðrÞ], the latter being time independent, as Coulomb interaction (the ionic motion is not considered we neglected ionic motion. In gases, the gradient of the here for the sake of simplicity). The exact equation of motion electron-nuclei potential is important, but the electronic for the total microscopic current, jðr;tÞ, can be rewritten density is low. In the case of solids, the electronic density as [33,34] is higher, but the potential is rather homogeneous, resulting in ∂ a smaller gradient of the potential than in the atomic case. In j r;t −n r;t ∇v r;t Πkin r;t Πint r;t ; fact, in the limit of a homogeneous electron gas, the gradient ∂t ð Þ¼ ð Þ ð Þþ ð Þþ ð Þ ð2Þ becomes 0, and no harmonics are generated, irrespective of where Πkinðr;tÞ and Πintðr;tÞ are the kinetic and the the value of the electronic density.
Recommended publications
  • Generation of Even and Odd High Harmonics in Resonant Metasurfaces Using Single and Multiple Ultra-Intense Laser Pulses
    Generation of even and odd high harmonics in resonant metasurfaces using single and multiple ultra-intense laser pulses Maxim Shcherbakov ( [email protected] ) Cornell University Haizhong Zhang Institute of Materials Research and Engineering, A*STAR (Agency for Science, Technology and Research) Michael Tripepi Department of Physics, The Ohio State University https://orcid.org/0000-0002-4657-9826 Giovanni Sartorello Cornell University Noah Talisa Department of Physics, The Ohio State University Abdallah AlShafey Department of Physics, The Ohio State University Zhiyuan Fan Cornell University Justin Twardowski Department of Physics, The Ohio State University Leonid Krivitsky Institute of Materials Research and Engineering, A*STAR (Agency for Science, Technology and Research) Arseniy Kuznetsov Institute of Materials Research and Engineering https://orcid.org/0000-0002-7622-8939 Enam Chowdhury Department of Physics, The Ohio State University Gennady Shvets Cornell University https://orcid.org/0000-0001-8542-6753 Article Keywords: High harmonic generation (HHG), resonant metasurfaces, ultra-intense laser pulses, Posted Date: October 26th, 2020 DOI: https://doi.org/10.21203/rs.3.rs-92835/v1 License: This work is licensed under a Creative Commons Attribution 4.0 International License. Read Full License Generation of even and odd high harmonics in resonant metasurfaces using single and multiple ultra-intense laser pulses Maxim R. Shcherbakov,1,* Haizhong Zhang,2 Michael Tripepi,3 Giovanni Sartorello,1 Noah Talisa,3 Abdallah AlShafey,3 Zhiyuan Fan,1 Justin Twardowski,4 Leonid A. Krivitsky,2 Arseniy I. Kuznetsov,2 Enam Chowdhury,3,4,5 Gennady Shvets1,* Affiliations: 1School of Applied and Engineering Physics, Cornell University, Ithaca, NY 14853, USA. 2Institute of Materials Research and Engineering, A*STAR (Agency for Science, Technology and Research), 138634, Singapore.
    [Show full text]
  • Light-Matter Interaction and Optical Spectroscopy from Infrared to X-Rays
    24th ETSF Workshop on Electronic Excitations Light-Matter Interaction and Optical Spectroscopy from Infrared to X-Rays Jena, Germany 16 – 20 September 2019 Welcome The workshop series of the European Theoretical Spectroscopy Facility (ETSF) provides a fo- rum for excited states and spectroscopy in condensed-matter physics, chemistry, nanoscience, materials science, and molecular physics attracting theoreticians, code developers, and experi- mentalists alike. Light-matter interaction will be at core of the 2019 edition of the ETSF workshop. Cutting-edge spectroscopy experiments allow to probe electrons, plasmons, excitons, and phonons across different energy and time scales with unprecedented accuracy. A deep physical understanding of the underlying quantum many-body effects is of paramount importance to analyze exper- imental observations and render theoretical simulations predictive. The workshop aims at discussing the most recent advances in the theoretical description of the interaction between light and matter focusing on first-principles methods. This broad subject will be covered in its diverse declinations, from core-level spectroscopy to collective low-energy excitations, dis- cussing also matter under extreme conditions, and systems driven out of equilibrium by strong laser pulses. Exchange between theorists and experimentalists is fostered to open new horizons towards the next generation of novel spectroscopy techniques. The workshop will also face the challenges posed by the formidable complexity of heterogeneous and nanostructured systems such as those of interest for light harvesting and energy generation, prompting to bridge the gap be- tween experimental and in silico spectroscopy. Workshop topics include: Linear and nonlinear optical spectroscopy • Core-level spectroscopies • Ultrafast excitation dynamics • Electron-phonon coupling • Light harvesting in natural and synthetic systems • We are glad to welcome you in Jena, the city of light, and wish you an inspiring workshop with lots of interesting science and fruitful discussions.
    [Show full text]
  • High Harmonic Generation by Short Laser Pulses
    Louisiana State University LSU Digital Commons LSU Doctoral Dissertations Graduate School 2006 High harmonic generation by short laser pulses: time-frequency behavior and applications to attophysics Mitsuko Murakami Louisiana State University and Agricultural and Mechanical College Follow this and additional works at: https://digitalcommons.lsu.edu/gradschool_dissertations Part of the Physical Sciences and Mathematics Commons Recommended Citation Murakami, Mitsuko, "High harmonic generation by short laser pulses: time-frequency behavior and applications to attophysics" (2006). LSU Doctoral Dissertations. 769. https://digitalcommons.lsu.edu/gradschool_dissertations/769 This Dissertation is brought to you for free and open access by the Graduate School at LSU Digital Commons. It has been accepted for inclusion in LSU Doctoral Dissertations by an authorized graduate school editor of LSU Digital Commons. For more information, please [email protected]. HIGH HARMONIC GENERATION BY SHORT LASER PULSES: TIME-FREQUENCY BEHAVIOR AND APPLICATIONS TO ATTOPHYSICS A Dissertation Submitted to the Graduate Faculty of the Louisiana State University and Agricultural and Mechanical College in partial fulfillment of the requirements for the degree of Doctor of Philosophy in The Department of Physics and Astronomy by Mitsuko Murakami B. Sc., McNeese State University, Lake Charles, LA, 1999 M. Sc., Louisiana State University, Baton Rouge, LA, 2001 May, 2006 To my grandmother ii Acknowledgements I am very grateful to my advisor, Dr. Mette Gaarde, who has been helpful and inspiring throughout my dissertation research. Her teaching is rigorous upon fundamental principles and yet extremely intuitive. I am also grateful to Dr. Ravi Rau, Dr. Dana Browne, Dr. Gabriela Gonzalez, and Dr. Xue-Bin Liang for serving as the members of my dissertation committee.
    [Show full text]
  • Programming Models for Quantum Many-Body Methods on Multicore and Manycore Processors
    Programming models for quantum many-body methods on multicore and manycore processors Jeff Hammond1 and Eugene DePrince2 1 Argonne 2 Georgia Tech 6 February 2011 Jeff Hammond Electronic Structure Calculation Methods on Accelerators Abstract The growing diversity in computer processor architectures poses a serious challenge to the computational chemistry community. This talk considers some of the key issues, including disjoint address spaces, non-standard architectures and execution models, and the different APIs required to use them. Specifically, we will describe our experiences in developing coupled-cluster methods for Intel multicore, Intel MIC, NVIDIA Fermi and Blue Gene/Q in both a clean-sheet implementation and NWChem. Of fundamental interest is the development of codes that scale not only within the node but across thousands of nodes; hence, the interaction between the processor and the network will be analyzed in detail. Jeff Hammond Electronic Structure Calculation Methods on Accelerators Software Automation Jeff Hammond Electronic Structure Calculation Methods on Accelerators The practical TCE { NWChem many-body codes What does it do? 1 GUI input quantum many-body theory e.g. CCSD. 2 Operator specification of theory. 3 Apply Wick's theory to transform operator expressions into array expressions. 4 Transform input array expression to operation tree using many types of optimization. 5 Produce Fortran+Global Arrays+NXTVAL implementation. Developer can intercept at various stages to modify theory, algorithm or implementation. Jeff Hammond Electronic Structure Calculation Methods on Accelerators The practical TCE { Success stories First parallel implementation of many (most) CC methods. First truly generic CC code (not string-based): fRHF,ROHF,UHF}×CCfSD,SDT,SDTQ}×{T /Λ,EOM,LR/QRg Most of the largest calculations of their kind employ TCE: CR-EOMCCSD(T), CCSD-LR α, CCSD-QR β, CCSDT-LR α Reduces implementation time for new methods from years to hours, TCE codes are easy to verify.
    [Show full text]
  • Propriétés Structurales, Opto-Électroniques, Élastiques Et Dynamiques Des Semi-Conducteurs Type II-VI
    REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE MINISTERE DE L’ENSEIGNEMENT SUPERIEUR ET DE LA RECHERCHE SCIENTIFIQUE UNIVERSITE FERHAT ABBAS – SETIF THESE Présentée à la faculté des Sciences Département de Physique Pour l’obtention du diplôme de DOCTORAT EN SCIENCES Option : physique du solide Par Mr. Benamrani Ammar Thème Propriétés structurales, opto-électroniques, élastiques et dynamiques des Semi-conducteurs type II-VI Soutenue le : 02/06/2012 Devant la commission d’examen : Dr. L. Louail Prof. U. F. A. SETIF Président Dr. K. Kassali Prof. U. F. A. SETIF Rapporteur Dr. Kh. Bouamama Prof. U. F. A. SETIF Co-Encadreur Dr. B. Bennecer Prof. U. Guelma Examinateur Dr. A. El Akrmi Prof. U. Annaba Examinateur Dr. H. Belkhir Prof. U. Annaba Examinateur À la mémoire de mon père, À ma mère, qui m’a enseigné le sens de la patience et de l’espoir, À ma femme, À toute ma famille À toute personne qui explore honnêtement les secrets de ce monde vaste, Je dédie ce modeste travail Remerciements REMERCIEMENTS Ce travail a été réalisé au laboratoire d’optoélectronique et des composants à l’université Ferhat Abbes à Sétif. Je tiens à remercier profondément Pr. Kamel Kassali, mon directeur de thèse pour m’avoir proposé ce thème intéressant et pour sa patience durant toutes les années de préparation de la thèse ainsi que ses conseils précieux et ses réponses à toutes mes questions reliées à mon sujet de thèse. J’exprime également toute ma reconnaissance à Mr. Kh. Bouamama, Professeur à l’université Ferhat Abbes de Sétif, qui m’a grandement fait profiter de ses fructueuses discussions en codes de calcul et dans mon sujet de thèse.
    [Show full text]
  • DFT-Absorption Spectra
    2054-9 Structure and Dynamics of Hydrogen-Bonded Systems 26 - 27 October 2009 Effect of Proton Disorder on the Excited State Properties of Ice Olivia PULCI Universita' di Roma II "Tor Vergata" Dipt. di Fisica, Via della Ricerca Scientfica 1, 00133 Rome Italy EffectEffect ofof protonproton disorderdisorder onon thethe excitedexcited statestate propertiesproperties ofof iceice V. Garbuio,Garbuio M. Cascella, R. Del Sole, O. Pulci OUTLINE: •Theoretical approaches •Ice Ic(bulk) • Ice Ih surface Theoretical approaches c MBPT c c EXC W hν hν ωcv hν ωcv ωcv v v v DFT GW BSE 1) 2) 3) Ground state properties Electronic band Optical spectra structure, I, A TDDFT (Step 2) Lars Hedin 1965 Σ = iGW G: single particle Green’s function −1 W: screened Coulomb interaction = ε VW Theoretical approaches c MBPT c c EXC W hν hν ωcv hν ωcv ωcv v v v DFT GW BSE 1) 2) 3) Ground state properties Electronic band Optical spectra structure, I, A TDDFT Step 3: calculation of optical spectra within the Bethe Salpeter Equation c Absorption spectra A photon excites an electron from an occupied state to a conduction state hν e 4 4 4 4 4 PPPP=IQP +IQP Ξ v h Bethe Salpeter Equation (BSE) Kernel: Ξ v= − W e-h exchange bound excitons GWBSE ApplicableApplicable to:to: Ab-initio: (NOT “one puts nothing in, one gets nothing out”!!) •Generality, transferability 0D-3D Biological systems •Detailed physical informations •Complex theory+large comp.cost 3-D 0-D 1-D 2-D Nanowires Surfaces Nanoclusters bulks OUTLINE: •Theoretical approaches •Ic Ice (bulk) • Ice Ih surface HH2OO phasephase diagramdiagram CubicCubic iceice ((IcIc)) • It is a metastable form of ice that can be formed, by condensation of water vapor, at ambient pressure but low temperatures Cubic ice (Ic) – diamond lattice 153 K down to 113 K Studied within DFT and Tight-binding G.
    [Show full text]
  • Enhanced High-Harmonic Generation from an All-Dielectric Metasurface
    LETTERS https://doi.org/10.1038/s41567-018-0233-6 Enhanced high-harmonic generation from an all-dielectric metasurface Hanzhe Liu 1,2*, Cheng Guo3,4, Giulio Vampa1, Jingyuan Linda Zhang3,4, Tomas Sarmiento 4,5, Meng Xiao 4, Philip H. Bucksbaum1,2,3,6, Jelena Vucković3,4,5 , Shanhui Fan3,4,5 and David A. Reis 1,3,6* The recent observation of high-harmonic generation from damage threshold of plasmonic-based metal structures is extremely solids creates a new possibility for engineering fundamen- sensitive to the fabrication quality and the high degree of confine- tal strong-field processes by patterning the solid target with ment in the vicinity of sharp features, such as edges or corners, sig- subwavelength nanostructures. All-dielectric metasurfaces nificantly limits the overall signal enhancement12,18,19. exhibit high damage thresholds and strong enhancement of A dielectric metasurface that hosts an optical resonance provides the driving field, making them attractive platforms to con- a promising solution to this issue. Electromagnetically induced trol high harmonics and other high-field processes at the transparency (EIT)20 can give rise to greatly enhanced nonlinear nanoscale. Here we report enhanced non-perturbative high- light–matter interactions in the spectral region of the resonance21. harmonic emission from a Fano-resonant Si metasurface Enhanced frequency conversion from EIT has been observed in that possesses a classical analogue of electromagnetically atomic systems21 and, very recently, enhanced perturbative third- induced transparency. The harmonic emission is enhanced by harmonic generation has been reported from photonic devices more than two orders of magnitude compared to unpatterned exhibiting a classical analogue of EIT14.
    [Show full text]
  • X-Ray Spectroscopy and Metrology with High-Harmonic Sources
    X-ray spectroscopy and metrology with high-harmonic sources Peter M. Kraus Advanced Research Center for Nanolithography (ARCNL), Amsterdam, The Netherlands. EUV Lithography Source Workshop, 06.11.2018 – HiLASE, Prague, Czech Republic. EUV lithography requires new XUV sources Resolve structures during lithography and Understand EUV (13.5 nm) interactions to after production with sub-nm precision improve masks, pellicles, mirrors, resists… L. Li et al., Chem. Soc. Rev., 2017, 46, 4855 Y. Zhang et al., J. Micro Nanolith. MEMS MOEMS 16, 023510 (2017) Broadband XUV and soft x-ray imaging Ulrafast XUV and soft x-ray spectroscopy We need an ultrafast, coherent, table-top XUV/soft x-ray source Solution: High-harmonic generation (HHG) High-harmonic generation argon atoms 1300 nm Space ~ 1 mJ ~1014 W/cm2 20 eV 50 eV 120 eV =^ 60 nm =^ 10 nm High-harmonic generation (HHG): Converting many IR-photons into one 3 XUV-photon Three-step model of HHG Typical parameters: E - 3 - Recombination 800 nm 3.17 Up - 2 - ~ 1 mJ Propagation ~ 1014 - 1015 W/cm2 XUV emission 0 r Active Figure from: electron P. M. Kraus, H. J. Woerner, Angew. Chem. 57, - IP 5228 (2018). - 1 - Three-step model of HHG: Core P. B. Corkum, PRL 71, 1994 (1993) Ionization J. L. Krause, K. J. Schafer, K. C. Kulander, PRL 68, 3535 (1992) Quantum theory of HHG: Electron M. Lewenstein et al., PRA 49, 2117 (1994) mean field Imaging needs in nanolithpgraphy Lensless imaging with x-rays from HHG Compared to a conventional microscope, the lenses are replaced with a computer algorithm Collaboration with Stefan Witte & Kjeld Eikema (ARCNL) Needs in lithography: Short wavelength, high average power Demonstrated HHG sources S.
    [Show full text]
  • New Frontiers in Quantum Chemistry Using Supercomputers
    New frontiers in quantum chemistry using supercomputers Jeff Hammond Leadership Computing Facility Argonne National Laboratory 4 May 2011 Jeff Hammond NREL Atomistic simulation in chemistry 1 classical molecular dynamics (MD) with empirical potentials 2 ab initio molecular dynamics based upon density-function theory (DFT) 3 quantum chemistry with wavefunctions e.g. perturbation theory (PT), Coupled-Cluster (CC) or Quantum Monte Carlo (QMC). Jeff Hammond NREL Classical molecular dynamics Solves Newton's equations of motion with empirical terms and classical electrostatics. Size: 100K-10M atoms Time: 1-10 ns/day Scaling: ∼ Natoms Data from K. Schulten, et al. \Biomolecular modeling in the era of petascale computing." In D. Bader, ed., Petascale Computing: Algorithms and Applications. Image courtesy of Beno^ıtRoux via ALCF. Jeff Hammond NREL Car-Parrinello molecular dynamics Forces obtained from solving an approximate single-particle Schr¨odingerequation; time-propagation via Lagrangian approach. Size: 100-1000 atoms Time: 0.01-1 ps/day x Scaling: ∼ Nel (x=1-3) F. Gygi, IBM J. Res. Dev. 52, 137 (2008); E. J. Bylaska et al. J. Phys.: Conf. Ser. 180, 012028 (2009). Image courtesy of Giulia Galli via ALCF. Jeff Hammond NREL Wavefunction theory MP2 is second-order PT and is accurate via magical cancellation of error. CC is infinite-order solution to many-body Schr¨odingerequation truncated via clusters. QMC is Monte Carlo integration applied to the Schr¨odingerequation. Size: 10-100 atoms, maybe 100-1000 atoms with MP2. Time: N/A x Scaling: ∼ Nbf (x=2-7) Image courtesy of Karol Kowalski and Niri Govind. Jeff Hammond NREL Quantum chemistry | standard model 1 Separate molecule(s) from environment (closed to both matter and energy) 2 Ignore relativity, QED, spin-orbit coupling 3 Separate electronic and nuclear degrees of freedom −! non-relativistic electronic Schr¨odingerequation in a vacuum at zero temperature.
    [Show full text]
  • Spin Conservation in High-Order-Harmonic Generation Using Bicircular fields
    Spin conservation in high-order-harmonic generation using bicircular fields Emilio Pisanty1,∗ Suren Sukiasyan1, and Misha Ivanov1;2;3† 1Blackett Laboratory, Imperial College London, South Kensington Campus, SW7 2AZ London, United Kingdom 2Department of Physics, Humboldt University, Newtonstrasse 15, 12489 Berlin, Germany 3Max Born Institute, Max Born Strasse 2a, 12489 Berlin, Germany (Dated: October 3, 2018) We present an alternative theoretical model for a recent experiment [A. Fleischer et al., Nature Photon. 8, 543 (2014)] which used bichromatic, counter-rotating high intensity laser pulses to probe the conservation of spin angular momentum in high harmonic generation. We separate elliptical polarizations into independent circular fields with definite angular momentum, instead of using the expectation value of spin for each photon in the conservation equation, and we find good agreement with the experimental results. In our description the generation of each individual harmonic con- serves spin angular momentum, in contrast to the model proposed by Fleischer et al. Our model also correctly describes analogous processes in standard perturbative optics. §I. INTRODUCTION that in general the production of each individual harmonic is not a closed process. The process of high harmonic generation [1] is the flag- This model, which we call Model 1, suggests that spin ship experiment of extreme nonlinear optics. It consists of angular momentum is not conserved for certain harmonic irradiating atoms or small molecules with long-wavelength lines, so that some harmonics must be assumed to be emit- laser pulses whose electric field is comparable to the inter- ted in correlated pairs for the overall process to be paramet- nal electric fields of the atoms, which results in the emission ric.
    [Show full text]
  • High Harmonic Generation Physics 208A Presentation October 11, 2004
    High Harmonic Generation Physics 208A Presentation October 11, 2004 Tom Allison UC Berkeley Perturbative Nonlinear Optics P = ! (1)E + ! (2)E 2 + ! (3)E 3 + ... •Accurately treated by treating the polarization as a power series in E. •With sufficiently intense laser fields, the higher order terms give rise to Fourier components of the polarization at harmonics of the laser frequency, creating radiation at harmonics of the laser frequency. • First demonstration of second harmonic generation by P.A. Franken et al. (1961). • Laser intensity of ~1×107 W/cm2 • E ~105 V/cm = 10-3 V/Angstrom. • For reference, the atomic field is e Eat ! 2 ! 50 V/Angstrom a0 • With careful phase matching, generation of low order harmonics can be very efficient! XUV or X-rays? • ALS costs ~$100,000/day in operating costs. • Large facility shared by many users ⇒ beamtime is scarce. • X-ray pulse durations typically ~100 picoseconds - too slow to resolve ultrafast dynamics. • Can we use ~100th harmonic generation as a table-top, high brightness, ultrafast x-ray source? How can we extend this process to higher harmonics? • Most obvious first try would be to successively double: How can we extend this process to higher harmonics? • Most obvious first try would be to successively double: UV absorption in doubling crystals?! Phase matching in region of strong dispersion?! High Order Harmonics • In order to get harmonics of very high order, we need to go to extremely intense electric fields, and use gas density targets. Harmonic Spectrum Harmonic spectrum is flat! Sharp cutoff! Old odd harmonics appear! Not well described By expansion of the Polarization! Ferray et al.
    [Show full text]
  • 0.1 Report on the 12Th Nanoquanta Workshop on Electronic Excitations
    0.1 Report on the 12th Nanoquanta Workshop on Electronic Excitations. Time-Dependent Density-Functional Theory: Advances and Prospects Aussois (France) 18-22 September 2007 SPONSORS NANOQUANTA Network of Excellence ESF Psi k Programme CNRS IdNano ORGANIZERS Valerio OLEVANO (LEPES-CNRS, Grenoble, France) John J. REHR (University of Washington, Seattle, USA) Gian-Marco RIGNANESE (Universit´eCatholique de Louvain, Louvain-la-Neuve, Belgium) Patrick RINKE (Fritz-Haber-Institut, Berlin, Germany) Francesco SOTTILE (Universidad del Pais Vasco, San Sebasti´an,Spain) Ludger WIRTZ (IEMN/ISEN, Villeneuve d’Ascq, France) WEBSITE http://lab-neel.grenoble.cnrs.fr/etsf/nanoquanta-workshop07 1 The workshop gathered 120 participants from leading international theory groups. There were 46 oral presen- tations (23 invited speakers and 23 contributed talks) and 32 posters. The workshop allowed the participants to discuss the advances in the theoretical and computational treatment of optical and dielectric spectroscopy, photoemission spectroscopy and quantum transport in the framework of time-dependent density-functional theory (TDDFT), many-body perturbation theory (MBPT) and non-equilibrium Green’s function (NEGF) theory. Invited speakers from leading international groups gave an in-depth overview of current research activities within these theories and placed recent results into context. Young researchers (Ph.D. students and Post- Docs) also had the opportunity to present their work (more than half of the oral presentations were given by non-permanent researchers). A special emphasis was placed on advances and perspectives in time-dependent density-functional theory. Indeed, the first session presented two review talks on TDDFT and DFT-like methods applied to quantum transport. A particular effort has been done by the two invited speakers to present the arguments in the most pedagogical manner, ad usum of the youngest participants to the workshop.
    [Show full text]