@article{18466, abstract = {{The transverse dynamic spin susceptibility is a correlation function that yields exact information about spin excitations in systems with a collinear magnetic ground state, including collective spin-wave modes. In an ab initio context, it may be calculated within many-body perturbation theory or time-dependent density-functional theory, but the quantitative accuracy is currently limited by the available functionals for exchange and correlation in dynamically evolving systems. To circumvent this limitation, the spin susceptibility is here alternatively formulated as the solution of an initial-value problem. In this way, the challenge of accurately describing exchange and correlation in many-electron systems is shifted to the stationary initial state, which is much better understood. The proposed scheme further requires the choice of an auxiliary basis set, which determines the speed of convergence but always allows systematic convergence in practical implementations.}}, author = {{Schindlmayr, Arno}}, issn = {{1687-9139}}, journal = {{Advances in Mathematical Physics}}, publisher = {{Hindawi}}, title = {{{Exact formulation of the transverse dynamic spin susceptibility as an initial-value problem}}}, doi = {{10.1155/2018/3732892}}, volume = {{2018}}, year = {{2018}}, } @article{10023, abstract = {{We perform a comprehensive theoretical study of the structural and electronic properties of potassium niobate (KNbO3) in the cubic, tetragonal, orthorhombic, monoclinic, and rhombohedral phase, based on density-functional theory. The influence of different parametrizations of the exchange-correlation functional on the investigated properties is analyzed in detail, and the results are compared to available experimental data. We argue that the PBEsol and AM05 generalized gradient approximations as well as the RTPSS meta-generalized gradient approximation yield consistently accurate structural data for both the external and internal degrees of freedom and are overall superior to the local-density approximation or other conventional generalized gradient approximations for the structural characterization of KNbO3. Band-structure calculations using a HSE-type hybrid functional further indicate significant near degeneracies of band-edge states in all phases which are expected to be relevant for the optical response of the material.}}, author = {{Schmidt, Falko and Landmann, Marc and Rauls, Eva and Argiolas, Nicola and Sanna, Simone and Schmidt, Wolf Gero and Schindlmayr, Arno}}, issn = {{1687-8442}}, journal = {{Advances in Materials Science and Engineering}}, publisher = {{Hindawi}}, title = {{{Consistent atomic geometries and electronic structure of five phases of potassium niobate from density-functional theory}}}, doi = {{10.1155/2017/3981317}}, volume = {{2017}}, year = {{2017}}, } @article{10021, abstract = {{The optical properties of pristine and titanium-doped LiNbO3 are modeled from first principles. The dielectric functions are calculated within time-dependent density-functional theory, and a model long-range contribution is employed for the exchange-correlation kernel in order to account for the electron-hole binding. Our study focuses on the influence of substitutional titanium atoms on lithium sites. We show that an increasing titanium concentration enhances the values of the refractive indices and the reflectivity.}}, author = {{Friedrich, Michael and Schmidt, Wolf Gero and Schindlmayr, Arno and Sanna, Simone}}, issn = {{2475-9953}}, journal = {{Physical Review Materials}}, number = {{3}}, publisher = {{American Physical Society}}, title = {{{Optical properties of titanium-doped lithium niobate from time-dependent density-functional theory}}}, doi = {{10.1103/PhysRevMaterials.1.034401}}, volume = {{1}}, year = {{2017}}, } @article{13416, abstract = {{The optical properties of congruent lithium niobate are analyzed from first principles. The dielectric function of the material is calculated within time-dependent density-functional theory. The effects of isolated intrinsic defects and defect pairs, including the NbLi4+ antisite and the NbLi4+−NbNb4+ pair, commonly addressed as a bound polaron and bipolaron, respectively, are discussed in detail. In addition, we present further possible realizations of polaronic and bipolaronic systems. The absorption feature around 1.64 eV, ascribed to small bound polarons [O. F. Schirmer et al., J. Phys.: Condens. Matter 21, 123201 (2009)], is nicely reproduced within these models. Among the investigated defects, we find that the presence of bipolarons at bound interstitial-vacancy pairs NbV−VLi can best explain the experimentally observed broad absorption band at 2.5 eV. Our results provide a microscopic model for the observed optical spectra and suggest that, besides NbLi antisites and Nb and Li vacancies, Nb interstitials are also formed in congruent lithium-niobate samples.}}, author = {{Friedrich, Michael and Schmidt, Wolf Gero and Schindlmayr, Arno and Sanna, Simone}}, issn = {{2475-9953}}, journal = {{Physical Review Materials}}, number = {{5}}, publisher = {{American Physical Society}}, title = {{{Polaron optical absorption in congruent lithium niobate from time-dependent density-functional theory}}}, doi = {{10.1103/PhysRevMaterials.1.054406}}, volume = {{1}}, year = {{2017}}, } @article{7481, abstract = {{The electronic band structures of hexagonal ZnO and cubic ZnS, ZnSe, and ZnTe compounds are determined within hybrid-density-functional theory and quasiparticle calculations. It is found that the band-edge energies calculated on the G0W0 (Zn chalcogenides) or GW (ZnO) level of theory agree well with experiment, while fully self-consistent QSGW calculations are required for the correct description of the Zn 3d bands. The quasiparticle band structures are used to calculate the linear response and second-harmonic-generation (SHG) spectra of the Zn–VI compounds. Excitonic effects in the optical absorption are accounted for within the Bethe–Salpeter approach. The calculated spectra are discussed in the context of previous experimental data and present SHG measurements for ZnO.}}, author = {{Riefer, Arthur and Weber, Nils and Mund, Johannes and Yakovlev, Dmitri R. and Bayer, Manfred and Schindlmayr, Arno and Meier, Cedrik and Schmidt, Wolf Gero}}, issn = {{1361-648X}}, journal = {{Journal of Physics: Condensed Matter}}, number = {{21}}, publisher = {{IOP Publishing}}, title = {{{Zn–VI quasiparticle gaps and optical spectra from many-body calculations}}}, doi = {{10.1088/1361-648x/aa6b2a}}, volume = {{29}}, year = {{2017}}, } @article{10024, abstract = {{The influence of electronic many-body interactions, spin-orbit coupling, and thermal lattice vibrations on the electronic structure of lithium niobate is calculated from first principles. Self-energy calculations in the GW approximation show that the inclusion of self-consistency in the Green function G and the screened Coulomb potential W opens the band gap far stronger than found in previous G0W0 calculations but slightly overestimates its actual value due to the neglect of excitonic effects in W. A realistic frozen-lattice band gap of about 5.9 eV is obtained by combining hybrid density functional theory with the QSGW0 scheme. The renormalization of the band gap due to electron-phonon coupling, derived here using molecular dynamics as well as density functional perturbation theory, reduces this value by about 0.5 eV at room temperature. Spin-orbit coupling does not noticeably modify the fundamental gap but gives rise to a Rashba-like spin texture in the conduction band.}}, author = {{Riefer, Arthur and Friedrich, Michael and Sanna, Simone and Gerstmann, Uwe and Schindlmayr, Arno and Schmidt, Wolf Gero}}, issn = {{2469-9969}}, journal = {{Physical Review B}}, number = {{7}}, publisher = {{American Physical Society}}, title = {{{LiNbO3 electronic structure: Many-body interactions, spin-orbit coupling, and thermal effects}}}, doi = {{10.1103/PhysRevB.93.075205}}, volume = {{93}}, year = {{2016}}, } @article{10025, abstract = {{The phonon dispersions of the ferro‐ and paraelectric phase of LiTaO3 are calculated within density‐functional perturbation theory. The longitudinal optical phonon modes are theoretically derived and compared with available experimental data. Our results confirm the recent phonon assignment proposed by Margueron et al. [J. Appl. Phys. 111, 104105 (2012)] on the basis of spectroscopical studies. A comparison with the phonon band structure of the related material LiNbO3 shows minor differences that can be traced to the atomic‐mass difference between Ta and Nb. The presence of phonons with imaginary frequencies for the paraelectric phase suggests that it does not correspond to a minimum energy structure, and is compatible with an order‐disorder type phase transition.}}, author = {{Friedrich, Michael and Schindlmayr, Arno and Schmidt, Wolf Gero and Sanna, Simone}}, issn = {{1521-3951}}, journal = {{Physica Status Solidi B}}, number = {{4}}, pages = {{683--689}}, publisher = {{Wiley-VCH}}, title = {{{LiTaO3 phonon dispersion and ferroelectric transition calculated from first principles}}}, doi = {{10.1002/pssb.201552576}}, volume = {{253}}, year = {{2016}}, } @article{10030, abstract = {{The vibrational properties of stoichiometric LiNbO3 are analyzed within density-functional perturbation theory in order to obtain the complete phonon dispersion of the material. The phonon density of states of the ferroelectric (paraelectric) phase shows two (one) distinct band gaps separating the high-frequency (~800 cm−1) optical branches from the continuum of acoustic and lower optical phonon states. This result leads to specific heat capacites in close agreement with experimental measurements in the range 0–350 K and a Debye temperature of 574 K. The calculated zero-point renormalization of the electronic Kohn–Sham eigenvalues reveals a strong dependence on the phonon wave vectors, especially near Γ. Integrated over all phonon modes, our results indicate a vibrational correction of the electronic band gap of 0.41 eV at 0 K, which is in excellent agreement with the extrapolated temperature-dependent measurements.}}, author = {{Friedrich, Michael and Riefer, Arthur and Sanna, Simone and Schmidt, Wolf Gero and Schindlmayr, Arno}}, issn = {{1361-648X}}, journal = {{Journal of Physics: Condensed Matter}}, number = {{38}}, publisher = {{IOP Publishing}}, title = {{{Phonon dispersion and zero-point renormalization of LiNbO3 from density-functional perturbation theory}}}, doi = {{10.1088/0953-8984/27/38/385402}}, volume = {{27}}, year = {{2015}}, } @article{18470, abstract = {{Using ab initio computational methods, we study the structural and electronic properties of strained silicon, which has emerged as a promising technology to improve the performance of silicon-based metal-oxide-semiconductor field-effect transistors. In particular, higher electron mobilities are observed in n-doped samples with monoclinic strain along the [110] direction, and experimental evidence relates this to changes in the effective mass as well as the scattering rates. To assess the relative importance of these two factors, we combine density-functional theory in the local-density approximation with the GW approximation for the electronic self-energy and investigate the effect of uniaxial and biaxial strains along the [110] direction on the structural and electronic properties of Si. Longitudinal and transverse components of the electron effective mass as a function of the strain are derived from fits to the quasiparticle band structure and a diagonalization of the full effective-mass tensor. The changes in the effective masses and the energy splitting of the conduction-band valleys for uniaxial and biaxial strains as well as their impact on the electron mobility are analyzed. The self-energy corrections within GW lead to band gaps in excellent agreement with experimental measurements and slightly larger effective masses than in the local-density approximation.}}, author = {{Bouhassoune, Mohammed and Schindlmayr, Arno}}, issn = {{1687-8124}}, journal = {{Advances in Condensed Matter Physics}}, publisher = {{Hindawi}}, title = {{{Ab initio study of strain effects on the quasiparticle bands and effective masses in silicon}}}, doi = {{10.1155/2015/453125}}, volume = {{2015}}, year = {{2015}}, } @inbook{18471, abstract = {{Collective spin excitations form a fundamental class of excitations in magnetic materials. As their energy reaches down to only a few meV, they are present at all temperatures and substantially influence the properties of magnetic systems. To study the spin excitations in solids from first principles, we have developed a computational scheme based on many-body perturbation theory within the full-potential linearized augmented plane-wave (FLAPW) method. The main quantity of interest is the dynamical transverse spin susceptibility or magnetic response function, from which magnetic excitations, including single-particle spin-flip Stoner excitations and collective spin-wave modes as well as their lifetimes, can be obtained. In order to describe spin waves we include appropriate vertex corrections in the form of a multiple-scattering T matrix, which describes the coupling of electrons and holes with different spins. The electron–hole interaction incorporates the screening of the many-body system within the random-phase approximation. To reduce the numerical cost in evaluating the four-point T matrix, we exploit a transformation to maximally localized Wannier functions that takes advantage of the short spatial range of electronic correlation in the partially filled d or f orbitals of magnetic materials. The theory and the implementation are discussed in detail. In particular, we show how the magnetic response function can be evaluated for arbitrary k points. This enables the calculation of smooth dispersion curves, allowing one to study fine details in the k dependence of the spin-wave spectra. We also demonstrate how spatial and time-reversal symmetry can be exploited to accelerate substantially the computation of the four-point quantities. As an illustration, we present spin-wave spectra and dispersions for the elementary ferromagnet bcc Fe, B2-type tetragonal FeCo, and CrO2 calculated with our scheme. The results are in good agreement with available experimental data.}}, author = {{Friedrich, Christoph and Şaşıoğlu, Ersoy and Müller, Mathias and Schindlmayr, Arno and Blügel, Stefan}}, booktitle = {{First Principles Approaches to Spectroscopic Properties of Complex Materials}}, editor = {{Di Valentin, Cristiana and Botti, Silvana and Cococcioni, Matteo}}, isbn = {{978-3-642-55067-6}}, issn = {{1436-5049}}, pages = {{259--301}}, publisher = {{Springer}}, title = {{{Spin excitations in solids from many-body perturbation theory}}}, doi = {{10.1007/128_2013_518}}, volume = {{347}}, year = {{2014}}, } @inbook{18472, abstract = {{Many-body perturbation theory is a well-established ab initio electronic-structure method based on Green functions. Although computationally more demanding than density functional theory, it has the distinct advantage that the exact expressions for all relevant observables, including the ground-state total energy, in terms of the Green function are known explicitly. The most important application, however, lies in the calculation of excited states, whose energies correspond directly to the poles of the Green function in the complex frequency plane. The accuracy of results obtained within this framework is only limited by the choice of the exchange-correlation self-energy, which must still be approximated in actual implementations. In this respect, the GW approximation has proved highly successful for systems governed by the Coulomb interaction. It yields band structures of solids, including the band gaps of semiconductors, as well as atomic and molecular ionization energies in very good quantitative agreement with experimental photoemission data.}}, author = {{Schindlmayr, Arno}}, booktitle = {{Many-Electron Approaches in Physics, Chemistry and Mathematics}}, editor = {{Bach, Volker and Delle Site, Luigi}}, isbn = {{978-3-319-06378-2}}, issn = {{2352-3905}}, pages = {{343--357}}, publisher = {{Springer}}, title = {{{The GW approximation for the electronic self-energy}}}, doi = {{10.1007/978-3-319-06379-9_19}}, volume = {{29}}, year = {{2014}}, } @article{18473, abstract = {{We investigate the band dispersion and related electronic properties of picene single crystals within the GW approximation for the electronic self-energy. The width of the upper highest occupied molecular orbital (HOMOu) band along the Γ–Y direction, corresponding to the b crystal axis in real space along which the molecules are stacked, is determined to be 0.60 eV and thus 0.11 eV larger than the value obtained from density-functional theory. As in our recent study of rubrene using the same methodology [S. Yanagisawa, Y. Morikawa, and A. Schindlmayr, Phys. Rev. B 88, 115438 (2013)], this increase in the bandwidth is due to the strong variation of the GW self-energy correction across the Brillouin zone, which in turn reflects the increasing hybridization of the HOMOu states of neighboring picene molecules from Γ to Y. In contrast, the width of the lower HOMO (HOMOl) band along Γ–Y remains almost unchanged, consistent with the fact that the HOMOl(Γ) and HOMOl(Y) states exhibit the same degree of hybridization, so that the nodal structure of the wave functions and the matrix elements of the self-energy correction are very similar.}}, author = {{Yanagisawa, Susumu and Morikawa, Yoshitada and Schindlmayr, Arno}}, issn = {{1347-4065}}, journal = {{Japanese Journal of Applied Physics}}, number = {{5S1}}, publisher = {{IOP Publishing and The Japan Society of Applied Physics}}, title = {{{Theoretical investigation of the band structure of picene single crystals within the GW approximation}}}, doi = {{10.7567/jjap.53.05fy02}}, volume = {{53}}, year = {{2014}}, } @inbook{18474, author = {{Friedrich, Christoph and Schindlmayr, Arno}}, booktitle = {{Computing Solids: Models, ab initio Methods and Supercomputing}}, editor = {{Blügel, Stefan and Helbig, Nicole and Meden, Volker and Wortmann, Daniel}}, isbn = {{978-3-89336-912-6}}, issn = {{1866-1807}}, location = {{Jülich}}, pages = {{A4.1--A4.21}}, publisher = {{Forschungszentrum Jülich}}, title = {{{Many-body perturbation theory: The GW approximation}}}, volume = {{74}}, year = {{2014}}, } @article{13525, abstract = {{The frequency-dependent dielectric function and the second-order polarizability tensor of ferroelectric LiNbO3 are calculated from first principles. The calculations are based on the electronic structure obtained from density-functional theory. The subsequent application of the GW approximation to account for quasiparticle effects and the solution of the Bethe-Salpeter equation for the stoichiometric material yield a dielectric function that slightly overestimates the absorption onset and the oscillator strength in comparison with experimental measurements. Calculations at the level of the independent-particle approximation indicate that these deficiencies are, at least, partially related to the neglect of intrinsic defects typical for the congruent material. The second-order polarizability calculated within the independent-particle approximation predicts strong nonlinear coefficients for photon energies above 1.5 eV. The comparison with measured data suggests that the inclusion of self-energy effects in the nonlinear optical response leads to a better agreement with experiments. The intrinsic defects of congruent samples reduce the optical nonlinearities, in particular, for the 21 and 31 tensor components, further improving the agreement between experiments and theory.}}, author = {{Riefer, Arthur and Sanna, Simone and Schindlmayr, Arno and Schmidt, Wolf Gero}}, issn = {{1550-235X}}, journal = {{Physical Review B}}, number = {{19}}, publisher = {{American Physical Society}}, title = {{{Optical response of stoichiometric and congruent lithium niobate from first-principles calculations}}}, doi = {{10.1103/PhysRevB.87.195208}}, volume = {{87}}, year = {{2013}}, } @inbook{18475, abstract = {{The frequency-dependent dielectric function and the second-order polarizability tensor of ferroelectric LiNbO3 are calculated from first principles. The calculations are based on the electronic structure obtained from density-functional theory. The subsequent application of the GW approximation to account for quasiparticle effects and the solution of the Bethe–Salpeter equation yield a dielectric function for the stoichiometric material that slightly overestimates the absorption onset and the oscillator strength in comparison with experimental measurements. Calculations at the level of the independent-particle approximation indicate that these deficiencies are at least partially related to the neglect of intrinsic defects typical for the congruent material. The second-order polarizability calculated within the independent-particle approximation predicts strong nonlinear coefficients for photon energies above 1.5 eV. The comparison with measured data suggests that self-energy effects improve the agreement between experiment and theory. The intrinsic defects of congruent samples reduce the optical nonlinearities, in particular for the 21 and 31 tensor components, further improving the agreement with measured data.}}, author = {{Riefer, Arthur and Rohrmüller, Martin and Landmann, Marc and Sanna, Simone and Rauls, Eva and Vollmers, Nora Jenny and Hölscher, Rebecca and Witte, Matthias and Li, Yanlu and Gerstmann, Uwe and Schindlmayr, Arno and Schmidt, Wolf Gero}}, booktitle = {{High Performance Computing in Science and Engineering ‘13}}, editor = {{Nagel, Wolfgang E. and Kröner, Dietmar H. and Resch, Michael M.}}, isbn = {{978-3-319-02164-5}}, pages = {{93--104}}, publisher = {{Springer}}, title = {{{Lithium niobate dielectric function and second-order polarizability tensor from massively parallel ab initio calculations}}}, doi = {{10.1007/978-3-319-02165-2_8}}, year = {{2013}}, } @article{18476, abstract = {{We investigate the band dispersion and relevant electronic properties of rubrene single crystals within the GW approximation. Due to the self-energy correction, the dispersion of the highest occupied molecular orbital (HOMO) band increases by 0.10 eV compared to the dispersion of the Kohn-Sham eigenvalues within the generalized gradient approximation, and the effective hole mass consequently decreases. The resulting value of 0.90 times the electron rest mass along the Γ-Y direction in the Brillouin zone is closer to experimental measurements than that obtained from density-functional theory. The enhanced bandwidth is explained in terms of the intermolecular hybridization of the HOMO(Y) wave function along the stacking direction of the molecules. Overall, our results support the bandlike interpretation of charge-carrier transport in rubrene.}}, author = {{Yanagisawa, Susumu and Morikawa, Yoshitada and Schindlmayr, Arno}}, issn = {{1550-235X}}, journal = {{Physical Review B}}, number = {{11}}, publisher = {{American Physical Society}}, title = {{{HOMO band dispersion of crystalline rubrene: Effects of self-energy corrections within the GW approximation}}}, doi = {{10.1103/PhysRevB.88.115438}}, volume = {{88}}, year = {{2013}}, } @article{18479, abstract = {{The GW approximation for the electronic self-energy is an important tool for the quantitative prediction of excited states in solids, but its mathematical exploration is hampered by the fact that it must, in general, be evaluated numerically even for very simple systems. In this paper I describe a nontrivial model consisting of two electrons on the surface of a sphere, interacting with the normal long-range Coulomb potential, and show that the GW self-energy, in the absence of self-consistency, can in fact be derived completely analytically in this case. The resulting expression is subsequently used to analyze the convergence of the energy gap between the highest occupied and the lowest unoccupied quasiparticle orbital with respect to the total number of states included in the spectral summations. The asymptotic formula for the truncation error obtained in this way, whose dominant contribution is proportional to the cutoff energy to the power −3/2, may be adapted to extrapolate energy gaps in other systems.}}, author = {{Schindlmayr, Arno}}, issn = {{1550-235X}}, journal = {{Physical Review B}}, number = {{7}}, publisher = {{American Physical Society}}, title = {{{Analytic evaluation of the electronic self-energy in the GW approximation for two electrons on a sphere}}}, doi = {{10.1103/PhysRevB.87.075104}}, volume = {{87}}, year = {{2013}}, } @article{18542, abstract = {{We present recent advances in numerical implementations of hybrid functionals and the GW approximation within the full-potential linearized augmented-plane-wave (FLAPW) method. The former is an approximation for the exchange–correlation contribution to the total energy functional in density-functional theory, and the latter is an approximation for the electronic self-energy in the framework of many-body perturbation theory. All implementations employ the mixed product basis, which has evolved into a versatile basis for the products of wave functions, describing the incoming and outgoing states of an electron that is scattered by interacting with another electron. It can thus be used for representing the nonlocal potential in hybrid functionals as well as the screened interaction and related quantities in GW calculations. In particular, the six-dimensional space integrals of the Hamiltonian exchange matrix elements (and exchange self-energy) decompose into sums over vector–matrix–vector products, which can be evaluated easily. The correlation part of the GW self-energy, which contains a time or frequency dependence, is calculated on the imaginary frequency axis with a subsequent analytic continuation to the real axis or, alternatively, by a direct frequency convolution of the Green function G and the dynamically screened Coulomb interaction W along a contour integration path that avoids the poles of the Green function. Hybrid-functional and GW calculations are notoriously computationally expensive. We present a number of tricks that reduce the computational cost considerably, including the use of spatial and time-reversal symmetries, modifications of the mixed product basis with the aim to optimize it for the correlation self-energy and another modification that makes the Coulomb matrix sparse, analytic expansions of the interaction potentials around the point of divergence at k=0, and a nested density and density-matrix convergence scheme for hybrid-functional calculations. We show CPU timings for prototype semiconductors and illustrative results for GdN and ZnO. }}, author = {{Friedrich, Christoph and Betzinger, Markus and Schlipf, Martin and Blügel, Stefan and Schindlmayr, Arno}}, issn = {{1361-648X}}, journal = {{Journal of Physics: Condensed Matter}}, number = {{29}}, publisher = {{IOP Publishing}}, title = {{{Hybrid functionals and GW approximation in the FLAPW method}}}, doi = {{10.1088/0953-8984/24/29/293201}}, volume = {{24}}, year = {{2012}}, } @inproceedings{4048, abstract = {{We present an ab-initio method for calculating nonlinear and nonlocal optical effects in metallic slabs with sub-wavelength thickness. We find a strong localization of the second-harmonic current at the metal-vacuum interface.}}, author = {{Wand, Mathias and Schindlmayr, Arno and Meier, Torsten and Förstner, Jens}}, booktitle = {{CLEO:2011 - Laser Applications to Photonic Applications }}, isbn = {{978-1-4577-1223-4}}, issn = {{2160-8989}}, keywords = {{tet_topic_shg}}, location = {{Baltimore, Maryland, United States}}, publisher = {{Optical Society of America}}, title = {{{Theoretical approach to the ultrafast nonlinear optical response of metal slabs}}}, doi = {{10.1364/CLEO_AT.2011.JTuI59}}, year = {{2011}}, } @article{4091, abstract = {{We present a nonequilibrium ab initio method for calculating nonlinear and nonlocal optical effects in metallic slabs with a thickness of several nanometers. The numerical analysis is based on the full solution of the time‐dependent Kohn–Sham equations for a jellium system and allows to study the optical response of metal electrons subject to arbitrarily shaped intense light pulses. We find a strong localization of the generated second‐harmonic current in the surface regions of the slabs. }}, author = {{Wand, Mathias and Schindlmayr, Arno and Meier, Torsten and Förstner, Jens}}, issn = {{1521-3951}}, journal = {{Physica Status Solidi B}}, keywords = {{tet_topic_shg}}, number = {{4}}, pages = {{887--891}}, publisher = {{Wiley-VCH}}, title = {{{Simulation of the ultrafast nonlinear optical response of metal slabs}}}, doi = {{10.1002/pssb.201001219}}, volume = {{248}}, year = {{2011}}, }