Optimization of effective atom centered potentials for london dispersion forces in density functional theory

We add an effective atom-centered nonlocal term to the exchange-correlation potential in order to cure the lack of London dispersion forces in standard density functional theory. Calibration of this long-range correction is performed using density functional perturbation theory and an arbitrary reference. Without any prior assignment of types and structures of molecular fragments, our corrected generalized gradient approximation density functional theory calculations yield correct equilibrium geometries and dissociation energies of argon-argon, benzene-benzene, graphite-graphite, and argon-benzene complexes.     

Quantum-mechanical calculations based on density functional theory

The basic concepts of density functional theory (DFT) together with the local density approximation (LDA) and the recent improvement in form of the generalized gradient approximations (GGA) are discussed. Band structure calculations using the full-potential linearized augmented plane wave (FP-LAPW) method are presented in relation to pseudopotential schemes both corresponding to T=0. For finite temperatures the most advanced technique is the Car-Parrinello (CP) molecular dynamics (MD) approach, e.g. in its projector augmented wave (PAW) implementation. In CP-MD simulations nuclear motion and the electronic degrees of freedom are treated within one formalism. Such DFT calculations are illustrated for selected examples, including the breathing mode of BaBiO₃. the phase transition in SrTiO₃ and VO₂ and the Li diffusion in the superionic conductor Li₃N studied by conventional and CP molecular dynamics.     

Geometric theory of columnar phases on curved substrates

We study thin self-assembled columns constrained to lie on a curved, rigid substrate. The curvature presents no local obstruction to equally spaced columns in contrast with curved crystals for which the crystalline bonds are frustrated. Instead, the vanishing compressional strain of the columns implies that their normals lie on geodesics which converge (diverge) in regions of positive (negative) Gaussian curvature, in analogy to the focusing of light rays by a lens. We show that the out of plane bending of the cylinders acts as an effective ordering field.     

Critical assessment of mean field models based on effective interactions

Mean field schemes, from simple Hartree—Fock plus random phase approximation calculations of the ground and excited states to more sophisticated approaches which include pairing as well, have been popular for quite a long time. In these models, the input is an effective interaction. We still lack a precise link between this interaction and a more fundamental theory; however, there have been various new recent attempts to correlate empirical pieces of evidence about nuclear (and neutron) matter, or experimental results, with the properties of the effective interactions. In this contribution, we claim that, while we have indeed made some progress in our understanding of certain features of the interactions, we are still missing a clue about its proper density dependence and about its isovector properties. The text was submitted by the author in English.     

Novel surface tension isotherm for surfactants based on local density functional theory

In this Letter we report a new general method for calculating of surface tension isotherms in the presence of surfactants, based on a local density functional. We illustrate this method by deriving the interfacial tension isotherm for nonionic surfactants at an air-water or oil-water interface by using the self-consistent field theory of polymer brushes. We consider a particular case of local density functional to calculate explicitly how the interfacial tension and the surfactant adsorption depend on the surfactant bulk concentration. Experimental data for the surface tension and the surfactant adsorption isotherm for nonionic surfactants were interpreted with the help of the new isotherm. Very good agreement between the adsorption of n-dodecyl pentaoxyethylene glycol ether (C12E5) at an air-water interface, calculated from the surface tension isotherm and small-angle neutron-scattering is obtained.     

First-principles simulations for attosecond photoelectron spectroscopy based on time-dependent density functional theory

We develop a first-principles simulation method for attosecond time-resolved photoelectron spectroscopy. This method enables us to directly simulate the whole experimental processes, including excitation, emission and detection on equal footing. To examine the performance of the method, we use it to compute the reconstruction of attosecond beating by interference of two-photon transitions (RABBITT) experiments of gas-phase Argon. The computed RABBITT photoionization delay is in very good agreement with recent experimental results from [Klünder et al., Phys. Rev. Lett. 106, 143002 (2011)] and [Guénot et al., Phys. Rev. A 85, 053424 (2012)]. This indicates the significance of a fully-consistent theoretical treatment of the whole measurement process to properly describe experimental observables in attosecond photoelectron spectroscopy. The present framework opens the path to unravel the microscopic processes underlying RABBITT spectra in more complex materials and nanostructures.     

Large-scale atomistic simulations of nanostructured materials based on divide-and-conquer density functional theory

A linear-scaling algorithm based on a divide-and-conquer (DC) scheme is designed to perform large-scale molecular-dynamics simulations, in which interatomic forces are computed quantum mechanically in the framework of the density functional theory (DFT). This scheme is applied to the thermite reaction at an Al/Fe₂O₃ interface. It is found that mass diffusion and reaction rate at the interface are enhanced by a concerted metal-oxygen flip mechanism. Preliminary simulations are carried out for an aluminum particle in water based on the conventional DFT, as a target system for large-scale DC-DFT simulations. A pair of Lewis acid and base sites on the aluminum surface preferentially catalyzes hydrogen production in a low activation-barrier mechanism found in the simulations.     

Efficient formalism for large-scale ab initio molecular dynamics based on time-dependent density functional theory

A new "on the fly" method to perform Born-Oppenheimer ab initio molecular dynamics (AIMD) simulations is presented. Inspired by Ehrenfest dynamics in time-dependent density functional theory, the electronic orbitals are evolved by a Schr?dinger-like equation, where the orbital time derivative is multiplied by a parameter. This parameter controls the time scale of the fictitious electronic motion and speeds up the calculations with respect to standard Ehrenfest dynamics. In contrast with other methods, wave function orthogonality needs not be imposed as it is automatically preserved, which is of paramount relevance for large-scale AIMD simulations.     

Development of electron-proton density functionals for multicomponent density functional theory

We present a strategy for the development of electron-proton density functionals in multicomponent density functional theory, treating electrons and selected nuclei quantum mechanically without the Born-Oppenheimer approximation. An electron-proton functional is derived using an explicitly correlated electron-proton pair density. This functional provides accurate hydrogen nuclear densities, thereby enabling reliable calculations of molecular properties. This approach is potentially applicable to relatively large molecular systems with key hydrogen nuclei treated quantum mechanically.     

Instability of a system and its estimation in terms of the hybrid density functional theory method: a magnetic effective density functional (MEDF) approach

Potential curves and high and low spin energy gaps for radical clusters were calculated by spin polarized molecular orbital methods. Through-space effective exchange integrals (J _ab) and relative energies of spin projected low spin states by post-Hartree-Fock (HF) calculation were reproduced by the hybrid density functional theory (DFT) method. The hybrid parameters that could reproduce post-HF values such as UCCSD(T)'s for each model had close relations with the instabilities of those systems. Information entropy and related chemical indices were used to estimate the magnitude of the instabilities. A magnetic effective density functional (MEDF) scheme for spin clusters was proposed for practical computation of J _ab values in molecular magnetic materials.     

Galloping of iced quad-conductors bundles based on curved beam theory

Galloping refers to wind-induced, low-frequency, large-amplitude oscillations that have been more frequently observed for a bundle conductor than for a single conductor. In the present work two different models are built to investigate the galloping of a bundle conductor: (1) a finite curved beam element method and (2) a hybrid model based on curved beam element theory. The finite curved beam element model is effective in dealing with the spacers between the bundled conductors and the joint between the conductors and spacers that can be simulated as a rigid joint or a hinge. Furthermore, the finite curved beam element model can be used to deal with large deformation. The hybrid model invokes the small deformation hypothesis and has a high computational efficiency. A hybrid model based on conventional cable element theory is also programmed to be compared with the aforementioned models based on curved beam element theory. Numerical examples are presented to assess the accuracy of the different models in predicting the equilibrium conductor position, natural frequencies and galloping amplitude. The results show that the curved beam element models, involving more degrees of freedom and coupling of translational and torsional motion, are more accurate at simulating the static and dynamic characters of an iced quad-conductor bundle. The use of hinges, rather than rigid connections, reduces the structural response amplitudes of a galloping conductor bundle.     

A continuum thermal stress theory for crystals based on interatomic potentials

This paper presents a new continuum thermal stress theory for crystals based on interatomic potentials.The effect of finite temperature is taken into account via a harmonic model.An EAM potential for copper is adopted in this paper and verified by computing the effect of the temperature on the specific heat,coefficient of thermal expansion and lattice constant.Then we calculate the elastic constants of copper at finite temperature.The calculation results are in good agreement with experimental data.The thermal stress theory is applied to an anisotropic crystal graphite,in which the Brenner potential is employed.Temperature dependence of the thermodynamic properties,lattice constants and thermal strains for graphite is calculated.The calculation results are also in good agreement with experimental data.     

Thirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionals

In the past 30 years, Kohn–Sham density functional theory has emerged as the most popular electronic structure method in computational chemistry. To assess the ever-increasing number of approximate exchange-correlation functionals, this review benchmarks a total of 200 density functionals on a molecular database (MGCDB84) of nearly 5000 data points. The database employed, provided as Supplemental Data, is comprised of 84 data-sets and contains non-covalent interactions, isomerisation energies, thermochemistry, and barrier heights. In addition, the evolution of non-empirical and semi-empirical density functional design is reviewed, and guidelines are provided for the proper and effective use of density functionals. The most promising functional considered is ωB97M-V, a range-separated hybrid meta-GGA with VV10 nonlocal correlation, designed using a combinatorial approach. From the local GGAs, B97-D3, revPBE-D3, and BLYP-D3 are recommended, while from the local meta-GGAs, B97M-rV is the leading choice, followed by MS1-D3 and M06-L-D3. The best hybrid GGAs are ωB97X-V, ωB97X-D3, and ωB97X-D, while useful hybrid meta-GGAs (besides ωB97M-V) include ωM05-D, M06-2X-D3, and MN15. Ultimately, today's state-of-the-art functionals are close to achieving the level of accuracy desired for a broad range of chemical applications, and the principal remaining limitations are associated with systems that exhibit significant self-interaction/delocalisation errors and/or strong correlation effects.     

Electron dynamics triggered by double attosecond pulses: Simulations based on time-dependent density functional theory

In order to observe the high-field effect, the external laser field must reach its peak intensity before the electron ionization. To this end, it is important to reduce pulse duration to typical attosecond timescale. In this paper, the interaction electron dynamics between attosecond pulses and dielectric is investigated within the time-dependent density functional theory. Taking the CaF₂ crystal as an example, we give a comparison of electron dynamics response between single and double pulses. Moreover, the nonlinear energy absorption and electron excitation processes are simulated by adjusting the polarization direction of the sub-pulse. Present results demonstrate that the double pulses show lower electron excitation and energy absorption than the single pulse, which is in accordance with experimental higher ablation threshold and smaller heat-affected zones of the double pulses. In addition, the curves of final excited electron number and energy absorption exhibit the quasi-symmetry about the axis of 180°, which has not been reported yet.     

A density functional theory study of formaldehyde adsorption on ceria

Molecular adsorption of formaldehyde on the stoichiometric CeO₂(1 1 1) and CeO₂(1 1 0) surfaces was studied using periodic density functional theory. Two adsorption modes (strong chemisorbed and weak physisorbed) were identified on both surfaces. This is consistent with recent experimental observations. On the (1 1 1) surface, formaldehyde strongly chemisorbs with an adsorption energy of 0.86 eV to form a dioxymethylene-like structure, in which a surface O lifts from the surface to bind with the C of formaldehyde. A weak physisorbed state with adsorption energy of 0.28 eV was found with the O of formaldehyde interacting with a surface Ce. On the (1 1 0) surface, dioxymethyelene formation was also observed, with an adsorption energy of 1.31 eV. The weakly adsorbed state of formaldehyde on the (1 1 0) surface was energetically comparable to the weak adsorption state on the (1 1 1) surface. Analysis of the local density of states and charge density differences after adsorption shows that strong covalent bonding occurs between the C of formaldehyde and surface O when dioxymethylene forms. Calculated vibrational frequencies also confirm dioxymethylene formation. Our results show that as the coverage increases, the adsorption of formaldehyde on the (1 1 1) surface becomes weak, but is nearly unaffected on the (1 1 0) surface.