Adiabatic density-functional perturbation theory

Electronic and vibrational properties of the highTcsuperconductor Bi2Sr2CaCu2O8: anab initiostudy

In this work,ab initiocalculations were performed in order to study the vibrational spectra of the Bi₂Sr₂CaCu₂O₈(Bi2212) compound. A structural modulation correction on some atomic positions, producing a distorted structure with lower symmetry, is used for the calculation. We argue that this correction allows to account for an average effect of the incommensurate superstructure, generating a more accurate representation of the real unit cell observed in this compound. A complete and conclusive vibrational assignment is performed, discussing the correspondences with previous experimental and theoretical reports. A brief analysis of the electronic density of states and band structure comparing the tetragonal and distorted unit cell is also included.

Analytic calculation and analysis of atomic polar tensors for molecules and materials using the Gaussian and plane waves approach

The evaluation of atomic polar tensors and Born Effective Charge (BEC) tensors from Density Functional Perturbation Theory (DFPT) has been implemented in the CP2K code package. This implementation is based on a combination of the Gaussian and plane wave approach for the description of basis functions and arising potentials. The presence of non-local pseudo-potentials has been considered, as well as contributions arising from the basis functions being centered on the atoms. Simulations of both periodic and non-periodic systems have been implemented and carried out. Dipole strengths and infrared absorption spectra have been calculated for two isomers of the tripeptide Ser-Pro-Ala using DFPT and are compared to the results of standard vibrational analyses using finite differences. The spectra are then decomposed into five subsets by employing localized molecular orbitals/maximally localized Wannier functions, and the results are discussed. Moreover, group coupling matrices are employed for visualization of results. Furthermore, the BECs and partial charges of the surface atoms of a periodic (101) anatase (TiO₂) slab have been investigated in a periodic framework.

Using High Multipolar Orders to Reconstruct the Sound Velocity in Piezoelectrics from Lattice Dynamics

Information over the phonon band structure is crucial to predicting many thermodynamic properties of materials, such as thermal transport coefficients. Highly accurate phonon dispersion curves can be, in principle, calculated in the framework of density-functional perturbation theory. However, well-established techniques can run into trouble (or even catastrophically fail) in the case of piezoelectric materials, where the acoustic branches hardly reproduce the physically correct sound velocity. Here we identify the culprit in the higher-order multipolar interactions between atoms and demonstrate an effective procedure that fixes the aforementioned issue. Our strategy drastically improves the predictive power of perturbative lattice-dynamical calculations in piezoelectric crystals and is directly implementable for high-throughput generation of materials databases.

Notes on density matrix perturbation theory

Density matrix perturbation theory (DMPT) is known as a promising alternative to the Rayleigh-Schrödinger perturbation theory, in which the sum-over-states (SOS) is replaced by algorithms with perturbed density matrices as the input variables. In this article, we formulate and discuss three types of DMPT, with two of them based only on density matrices: the approach of Kussmann and Ochsenfeld [J. Chem. Phys. 127, 054103 (2007)] is reformulated via the Sylvester equation and the recursive DMPT of Niklasson and Challacombe [Phys. Rev. Lett. 92, 193001 (2004)] is extended to the hole-particle canonical purification (HPCP) from Truflandier et al. [J. Chem. Phys. 144, 091102 (2016)]. A comparison of the computational performances shows that the aforementioned methods outperform the standard SOS. The HPCP-DMPT demonstrates stable convergence profiles but at a higher computational cost when compared to the original recursive polynomial method.

Magnon-phonon interaction in antiferromagnetic two-dimensional MXenes

Antiferromagnetic material possesses excellent robustness to an external magnetic field perturbation, which makes it promising in application of spintronic devices. The magnon-phonon interaction plays a vital role in spintronic devices. In this work, we performed first-principles calculation to study the effect of magnon-phonon interaction on magnon spectra of the antiferromagnetic MXenes Cr₂TiC₂FCl, and calculated the phonon dominated magnon relaxation time based on the magnon spectra broadening. Due to the large exchange constants across Cr-Cr pairs, high magnon energy is found in Cr₂TiC₂FCl. We find that compared with the acoustic magnons, the optical magnons have stronger interaction with phonon modes. Moreover, relaxation time of optical magnons and acoustic magnons have quite different wavevector dependence. Our results about spin coupling to specific phonon polarizations can shed light on the understanding of magnon damping and energy dissipation in two-dimensional antiferromagnetic materials.

Analytical approach to phonon calculations in the SCC-DFTB framework

Detailed derivation of the analytical, reciprocal-space approach of Hessian calculation within the self-consistent-charge density-functional based tight-binding framework (SCC-DFTB) is presented. This approach provides an accurate and efficient way for obtaining the SCC-DFTB Hessian of periodic systems. Its superiority with respect to the traditional numerical force differentiation method is demonstrated for doped graphene, graphene nanoribbons, boron-nitride nanotubes, bulk zinc-oxide, and other systems.

Effects of molecular adsorption on the spin-wave spectrum and magnon relaxation in two-dimensional Cr2Ge2Te6

In this work, we performed detailed first-principles calculation and theoretical analysis to investigate the effect of molecular adsorption on the spin-wave spectrum and magnon relaxation in a Cr2Ge2Te6 (CGT) monolayer. It is found that NH3, NO, and NO2 adsorption can enhance the exchange constant of CGT, which can result in a blue-shift in the spin-wave spectrum. At 30 K, by means of a thorough investigation of many possible lattice configurations excited by thermal fluctuation, we identify the magnon scattering rate from the intrinsic lattice vibrational modes, and find that the relaxation of optical and acoustic magnons exhibits a completely different wave vector dependence. Moreover, although the adsorption of NO2 and NH3 molecules has a negligible influence on the magnon-phonon interaction, the adsorption of NO molecules results in a significant increase in magnon scattering strength. In the long-wavelength limit, the interlayer vibrational modes induced by NO adsorption increase the magnon-phonon scattering strength by ∼12.7%. The remarkable interlayer magnon-phonon interaction is ascribed to the strong CGT-NO coupling and large molecular vibration amplitude. Considering the importance of magnon relaxation time in the application of spin devices, we suggest that both the impacts on the exchange interaction and scattering rate must be considered when manipulating two-dimensional magnets by surface functionalization.

CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations

CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.

Quantum ESPRESSO toward the exascale

Quantum ESPRESSO is an open-source distribution of computer codes for quantum-mechanical materials modeling, based on density-functional theory, pseudopotentials, and plane waves, and renowned for its performance on a wide range of hardware architectures, from laptops to massively parallel computers, as well as for the breadth of its applications. In this paper, we present a motivation and brief review of the ongoing effort to port Quantum ESPRESSO onto heterogeneous architectures based on hardware accelerators, which will overcome the energy constraints that are currently hindering the way toward exascale computing.

In silico Raman spectroscopy of YAlO3 single-crystalline film

The Raman response of the YAlO₃ (YAP) perovskite is modeled by means of periodic density functional theory. A number of different approximations to the exchange-correlation functional are benchmarked against the structural and spectroscopic data as imposing all-electron Gaussian-type basis sets. The WC1LYP functional was found to be superior, particularly outperforming other tested approaches in the prediction of the local structure of the AlO subunits, which reflects in the observed lattice-dynamics. The Raman response is further decomposed into the directional spectra, which are due to different components of the polarizability tensor, and confronted with the experimental Raman spectra, recorded in different scattering geometries of the single-crystalline film of YAP. The in silico lattice dynamics provides the unequivocal assignment of the observed bands with an excellent match to the experimental spectra, allowing for a complete analysis of the underlying phonon modes in terms of their energy, symmetry and the directional activity. The presented analysis serves as a high-quality reference, potentially useful in the future studies of other YAP materials, where Raman spectroscopy along with the X-Ray diffraction is the first method of choice.

Direct evaluation of the isotope effect within the framework of density functional theory for superconductors

Within recent developments of density functional theory, its numerical implementation and of the superconducting density functional theory is nowadays possible to predict the superconducting critical temperature, [Formula: see text], with sufficient accuracy to anticipate the experimental verification. In this paper we present an analytical derivation of the isotope coefficient within the superconducting density functional theory. We calculate the partial derivative of [Formula: see text] with respect to atomic masses. We verified the final expression by means of numerical calculations of isotope coefficient in monatomic superconductors (Pb) as well as polyatomic superconductors (CaC₆). The results confirm the validity of the analytical derivation with respect to the finite difference methods, with considerable improvement in terms of computational time and calculation accuracy. Once the critical temperature is calculated (at the reference mass(es)), various isotope exponents can be simply obtained in the same run. In addition, we provide the expression of interesting quantities like partial derivatives of the deformation potential, phonon frequencies and eigenvectors with respect to atomic masses, which can be useful for other derivations and applications.

Toward First-Principles-Level Polarization Energies in Force Fields: A Gaussian Basis for the Atom-Condensed Kohn-Sham Method

The last 20 years of force field development have shown that even well parametrized classical models need to at least approximate the dielectric response of molecular systems-based, e.g., on atomic polarizabilities-in order to correctly render their structural and dynamic properties. Yet, despite great advances most approaches tend to be based on ad hoc assumptions and often insufficiently capture the dielectric response of the system to external perturbations, such as, e.g., charge carriers in semiconducting materials. A possible remedy was recently introduced with the atom-condensed Kohn-Sham density-functional theory approximated to second order (ACKS2), which is fully derived from first principles. Unfortunately, specifically its reliance on first-principles derived parameters so far precluded the widespread adoption of ACKS2. Opening up ACKS2 for general use, we here present a reformulation of the method in terms of Gaussian basis functions, which allows us to determine many of the ACKS2 parameters analytically. Two sets of parameters depending on exchange-correlation interactions are still calculated numerically, but we show that they could be straightforwardly parametrized owing to the smoothness of the new basis. Our approach exhibits three crucial benefits for future applications in force fields: i) efficiency, ii) accuracy, and iii) transferability. We numerically validate our Gaussian augmented ACKS2 model for a set of small hydrocarbons which shows a very good agreement with density-functional theory reference calculations. To further demonstrate the method's accuracy and transferability for realistic systems, we calculate polarization responses and energies of anthracene and tetracene, two major building blocks in organic semiconductors.

The GW Compendium: A Practical Guide to Theoretical Photoemission Spectroscopy

The GW approximation in electronic structure theory has become a widespread tool for predicting electronic excitations in chemical compounds and materials. In the realm of theoretical spectroscopy, the GW method provides access to charged excitations as measured in direct or inverse photoemission spectroscopy. The number of GW calculations in the past two decades has exploded with increased computing power and modern codes. The success of GW can be attributed to many factors: favorable scaling with respect to system size, a formal interpretation for charged excitation energies, the importance of dynamical screening in real systems, and its practical combination with other theories. In this review, we provide an overview of these formal and practical considerations. We expand, in detail, on the choices presented to the scientist performing GW calculations for the first time. We also give an introduction to the many-body theory behind GW, a review of modern applications like molecules and surfaces, and a perspective on methods which go beyond conventional GW calculations. This review addresses chemists, physicists and material scientists with an interest in theoretical spectroscopy. It is intended for newcomers to GW calculations but can also serve as an alternative perspective for experts and an up-to-date source of computational techniques.

Controlling the Infrared Dielectric Function through Atomic-Scale Heterostructures

Surface phonon polaritons (SPhPs), the surface-bound electromagnetic modes of a polar material resulting from the coupling of light with optic phonons, offer immense technological opportunities for nanophotonics in the infrared (IR) spectral region. However, once a particular material is chosen, the SPhP characteristics are fixed by the spectral positions of the optic phonon frequencies. Here, we provide a demonstration of how the frequency of these optic phonons can be altered by employing atomic-scale superlattices (SLs) of polar semiconductors using AlN/GaN SLs as an example. Using second harmonic generation (SHG) spectroscopy, we show that the optic phonon frequencies of the SLs exhibit a strong dependence on the layer thicknesses of the constituent materials. Furthermore, new vibrational modes emerge that are confined to the layers, while others are centered at the AlN/GaN interfaces. As the IR dielectric function is governed by the optic phonon behavior in polar materials, controlling the optic phonons provides a means to induce and potentially design a dielectric function distinct from the constituent materials and from the effective-medium approximation of the SL. We show that atomic-scale AlN/GaN SLs instead have multiple Reststrahlen bands featuring spectral regions that exhibit either normal or extreme hyperbolic dispersion with both positive and negative permittivities dispersing rapidly with frequency. Apart from the ability to engineer the SPhP properties, SL structures may also lead to multifunctional devices that combine the mechanical, electrical, thermal, or optoelectronic functionality of the constituent layers. We propose that this effort is another step toward realizing user-defined, actively tunable IR optics and sources.

Galvanic replacement of liquid metal galinstan with Pt for the synthesis of electrocatalytically active nanomaterials

The galvanic replacement reaction is a verstile method for the fabrication of bimetallic nanomaterials which is usually limited to solid precursors. Here we report on the galvanic replacement of liquid metal galinstan with Pt which predominantly results in the formation of a Pt5Ga1 material. During the galvanic replacement process an interesting phenomenon was observed whereby a plume of nanomaterial is ejected upwards from the centre of the liquid metal droplet into solution which is due to surface tension gradients on the liquid metal surface that induces surface convection. It was also found that hydrogen gas was liberated during the process facilitated by the formation of the Pt rich nanomaterial which is a highly effective catalyst for the hydrogen evolution reaction (HER). The material was characterised by scanning electron microscopy, transmission electron microscopy, X-ray photoelectron spectroscopy, X-ray diffraction and dynamic light scattering measurements. It was found that Pt5Ga1 was highly effective for the electrochemical oxidation of methanol and ethanol and outperformed a commercial Pt/C catalyst. Density functional theory calculations confirmed that the increased activity is due to the anti poisoning properties of the surface towards CO upon the incorporation of Ga atoms into a Pt catalyst. The use of liquid metals and galvanic replacement offers a simple approach to fabricating Ga based alloy nanomaterials that may have use in many other types of applications.

Piezoelectric Effects in Surface-Engineered Two-Dimensional Group III Nitrides

Piezoelectric effects of two-dimensional (2D) group III-V compounds have received considered attention in recent years because of their wide applications in semiconductor devices. However, they face a problem that only metastable or unstable structures are noncentrosymmetric with piezoelectricity, thus leading to the difficulty in experimental observation. Motivated by the recent advances in the synthesis of 2D group III nitrides, in this paper, for the first time, we study the piezoelectric properties of the 2D group III nitrides (XN, X = Al, Ga, and In) with buckled hexagonal configurations by surface passivation, which are thermodynamically stable. Unlike the previously reported planar graphitic structure, we demonstrate that the hydrogenated 2D nitrides (H-XN-H, X = Al, Ga, and In) exhibit both the in-plane and out-of-plane piezoelectric effects in their monolayer and multilayer structures under an external strain in the basal plane. We further elucidate the underlying mechanism of the piezoelectricity by analyzing the correlations between the piezoelectric coefficients and their structural, electronic, and chemical properties. In addition, we show that H-F cofunctionalization not only enhances the stability, but also significantly improves the ionic polarization because of the charge redistribution, thus leading to large in-plane piezoelectric coefficients in F-XN-H. Our study advances the research in 2D piezoelectric materials and would stimulate more theoretical and experimental efforts in developing effective piezoelectric materials for device applications.

Theoretical prediction of some layered Pa2O5 phases: structure and properties

Density functional theory (DFT) was used to predict and study protactinium pentoxide (Pa₂O₅), which presents a fluorite and layered protactinium oxide-type structure. Although the layered structure has been observed with the isostructural transition Nb and Ta metal pentoxides experimentally, the detailed structure and properties of the layered Pa₂O₅ are not clear and understandable. Our theoretical prediction explored some possible stable structures of the Pa₂O₅ stoichiometry according to the existing M₂O₅ structures (where M is an actinide Np or transition Nb, Ta, and V metal) and replacing the M ions with protactinium ions. The structural, mechanical, thermodynamic and electronic properties including lattice parameters, bulk moduli, elastic constants, entropy and band gaps were predicted for all the simulated structures. Pa₂O₅ in the β-V₂O₅ structure was found to be a competitive structure in terms of stability, whereas Pa₂O₅ in the ζ-Nb₂O₅ structure was found to be the most stable overall. This is consistent with Sellers's experimental observations. In particular, Pa₂O₅ in the ζ-Nb₂O₅ structure is predicted to be charge-transfer insulators. Furthermore, we predict that ζ-Nb₂O₅-structured Pa₂O₅ is the most thermodynamically stable under ambient conditions and pressure.

Density functional theory (DFT) was used to predict and study protactinium pentoxide (Pa₂O₅), which presents a fluorite and layered protactinium oxide-type structure.

Vibrational spectra of Pb2Bi2Te3, PbBi2Te4, and PbBi4Te7 topological insulators: temperature-dependent Raman and theoretical insights from DFT simulations

Insertion of lead and lead telluride in Bi₂Te₃ leads to a change in the thermal conductivity, frequency shift, and the broadening of phonon modes.

Metal and ligand effects on the stability and electronic properties of crystalline two-dimensional metal-benzenehexathiolate coordination compounds

The cohesive energy, phonon spectrum and quantum molecular-dynamic simulation have been used successively to determine whether the crystalline two-dimensional (2D) metal-benzenehexathiolate (M-BHT) coordination compounds are stable or not. The electronic structures of stable M-BHTs and the corresponding inorganic semiconducting materials have been compared. From the point of view of satisfying stoichiometric ratios and saturation of chemical bonds, we designed possible planar molecular structures and demonstrated that there may be two different 2D M-BHTs, i.e. group II-[Formula: see text] and group IV-[Formula: see text]. However, the cohesive energy calculation indicates that the group IV-[Formula: see text] coordination compound cannot be obtained by thermodynamic equilibrium growth. In contrast, [Formula: see text] and [Formula: see text] from the group II-[Formula: see text] have not only thermodynamic stability, but also dynamic stability due to their phonon spectrum with no imaginary frequency. Moreover, they are still the two most stable ones when the bridge atom S of ligand BHT is replaced by the other chalcogens of O, Se and Te. Further studies indicated that [Formula: see text] and [Formula: see text] both have room temperature dynamic stability and exhibit semiconducting. The exceptional stability and relatively narrow band gap make them advantageous over their inorganic counterparts. Our findings open opportunities to search for new 2D planar conducting coordination compound for organic electronic applications.

Isotopic Heft on the B1 l Silent Mode in Ultra-Narrow Gallium Nitride Nanowires

Wurtzite semiconductor compounds have two silent modes, B_{1 l} and B_{1 h}. A silent mode is a vibrational mode that carries neither a dipole moment nor Raman polarizability. Thus, they are forbidden in both infrared reflectivity and Raman spectroscopy. Astonishingly, we detected the B_{1 l} mode in high-quality, ultra-narrow GaN nanowires using resonant Raman scattering, although the B_{1 h} was not observed, and there is no immediate explanation for this asymmetric finding. The Raman experiments were performed using several laser lines from 647 to 325 nm; the latter is a wavelength in which Raman becomes resonant. Actually, we observed the B_{1 l} mode only in resonance, indicating that the appearance of this mode is related to Fröhlich electron-phonon interactions; i.e., a dipole moment emerging in the B_{1 l} silent mode may not be present in the B_{1 h} mode. To shed light onto the physical origin of these observations, we performed density functional theory calculations of the lattice dynamics in GaN. We performed a careful analysis of the different physical mechanisms that allow the forbidden mode to appear to explain the physics underlying the nonzero dipole moment in the B_{1 l} mode, and the reason why this dipole moment is not present in the B_{1 h} mode.

First-principles analysis of a molecular piezoelectric meta-nitroaniline

The piezoelectric and elastic properties of a molecular piezoelectric meta-nitroaniline (mNA) in its single-crystal form were investigated in the framework of first-principles density functional perturbation theory (DFPT). Results support the recent experimental findings those despite being soft and flexible, mNA's piezoelectric coefficients are an order of magnitude greater than that of ZnO and LiNbO₃. A molecular-level insight into the piezoelectric properties of mNA is provided. These results are helpful not only for better understanding mNA, but also for developing new piezoelectric materials.

Computational methods for 2D materials: discovery, property characterization, and application design

The discovery of two-dimensional (2D) materials comes at a time when computational methods are mature and can predict novel 2D materials, characterize their properties, and guide the design of 2D materials for applications. This article reviews the recent progress in computational approaches for 2D materials research. We discuss the computational techniques and provide an overview of the ongoing research in the field. We begin with an overview of known 2D materials, common computational methods, and available cyber infrastructures. We then move onto the discovery of novel 2D materials, discussing the stability criteria for 2D materials, computational methods for structure prediction, and interactions of monolayers with electrochemical and gaseous environments. Next, we describe the computational characterization of the 2D materials' electronic, optical, magnetic, and superconducting properties and the response of the properties under applied mechanical strain and electrical fields. From there, we move on to discuss the structure and properties of defects in 2D materials, and describe methods for 2D materials device simulations. We conclude by providing an outlook on the needs and challenges for future developments in the field of computational research for 2D materials.

Advanced capabilities for materials modelling with Quantum ESPRESSO

Quantum EXPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the-art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudopotential and projector-augmented-wave approaches. Quantum EXPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement their ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.

Temperature-Induced Large Broadening and Blue Shift in the Electronic Band Structure and Optical Absorption of Methylammonium Lead Iodide Perovskite

The power conversion efficiency of hybrid halide perovskite solar cells is profoundly influenced by the operating temperature. Here we investigate the temperature influence on the electronic band structure and optical absorption of cubic CH₃NH₃PbI₃ from first-principles by accounting for both the electron-phonon interaction and thermal expansion. Within the framework of density functional perturbation theory, the electron-phonon coupling induces slightly enlarged band gap and strongly broadened electronic relaxation time as temperature increases. The large broadening effect is mainly due to the presence of cation organic atoms. Consequently, the temperature-dependent absorption peak exhibits blue-shift position, decreased amplitude, and broadened width. This work uncovers the atomistic origin of temperature influence on the optical absorption of cubic CH₃NH₃PbI₃ and can provide guidance to design high-performance hybrid halide perovskite solar cells at different operating temperatures.

Raman and DRIFT spectra, vibrational assignments and quantum mechanical calculations of centrosymmetric meso-2,3-Dimercaptosuccinic acid

The Raman spectrum (3700-100cm^-1) of meso-2,3-Dimercaptosuccinic acid (meso-DMSA; C₄H₆O₄S₂) was recorded in the solid phase using 514.5 and 785nm excitation lines. Whereas, the DRIFT spectrum (4000-400cm^-1) of the sample powdered in KBr was obtained. Moreover, DFT-B3LYP/6-31G(d) geometry optimization and frequency calculations were carried out for centrosymmetric trans (C_i), gauche (C₁; G⁺/G^-) and eclipsed (C_s; E_f and C₁; E⁺/E^-) rotational isomers in favor of a trans conformation, the least energy with real frequencies. However, other conformers were found at either local minima or local maxima as a result of the rotation of carboxyl, hydroxyl and thiol groups according to a potential energy surface scan. Moreover, an imaginary wavenumber was predicted; therefore, they are considered transition states. On the other hand, the mass spectrum of the sample dissolved in an acetonitrile/methanol mixture reveal 4-6% dimer through intermolecular hydrogen bonding interactions via the dicarboxylic groups. Therefore, we have modeled the complex structure obeying C_i restricted symmetry for an isolated dimer unit using DFT-B3LYP/6-31G(d) and for two molecules per unit cell in the solid phase implementing DFT-PBE functional. Thus, the meso-DMSA forms long strands in which individual molecules are bonded together at each termini through hydrogen bonding. Aided by normal coordinate analysis, complete vibrational assignments were provided herein which support C_i configuration of meso-DMSA in the solid state which found consistent with the observed broadening, composite, split bands, and the mutual exclusion rule.

Ab Initio Approach to Second-order Resonant Raman Scattering Including Exciton-Phonon Interaction

Raman spectra obtained by the inelastic scattering of light by crystalline solids contain contributions from first-order vibrational processes (e.g. the emission or absorption of one phonon, a quantum of vibration) as well as higher-order processes with at least two phonons being involved. At second order, coupling with the entire phonon spectrum induces a response that may strongly depend on the excitation energy, and reflects complex processes more difficult to interpret. In particular, excitons (i.e. bound electron-hole pairs) may enhance the absorption and emission of light, and couple strongly with phonons in resonance conditions. We design and implement a first-principles methodology to compute second-order Raman scattering, incorporating dielectric responses and phonon eigenstates obtained from density-functional theory and many-body theory. We demonstrate our approach for the case of silicon, relating frequency-dependent relative Raman intensities, that are in excellent agreement with experiment, to different vibrations and regions of the Brillouin zone. We show that exciton-phonon coupling, computed from first principles, indeed strongly affects the spectrum in resonance conditions. The ability to analyze second-order Raman spectra thus provides direct insight into this interaction.

High-quality PVD graphene growth by fullerene decomposition on Cu foils

We present a new protocol to grow large-area, high-quality single-layer graphene on Cu foils at relatively low temperatures. We use C₆₀ molecules evaporated in ultra high vacuum conditions as carbon source. This clean environment results in a strong reduction of oxygen-containing groups as depicted by X-ray photoelectron spectroscopy (XPS). Unzipping of C₆₀ is thermally promoted by annealing the substrate at 800ºC during evaporation. The graphene layer extends over areas larger than the Cu crystallite size, although it is changing its orientation with respect to the surface in the wrinkles and grain boundaries, producing a modulated ring in the low energy electron diffraction (LEED) pattern. This protocol is a self-limiting process leading exclusively to one single graphene layer. Raman spectroscopy confirms the high quality of the grown graphene. This layer exhibits an unperturbed Dirac-cone with a clear n-doping of 0.77 eV, which is caused by the interaction between graphene and substrate. Density functional theory (DFT) calculations show that this interaction can be induced by a coupling between graphene and substrate at specific points of the structure leading to a local sp³ configuration, which also contribute to the D-band in the Raman spectra.

Influence of point defects and impurities on the dynamical stability of δ-plutonium

We use first-principles calculations to provide direct evidence of the effect of aluminum, gallium, iron and uranium on the dynamical stability of δ-plutonium. We first show that the δ phase is dynamically unstable at low temperature, as seen in experiments, and that this stability directly depends on the plutonium 5f orbital occupancies. Then, we demonstrate that both aluminum and gallium stabilize the δ phase, contrary to iron. As for uranium, which is created during self-irradiation and whose effect on plutonium has yet to be understood, we show that it leaves a few unstable vibrational modes and that higher concentrations lead to an almost complete stabilization. Finally, we provide an attempt at a consistent analysis of the experimental Pu-Ga phonon density of states. We show that the presence of gallium can reproduce only partially the experimental measurements, and we investigate how point defects, such as interstitials and vacancies, affect the calculated phonon density of states.

Uncovering the Thermo-Kinetic Origins of Phase Ordering in Mixed-Valence Antimony Tetroxide by First-Principles Modeling

Phase ordering in the mixed-valence oxide Sb₂O₄ has been examined by density functional theory (DFT) calculations. We find that the ground-state total energies of the two phases (α and β) are almost degenerate and are highly sensitive to the choice of the approximation to the exchange correlation (xc) functional used in our calculations. Interestingly, with the inclusion of the zero-point energy corrections, the α phase is predicted to be the ground state polymorph for most xc functionals used. We also illustrate the pronounced stereochemical activity of Sb in these polymorphs of Sb₂O₄, setting an exception to the Keve and Skapski rule. Here, we find that the actual bonding in the α phase is more asymmetric, while the anomalous stability of the β phase could be rationalized from kinetic considerations. We find a non-negligible activation barrier for this α-β phase transition, and the presence of a saddle point (β phase) supports the separation of Sb(III) over a continuous phase transition, as observed in experiments.

Nuclear dynamics and phase polymorphism in solid formic acid

We apply a unique sequence of structural and dynamical neutron-scattering techniques, augmented with density-functional electronic-structure calculations, to establish the degree of polymorphism in an archetypal hydrogen-bonded system - crystalline formic acid. Using this combination of experimental and theoretical techniques, the hypothesis by Zelsmann on the coexistence of the β₁ and β₂ phases above 220 K is tested. Contrary to the postulated scenario of proton-transfer-driven phase coexistence, the emerging picture is one of a quantitatively different structural change over this temperature range, whereby the loosening of crystal packing promotes temperature-induced shearing of the hydrogen-bonded chains. The presented work, therefore, solves a fifty-year-old puzzle and provides a suitable framework for the use neutron-Compton-scattering techniques in the exploration of phase polymorphism in condensed matter.

Temperature dependence of X-ray absorption and nuclear magnetic resonance spectra: probing quantum vibrations of light elements in oxides

Probing the quantum thermal fluctuations of nuclei in light-element oxides using XANES and NMR spectroscopies.

Modeling non-harmonic behavior of materials from experimental inelastic neutron scattering and thermal expansion measurements

Based on thermodynamic principles, we derive expressions quantifying the non-harmonic vibrational behavior of materials, which are rigorous yet easily evaluated from experimentally available data for the thermal expansion coefficient and the phonon density of states. These experimentally-derived quantities are valuable to benchmark first-principles theoretical predictions of harmonic and non-harmonic thermal behaviors using perturbation theory, ab initio molecular-dynamics, or Monte-Carlo simulations. We illustrate this analysis by computing the harmonic, dilational, and anharmonic contributions to the entropy, internal energy, and free energy of elemental aluminum and the ordered compound [Formula: see text] over a wide range of temperature. Results agree well with previous data in the literature and provide an efficient approach to estimate anharmonic effects in materials.

Unveiling the Room-Temperature Magnetoelectricity of Troilite FeS

We report on a first-principles study of the troilite phase of iron sulfide (FeS). We show that even if, a few decades ago, this material was thought to be ferroelectric, the structural transition from the high P6_{3}/mmc to the low P6[over ¯]2c symmetry phase does not involve polar instabilities, though the space inversion center symmetry is broken. Our calculations and symmetry analysis nevertheless reveal that FeS is magnetoelectric at room temperature with a response larger than the prototypical room-temperature magnetoelectric crystal Cr_{2}O_{3}. We also show that the spin channel decomposition of the polarization exhibits nonzero values in the opposite direction in FeS, which is actually a general hint of the presence of a magnetoelectric monopole in diagonal magnetoelectrics.