Search

With the aim of progressing toward a practical implementation of an effective quantum-electrodynamics (QED) theory of atoms and molecules, which includes the effects of vacuum polarization through the creation of virtual electron-positron pairs but without the explicit photon degrees of freedom, we study a one-dimensional effective QED model of the hydrogen-like atom with delta-potential interactions. This model resembles the three-dimensional effective QED theory with Coulomb interactions while being substantially simpler. We provide some mathematical details about the definition of this model, calculate the vacuum-polarization density, and the Lamb-type shift of the bound-state energy, correcting and extending results of previous works. We also study the approximation of the model in a finite plane-wave basis, and in particular we discuss the basis convergence of the bound-state energy and eigenfunction, of the vacuum-polarization density, and of the Lamb-type shift of the bound-state energy. We highlight the difficulty of converging the vacuum-polarization density in a finite basis and we propose a way to improve it. The present work could give hints on how to perform similar calculations for the three-dimensional effective QED theory of atoms and molecules.

Strain engineering is a powerful tool routinely used to control and enhance properties such as ferroelectricity, magnetic ordering, or metal-insulator transitions. Epitaxial strain in thin films allows manipulation of in-plane lattice parameters, achieving strain values generally up to 4\%, and above in some specific cases. In polycrystalline films, which are more suitable for functional applications due to their lower fabrication costs, strains above 1% often cause cracking. This poses challenges for functional property tuning by strain engineering. Helium implantation has been shown to induce negative pressure through interstitial implantation, which increases the unit cell volume and allows for continuous strain tuning with the implanted dose in epitaxial monocrystalline films. However, there have been no studies on the transferability of helium implantation as a strain-engineering technique to polycrystalline films. Here, we demonstrate the technique's applicability for strain engineering beyond epitaxial monocrystalline samples. Helium implantation can trigger an unprecedented lattice parameter expansion of up to 3.2% in polycrystalline BiFeO3 films without causing structural cracks. The film maintains stable ferroelectric properties with doses up to 1E15 He/cm2. This finding underscores the potential of helium implantation in strain engineering polycrystalline materials, enabling cost-effective and versatile applications., Comment: 6 pages, 4 Figures

Traditionally, inapproximability results for $\mathsf{Max\,k\,CSP\!-\!q}$ have been established using balanced $t$-wise independent distributions, which are related to orthogonal arrays of strength $t$. We contribute to the field by exploring such combinatorial designs in the context of the differential approximation measure. First, we establish a connection between the average differential ratio and orthogonal arrays. We deduce a differential approximation ratio of $1/q^k$ on $(k +1)$-partite instances of $\mathsf{k\,CSP\!-\!q}$, $\Omega(1/n^{k/2})$ on Boolean instances, $\Omega(1/n)$ when $k =2$, and $\Omega(1/\nu^{k -\lceil\log_{p^\kappa} k\rceil})$ when $k\geq 3$ and $q\geq 3$, where $p^\kappa$ is the smallest prime power greater than $q$. Secondly, by considering pairs of arrays related to balanced $k$-wise independence, we establish a reduction from $\mathsf{k\,CSP\!-\!q}$ to $\mathsf{k\,CSP\!-\!k}$ (with $q>k$), with an expansion factor of $1/(q -k/2)^k$ on the differential approximation guarantee. This, combined with the result of Yuri Nesterov, implies a lower differential approximability bound of $0.429/(q -1)^2$ for $\mathsf{2\,CSP\!-\!q}$. Finally, we prove that every Hamming ball with radius $k$ provides a $\Omega(1/n^k)$-approximation of the instance diameter. Our work highlights the utility of combinatorial designs in establishing approximation results., Comment: Preliminary versions of this work have been presented or published at the ISCO 2012, ISCO 2018 and IWOCA 2018 conferences

CO$_2$ is a major component of the icy mantles surrounding dust grains in planet and star formation regions. Understanding its photodesorption is crucial for explaining gas phase abundances in the coldest environments of the interstellar medium irradiated by vacuum-UV (VUV) photons. Photodesorption yields determined experimentally from CO$_2$ samples grown at low temperatures (T=15~K) have been found to be very sensitive to experimental methods and conditions. Several mechanisms have been suggested for explaining the desorption of CO$_2$, O$_2$ and CO from CO$_2$ ices. In the present study, the cross sections characterizing the dynamics of photodesorption as a function of photon fluence (determined from released molecules in the gas phase) and of ice composition modification (determined in situ in the solid phase) are compared for the first time for different photon flux conditions (from 7.3$\times 10^{12}$~photon/s/cm$^2$ to 2.2$\times 10^{14}$~photon/s/cm$^2$) using monochromatic synchrotron radiation in the VUV range (on the DESIRS beamline at SOLEIL). This approach reveals that CO and O$_2$ desorption are decorrelated from that of CO$_2$. CO and O$_2$ photodesorption yields depend on photon flux conditions and can be linked to surface chemistry. By contrast, the phodesorption yield of CO$_2$ is independent of the photon flux conditions and can be linked to bulk ice chemical modification, consistently with an indirect desorption induced by electronic transition (DIET) process., Comment: 11 pages, 7 figures

In [1], we introduced a family of combinatorial designs, which we call "alphabet reduction pairs of arrays", ARPAs for short. These designs depend on three integer parameters $q, p \leq q, k\leq p$: $q$ is the size of the symbol set $\{0, 1 ,\ldots, q -1\}$ in which the coefficients of the arrays take their values; $p$ is the maximum number of distinct symbols allowed in a row of the second array of the pair; $k$ is the larger integer for which the two arrays of the pair coincide -- up to the order of their rows -- on any $k$-ary subset of their columns. The first array must contain at least one occurrence of the word $0\ 1\ \ldots\ q -1$. Intuitively, the idea is to cover "as many as possible" occurrences of this word of $q$ symbols with "as few as possible" words of at most $p$ different symbols. These designs are related to the approximability of "Constraint Satisfaction Problems with bounded constraint arity", known as $k\,$CSPs. In this context, we are particularly interested in ARPAs in which the frequency of the word $0\ 1\ \ldots\ q -1$ is maximal. We introduce a seemingly simpler family of combinatorial designs as "Cover pairs of arrays" (CPAs). The arrays of a CPA take Boolean coefficients, and must still coincide on any $k$-ary subset of their columns. The purpose is, as it were, to cover "as many as possible" occurrences of the word of $q$ ones using "as few as possible" $q$-length Boolean words of weight at most $p$. We show that, when it comes to maximizing the frequency of the words $0\ 1\ \ldots\ q -1$ in ARPAs and $1\ 1\ \ldots\ 1$ in CPAs, ARPAs and CPAs are equivalent. We prove the optimality of the ARPAs given in [1] for the case $p =k$. In addition, we provide optimal ARPAs for the cases $k =1$ and $k =2$. We emphasize the fact that both families of combinatorial designs are related to the approximability of $k\,$CSPs.

Using GPU-accelerated state-vector emulation, we propose to embed a quantum computing ansatz into density-functional theory via density-based basis-set corrections (DBBSC) to obtain quantitative quantum-chemistry results on molecules that would otherwise require brute-force quantum calculations using hundreds of logical qubits. Indeed, accessing a quantitative description of chemical systems while minimizing quantum resources is an essential challenge given the limited qubit capabilities of current quantum processors. We provide a shortcut towards chemically accurate quantum computations by approaching the complete-basis-set limit through coupling the DBBSC approach, applied to any given variational ansatz, to an on-the-fly crafting of basis sets specifically adapted to a given system and user-defined qubit budget. The resulting approach self-consistently accelerates the basis-set convergence, improving electronic densities, ground-state energies, and first-order properties (e.g. dipole moments), but can also serve as a classical, a posteriori, energy correction to quantum hardware calculations with expected applications in drug design and materials science.

It is the position of the National Association of School Nurses (NASN) to actively support and protect the right of all students, including those with healthcare needs or disabilities, to participate in school-sponsored trips, in accordance with federal regulations (U.S. Department of Education [USDOE]/Office for Civil Rights [OCR], 2020). The registered professional school nurse (hereinafter referred to as school nurse) has the expertise necessary to prepare, plan, educate, and coordinate with staff, families, students, and other healthcare professionals in arranging equitable, inclusive, and appropriate care for all students to be able to attend school-sponsored trips. This statement provides the background and rationale for NASN's position. [This Position Statement was initially adopted in June 2013 and revised in January 2024.]

We present the first application to real molecular systems of the recently proposed linear-response theory for the density-based basis-set correction method [J. Chem. Phys. 158, 234107 (2023)]. We apply this approach to accelerate the basis-set convergence of excitation energies in the equation-of-motion coupled-cluster singles doubles (EOM-CCSD) method. We use an approximate linear-response framework which neglects the second-order derivative of the basis-set correction density functional and consists in simply adding to the usual Hamiltonian the one-electron potential generated by the first-order derivative of the functional. This additional basis-set correction potential is evaluated at the Hartree-Fock density, leading to a very computationally cheap basis-set correction. We tested this approach over a set of about 30 excitation energies computed for five small molecular systems and found that the excitation energies from the ground state to Rydberg states are the main source of basis-set error. These excitation energies systematically increase when the size of the basis set is increased, suggesting a biased description in favour of the excited state. Despite the simplicity of the present approach, the results obtained with the basis-set corrected EOM-CCSD method are encouraging as they yield to a mean absolute deviation of 0.02 eV for the aug-cc-pVTZ basis set, while it is of 0.04 eV using the standard EOM-CCSD method. This might open the path to an alternative to explicitly correlated approaches to accelerate the basis-set convergence of excitation energies.

This work reports an efficient density-fitting implementation of the density-based basis-set correction (DBBSC) method in the MOLPRO software. This method consists in correcting the energy calculated by a wave-function method with a given basis set by an adapted basis-set correction density functional incorporating the short-range electron correlation effects missing in the basis set, resulting in an accelerated convergence to the complete-basis-set limit. Different basis-set correction density-functional approximations are explored and the complementary-auxiliary-basis-set single-excitation correction is added. The method is tested on a benchmark set of reaction energies at the second-order M{\o}ller-Plesset (MP2) level and a comparison with the explicitly correlated MP2-F12 method is provided. The results show that the DBBSC method greatly accelerates the basis convergence of MP2 reaction energies, without reaching the accuracy of the MP2-F12 method but with a lower computational cost., Comment: 6 pages, 2 figures, 1 table

The basis-set correction method based on density-functional theory consists in correcting the energy calculated by a wave-function method with a given basis set by a density functional. This basis-set correction density functional incorporates the short-range electron correlation effects missing in the basis set. This results in accelerated basis convergences of ground-state energies to the complete-basis-set limit. In this work, we extend the basis-set correction method to a linear-response formalism for calculating excited-state energies. We give the general linear-response equations, as well as the more specific equations for configuration-interaction wave functions. As a proof of concept, we apply this approach to the calculations of excited-state energies in a one-dimensional two-electron model system with harmonic potential and a Dirac-delta electron-electron interaction. The results obtained with full-configuration-interaction wave functions expanded in a basis of Hermite functions and a local-density-approximation basis-set correction functional show that the present approach does not help in accelerating the basis convergence of excitation energies. However, we show that it significantly accelerates basis convergences of excited-state total energies.

We consider a one-dimensional effective quantum electrodynamics (QED) model of the relativistic hydrogen-like atom using delta-potential interactions. We discuss the general exact theory and the Hartree-Fock approximation. The present one-dimensional effective QED model shares the essential physical feature of the three-dimensional theory: the nuclear charge polarizes the vacuum state (creation of electron-positron pairs) which results in a QED Lamb-type shift of the bound-state energy. Yet, this 1D effective QED model eliminates some of the most serious technical difficulties of the three-dimensional theory coming from renormalization. We show how to calculate the vacuum-polarization density at zeroth order in the two-particle interaction and the QED Lamb-type shift of the bound-state energy at first order in the two-particle interaction. The present work may be considered as a step toward the development of a quantum-chemistry effective QED theory of atoms and molecules.

We develop semiclassical approximations for calculating photoabsorption cross sections beyond the continuum threshold in quantum many-body systems. These approximations use the fully quantum-mechanical Wigner function of the ground state and semiclassical expansions only for the part of the cross section depending on the continuum states, thus avoiding the difficult explicit calculation of the continuum states. Even though the approach is general, we test it in electronic-structure theory for the photoionization cross sections of the hydrogen and helium atoms. The results suggest that these semiclassical approximations can be used to obtain good estimates of cross sections at high energy.

Barium zirconate (BaZrO$_3$) is one of the very few perovskites that is claimed to retain an average cubic structure down to \SI{0}{\K}, while being energetically very close to an antiferrodistortive phase obtained by condensation of a soft phonon mode at the R point of the Brillouin zone boundary. In this work, we report a combined experimental and theoretical study of the temperature dependence of this soft phonon mode. Inelastic neutron and x-ray scattering measurements on single crystals show that it softens substantially from \SI{9.4}{\meV} at room temperature to \SI{5.6}{\meV} at \SI{2}{\K}. In contrast, the acoustic mode at the same R point is nearly temperature independent. The effect of the anharmonicity on the lattice dynamics is investigated non-perturbatively using direct dynamic simulations as well as a first-principles based self-consistent phonon theory, including quantum fluctuations of the atomic motion. By adding cubic and quartic anharmonic force constants, quantitative agreement with the neutron data for the temperature dependence of the antiferrodistortive mode is obtained. The quantum fluctuations of the atomic motion are found to be important to obtain the proper temperature dependence at low temperatures. The mean squared displacements of the different atoms are determined as function of temperature and are shown to be consistent with available experimental data. Adding anharmonicity to the computed fluctuations of the Ba-O distances also improves the comparison with available EXAFS data at \SI{300}{\K}.

We consider the calculations of photoionization spectra and core resonances of open-shell systems using range-separated time-dependent density-functional theory. Specifically, we use the time-dependent range-separated hybrid (TDRSH) scheme, combining a long-range Hartree-Fock (HF) exchange potential and kernel with a short-range potential and kernel from a local density-functional approximation, and the time-dependent locally range-separated hybrid (TDLRSH) scheme, which uses a local range-separation parameter. To efficiently perform the calculations, we formulate a spin-unrestricted linear-response Sternheimer approach in a non-orthogonal B-spline basis set and using appropriate frequency-dependent boundary conditions. We illustrate this approach on the Li atom, which suggests that TDRSH and TDLRSH are adequate simple methods for estimating single-electron photoionization spectra of open-shell systems., Comment: arXiv admin note: substantial text overlap with arXiv:2203.05951

The VB-QMC method is presented in this chapter. It consists of using in quantum Monte Carlo (QMC) approaches with a wave function expressed as a usually short expansion of classical Valence-Bond (VB) structures supplemented by a Jastrow factor to account for dynamical correlation. Two variants exist: the VB-VMC (using variational Monte Carlo) and VB-DMC (using diffusion Monte Carlo) methods. QMC algorithms circumvent the notorious non-orthogonality issue of classical VB approaches, and allow highly efficient calculations on massively parallel machines. Calculation of VB weights and resonance energies are possible at the VB-VMC level, which makes VB-VMC a correlated method retaining all the interpretative capabilities of classical VB methods. Several recent applications are shown to illustrate the potential of this method as a modern alternative to classical VB methods to study ground and excited states of molecules., Comment: book chapter submitted to Comprehensive Computational Chemistry

Public blockchains are decentralized networks where each participating node executes the same decision-making process. This form of decentralization does not scale well because the same data are stored on each network node, and because all nodes must validate each transaction prior to their confirmation. One solution approach decomposes the nodes of a blockchain network into subsets called "shards", each shard processing and storing disjoint sets of transactions in parallel. To fully benefit from the parallelism of sharded blockchains, the processing load of shards must be evenly distributed. However, the problem of computing balanced workloads is theoretically hard and further complicated in practice as transaction processing times are unknown prior to be assigned to shards. In this paper we introduce a dynamic workload-balancing algorithm where the allocation strategy of transactions to shards is periodically adapted based on the recent workload history of shards. Our algorithm is an adaptation to sharded blockchains of a consensus-based load-balancing algorithm. It is a fully distributed algorithm inline with network based applications such as blockchains. Some preliminary results are reported based on simulations that shard transactions of three well-known blockchain platforms.

In this paper, we revisit the high-pressure behavior of BaZrO3 by a combination of first-principle calculations, Raman spectroscopy, and x-ray diffraction under high pressure. We confirm experimentally the cubic-to-tetragonal transition at 10 GPa and find no evidence for any other phase transition up to 45 GPa, the highest pressures investigated, at variance with past reports. We re-investigate phase stability with density functional theory considering not only the known tetragonal (I4/mcm) phase but also other potential antiferrodistortive candidates. This shows that the tetragonal phase becomes progressively more stable upon increasing pressure as compared to phases with more complex tilt systems. The possibility for a second transition to another tilted phase at higher pressures, and in particular to the very common orthorhombic Pnma structure, is therefore ruled out.

We propose to compute physical properties by Monte Carlo calculations using conditional expectation values. The latter are obtained on top of the usual Monte Carlo sampling by partitioning the physical space in several subspaces or fragments, and subsampling each fragment (i.e., performing side walks) while freezing the environment. No bias is introduced and a zero-variance principle holds in the limit of separability, i.e. when the fragments are independent. In practice, the usual bottleneck of Monte Carlo calculations -- the scaling of the statistical fluctuations as a function of the number of particles N -- is relieved for extensive observables. We illustrate the method in variational Monte Carlo on the 2D Hubbard model and on metallic hydrogen chains using Jastrow-Slater wave functions. A factor O(N) is gained in numerical efficiency.

We explore the merits of linear-response range-separated time-dependent density-functional theory (TDDFT) for the calculation of photoionization spectra. We consider two variants of range-separated TDDFT, namely the timedependent range-separated hybrid (TDRSH) scheme which uses a global range-separation parameter and the timedependent locally range-separated hybrid (TDLRSH) which uses a local range-separation parameter, and compare with standard time-dependent local-density approximation (TDLDA) and time-dependent Hartree-Fock (TDHF). We show how to calculate photoionization spectra with these methods using the Sternheimer approach formulated in a nonorthogonal B-spline basis set with appropriate frequency-dependent boundary conditions. We illustrate these methods on the photoionization spectrum of the Be atom, focusing in particular on the core resonances. Both the TDRSH and TDLRSH photoionization spectra are found to constitute a large improvement over the TDLDA photoionization spectrum and a more modest improvement over the TDHF photoionization spectrum.

The present work proposes an approach to obtain a basis-set correction based on density-functional theory (DFT) for the computation of molecular properties in wave-function theory (WFT). This approach allows one to accelerate the basis-set convergence of any energy derivative of a non-variational WFT method, generalizing previous works on the DFT-based basis-set correction where either only ground-state energies could be computed with non-variational wave functions [J. Phys. Chem. Lett. 10, 2931 (2019)] or properties where computed as expectation values over variational wave functions [J. Chem. Phys. 155, 044109 (2021)]. This work focuses on the basis-set correction of dipole moments in coupled-cluster with single, double, and perturbative triple excitations (CCSD(T)), which is numerically tested on a set of fourteen molecules with dipole moments covering two orders of magnitude. As the basis-set correction relies only on Hartree-Fock densities, its computational cost is marginal with respect to the one of the CCSD(T) calculations. Statistical analysis of the numerical results shows a clear improvement of the basis convergence of the dipole moment with respect to the usual CCSD(T) calculations.

We reexamine the recently introduced basis-set correction theory based on density-functional theory consisting in correcting the basis-set incompleteness error of wave-function methods using a density functional. We use a one-dimensional model Hamiltonian with delta-potential interactions which has the advantage of making easier to perform a more systematic analysis than for three-dimensional Coulombic systems while keeping the essence of the slow basis convergence problem of wave-function methods. We provide some mathematical details about the theory and propose a new variant of basis-set correction which has the advantage of being suited to the development of an adapted local-density approximation. We show indeed how to develop a local-density approximation for the basis-set correction functional which is automatically adapted to the basis set employed, without resorting to range-separated density-functional theory as in previous works, but using instead a finite uniform electron gas whose electron-electron interaction is projected on the basis set. The work puts the basis-set correction theory on firmer grounds and provides an interesting strategy for the improvement of this approach.

The semistochastic heat-bath configuration interaction (SHCI) method is a selected configuration interaction plus perturbation theory method that has provided near-full configuration interaction (FCI) levels of accuracy for many systems with both single- and multi-reference character. However, obtaining accurate energies in the complete basis set limit is hindered by the slow convergence of the FCI energy with respect to basis size. Here we show that the recently developed basis-set correction method based on range-separated density-functional theory can be used to significantly speed up basis-set convergence in SHCI calculations. In particular, we study two such schemes that differ in the functional used, and apply them to transition metal atoms and monoxides to obtain total, ionization, and dissociation energies well converged to the complete-basis-set limit within chemical accuracy., Comment: 5 figures, supplemental material

Altair is the fastest rotating star at less than 10 parsecs from the Sun. Its precise modelling is a landmark for our understanding of stellar evolution with fast rotation, and all observational constraints are most welcome to better determine the fundamental parameters of this star. We wish to improve the seismic spectrum of Altair and confirm the $\delta$-Scuti nature of this star. We used the photometric data collected by the Microvariability and Oscillations of STars (MOST) satellite in the form of a series of Fabry images to derive Altair light curves at four epochs, namely in 2007, 2011, 2012, and 2013. We first confirm the presence of $\delta$-Scuti oscillations in the light curves of Altair. We extend the precision of some eigenfrequencies and add new ones to the spectrum of Altair, which now has 15 detected eigenmodes. The rotation period, which is expected at $\sim$7h46min from models reproducing interferometric data, seems to appear in the 2012 data set, but it still needs confirmation. Finally, Altair modal oscillations show noticeable amplitude variations on a timescale of 10 to 15 days, which may be the signature of a coupling between oscillations and thermal convection in the layer where the kappa-mechanism is operating.The Altair oscillation spectrum does not contain a large number of excited eigenmodes, which is similar to the fast rotating star HD220811. This supports the idea that fast rotation hinders the excitation of eigenmodes as already pointed out by theoretical investigations., Comment: 9 pages, 13 figures, to appear in Astronomy & Astrophysics

This work provides a self-consistent extension of the recently proposed density-based basis-set correction method for wave-function electronic-structure calculations [J. Chem. Phys. 149, 194301 (2018)]. In contrast to the previously used approximation where the basis-set correction density functional was a posteriori added to the energy from a wave-function calculation, here the energy minimization is performed including the basis-set correction. Compared to the non-self-consistent approximation, this allows one to lower the total energy and change the wave function under the effect of the basis-set correction. This work addresses two main questions: i) What is the change in total energy compared to the non-self-consistent approximation, and ii) can we obtain better properties, namely dipole moments, with the basis-set corrected wave functions? We implement the present formalism with two different basis-set correction functionals and test it on different molecular systems. The main results of the study are that i) the total energy lowering obtained by the self-consistent approach is extremely small, which justifies the use of the non-self-consistent approximation, and ii) the dipole moments obtained from the basis-set corrected wave functions are improved, being already close to their complete-basis-set values with triple-zeta basis sets. Thus, the present study further confirms the soundness of the density-based basis-set correction scheme.

In this chapter, we provide a review of ground-state Kohn-Sham density-functional theory of electronic systems and some of its extensions, we present exact expressions and constraints for the exchange and correlation density functionals, and we discuss the main families of approximations for the exchange-correlation energy: semilocal approximations, single-determinant hybrid approximations, multideterminant hybrid approximations, dispersion-corrected approximations, as well as orbital-dependent exchange-correlation density functionals. The chapter aims at providing both a consistent bird's-eye view of the field and a detailed description of some of the most used approximations. It is intended to be readable by chemists/physicists and applied mathematicians., Comment: Chapter in the book "Density Functional Theory" edited by Eric Canc\`es and Gero Friesecke

A relativistic density-functional theory based on a Fock-space effective quantum-electrodynamics (QED) Hamiltonian using the Coulomb or Coulomb-Breit two-particle interaction is developed. This effective QED theory properly includes the effects of vacuum polarization through the creation of electron-positron pairs but does not include explicitly the photon degrees of freedom. It is thus a more tractable alternative to full QED for atomic and molecular calculations. Using the constrained-search formalism, a Kohn-Sham scheme is formulated in a quite similar way to non-relativistic density-functional theory, and some exact properties of the involved density functionals are studied, namely charge-conjugation symmetry and uniform coordinate scaling. The usual no-pair Kohn-Sham scheme is obtained as a well-defined approximation to this relativistic density-functional theory., Comment: version 2, with a few improvements, to appear in SciPost Chemistry

Helium implantation in epitaxial thin films is a way to control the out-of-plane deformation independently from the in-plane strain controlled by epitaxy. In particular, implantation by means of a helium microscope allows for local implantation and patterning down to the nanometer resolution, which is of interest for device applications. We present here a study of bismuth ferrite (BiFeO3) films where strain was patterned locally by helium implantation. Our combined Raman, XRD and TEM study shows that the implantation causes an elongation of the BiFeO3 unit cell and ultimately a transition towards the so-called super-tetragonal polymorph via states with mixed phases. In addition, TEM reveals the onset of amorphization at a threshold dose that does not seem to impede the overall increase in tetragonality. The phase transition from the R-like to T-like BiFeO3 appears as first-order in character, with regions of phase coexistence and abrupt changes in lattice parameters.

We construct the complementary short-range correlation relativistic local-density-approximation functional to be used in relativistic range-separated density-functional theory based on a Dirac-Coulomb Hamiltonian in the no-pair approximation. For this, we perform relativistic random-phase-approximation calculations of the correlation energy of the relativistic homogeneous electron gas with a modified electron-electron interaction, we study the high-density behavior, and fit the results to a parametrized expression. The obtained functional should eventually be useful for electronic-structure calculations of strongly correlated systems containing heavy elements.

The Double Travelling Salesman Problem with Multiple Stacks, DTSPMS, deals with the collect and delivery of n commodities in two distinct cities, where the pickup and the delivery tours are related by LIFO constraints. During the pickup tour, commodities are loaded into a container of k rows, or stacks, with capacity c. This paper focuses on computational aspects of the DTSPMS, which is NP-hard. We first review the complexity of two critical subproblems: deciding whether a given pair of pickup and delivery tours is feasible and, given a loading plan, finding an optimal pair of pickup and delivery tours, are both polynomial under some conditions on k and c. We then prove a (3k)/2 standard approximation for the MinMetrickDTSPMS, where k is a universal constant, and other approximation results for various versions of the problem. We finally present a matching-based heuristic for the 2DTSPMS, which is a special case with k=2 rows, when the distances are symmetric. This yields a 1/2-o(1), 3/4-o(1) and 3/2+o(1) standard approximation for respectively Max2DTSPMS, its restriction Max2DTSPMS-(1,2) with distances 1 and 2, and Min2DTSPMS-(1,2), and a 1/2-o(1) differential approximation for Min2DTSPMS and Max2DTSPMS.

The recently developed semistochastic heat-bath configuration interaction (SHCI) method is a systematically improvable selected configuration interaction plus perturbation theory method capable of giving essentially exact energies for larger systems than is possible with other such methods. We compute SHCI atomization energies for 55 molecules which have been used as a test set in prior studies because their atomization energies are known from experiment. Basis sets from cc-pVDZ to cc-pV5Z are used, totaling up to 500 orbitals and a Hilbert space of $10^{32}$ Slater determinants for the largest molecules. For each basis, an extrapolated energy well within chemical accuracy (1 kcal/mol or 1.6 mHa/mol) of the exact energy for that basis is computed using only a tiny fraction of the entire Hilbert space. We also use our almost exact energies to benchmark coupled-cluster [CCSD(T)] energies. The energies are extrapolated to the complete basis set limit and compared to the experimental atomization energies. The extrapolations are done both without and with a basis-set correction based on density-functional theory. The mean absolute deviations from experiment for these extrapolations are 0.46 kcal/mol and 0.51 kcal/mol, respectively. Orbital optimization methods used to obtain improved convergence of the SHCI energies are also discussed., Comment: 6 figures, 4 supplemental files

We develop relativistic short-range exchange energy functionals for four-component relativistic range-separated density-functional theory using a Dirac-Coulomb Hamiltonian in the no-pair approximation. We show how to improve the short-range local-density approximation exchange functional for large range-separation parameters by using the on-top exchange pair density as a new variable. We also develop a relativistic short-range generalized-gradient approximation exchange functional which further increases the accuracy for small range-separation parameters. Tests on the helium, beryllium, neon, and argon isoelectronic series up to high nuclear charges show that this latter functional gives exchange energies with a maximal relative percentage error of 3 %. The development of this exchange functional represents a step forward for the application of four-component relativistic range-separated density-functional theory to chemical compounds with heavy elements., Comment: Journal of Chemical Physics, American Institute of Physics, In press

We extend to strongly correlated molecular systems the recently introduced basis-set incompleteness correction based on density-functional theory (DFT) [E. Giner et al., J. Chem. Phys. 149, 194301 (2018)]. This basis-set correction relies on a mapping between wave-function calculations in a finite basis set and range-separated DFT (RSDFT) through the definition of an effective non-divergent interaction corresponding to the electron-electron Coulomb interaction projected in the finite basis set. This enables the use of RSDFT-type complementary density functionals to recover the dominant part of the short-range correlation effects missing in this finite basis set. To study both weak and strong correlation regimes we consider the potential energy curves of the H10, N2, O2, and F2 molecules up to the dissociation limit, and we explore various approximations of complementary functionals fulfilling two key properties: spin-multiplet degeneracy (i.e., independence of the energy with respect to the spin projection Sz) and size consistency. Specifically, we investigate the dependence of the functional on different types of on-top pair densities and spin polarizations. The key result of this study is that the explicit dependence on the on-top pair density allows one to completely remove the dependence on any form of spin polarization without any significant loss of accuracy. Quantitatively, we show that the basis-set correction reaches chemical accuracy on atomization energies with triple-zeta quality basis sets for most of the systems studied here. Also, the present basis-set incompleteness correction provides smooth potential energy curves along the whole range of internuclear distances.

We extend the range-separated double-hybrid RSH+MP2 method [J. G. Angyan et al., Phys. Rev. A 72, 012510 (2005)], combining long-range HF exchange and MP2 correlation with a short-range density functional, to a fully self-consistent version using the optimized-effective-potential technique in which the orbitals are obtained from a local potential including the long-range HF and MP2 contributions. We test this approach, that we name RS-OEP2, on a set of small closed-shell atoms and molecules. For the commonly used value of the range-separation parameter $\mu=0.5$ bohr$^{-1}$, we find that self-consistency does not seem to bring any improvement for total energies, ionization potentials, and electronic affinities. However, contrary to the non-self-consistent RSH+MP2 method, the present RS-OEP2 method gives a LUMO energy which physically corresponds to a neutral excitation energy and gives local exchange-correlation potentials which are reasonably good approximations to the corresponding Kohn-Sham quantities. At a finer scale, we find that RS-OEP2 gives largely inaccurate correlation potentials and correlated densities, which points to the need of further improvement of this type of range-separated double hybrids., Comment: Journal of Chemical Theory and Computation, American Chemical Society, In press

Similar to other electron correlation methods, many-body perturbation theory methods based on Green functions, such as the so-called $GW$ approximation, suffer from the usual slow convergence of energetic properties with respect to the size of the one-electron basis set. This displeasing feature is due to lack of explicit electron-electron terms modeling the infamous Kato electron-electron cusp and the correlation Coulomb hole around it. Here, we propose a computationally efficient density-based basis set correction based on short-range correlation density functionals which significantly speeds up the convergence of energetics towards the complete basis set limit. The performance of this density-based correction is illustrated by computing the ionization potentials of the twenty smallest atoms and molecules of the GW100 test set at the perturbative $GW$ (or $G_0W_0$) level using increasingly large basis sets. We also compute the ionization potentials of the five canonical nucleobases (adenine, cytosine, thymine, guanine, and uracil) and show that, here again, a significant improvement is obtained., Comment: 11 pages, 2 figures (supporting information available)

Model materials are precious test cases for elementary theories and provide building blocks for the understanding of more complex cases. Here, we describe the lattice dynamics of the structural phase transition in francisite Cu3Bi(SeO3)2O2Cl at 115 K and show that it provides a rare archetype of a transition driven by a soft antipolar phonon mode. In the high-symmetry phase at hightemperatures, the soft mode is found at (0,0,0.5) at the Brillouin zone boundary and is measured by inelastic X-ray scattering and thermal diffuse scattering. In the low-symmetry phase, this softmode is folded back onto the center of the Brillouin zone as a result of the doubling of the unit cell, and appears as a fully symmetric mode that can be tracked by Raman spectroscopy. On both sides of the transition, the mode energy squared follows a linear behaviour over a large temperature range. First-principles calculations reveal that, surprisingly, the flat phonon band calculated for the high-symmetry phase seems incompatible with the displacive character found experimentally. We discuss this unusual behavior in the context of an ideal Kittel model of an antiferroelectric transition., Comment: 6 pages, 5 Figures, supplemental material included

BaZrO3 is a perovskite that remains in the simple cubic phase at all temperatures, hence with no first-order Raman-active phonon mode allowed by symmetry. Yet, it exhibits an intense Raman spectrum with sharp and well-defined features. Here, we report the evolution of the Raman spectrum of BaZrO3 single crystals in a broad temperature range (4--1200~K) and discuss its origin with the support of detailed first-principle calculations of the lattice dynamics. Phonon calculations are performed not only for the cubic phase of BaZrO3, but also for the low-symmetry phases with octahedra tilts that have been suspected to exist at the nanoscale. We show that the Raman spectrum shows no direct evidence for these nanodomains, but can instead be explained by classical second-order Raman scattering. We provide an assignment of the dominant features to phonon modes combinations. In particular, we show that the high frequency range of the spectrum is dominated by overtones and shows an image of the phonon density of states corresponding to the stretching modes of the oxygen octahedra.

By combining extrapolated selected configuration interaction (sCI) energies obtained with the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) algorithm with the recently proposed short-range density-functional correction for basis-set incompleteness [Giner et al.,J. Chem. Phys. 2018, 149, 194301], we show that one can get chemically accurate vertical and adiabatic excitation energies with, typically, augmented double-$\zeta$ basis sets. We illustrate the present approach on various types of excited states (valence, Rydberg, and double excitations) in several small organic molecules (methylene, water, ammonia, carbon dimer and ethylene). The present study clearly evidences that special care has to be taken with very diffuse excited states where the present correction does not catch the radial incompleteness of the one-electron basis set., Comment: 9 pages, 6 figures

We construct range-separated double-hybrid schemes which combine coupled-cluster or random-phase approximations with a density functional based on a two-parameter Coulomb-attenuating-method-like decomposition of the electron-electron interaction. We find that the addition of a fraction of short-range electron-electron interaction in the wave-function part of the calculation is globally beneficial for the range-separated double-hybrid scheme involving a variant of the random-phase approximation with exchange terms. Even though the latter scheme is globally as accurate as the corresponding scheme employing only second-order M{{\o}}ller-Plesset perturbation theory for atomization energies, reaction barrier heights, and weak intermolecular interactions of small molecules, it is more accurate for the more complicated case of the benzene dimer in the stacked configuration. The present range-separated double-hybrid scheme employing a random-phase approximation thus represents a new member in the family of double hybrids with minimal empiricism which could be useful for general chemical applications., Comment: arXiv admin note: text overlap with arXiv:1804.03373

We report a universal density-based basis-set incompleteness correction that can be applied to any wave function method. The present correction, which appropriately vanishes in the complete basis set (CBS) limit, relies on short-range correlation density functionals (with multi-determinant reference) from range-separated density-functional theory (RS-DFT) to estimate the basis-set incompleteness error. Contrary to conventional RS-DFT schemes which require an \textit{ad hoc} range-separation \textit{parameter} $\mu$, the key ingredient here is a range-separation \textit{function} $\mu(\bf{r})$ that automatically adapts to the spatial non-homogeneity of the basis-set incompleteness error. As illustrative examples, we show how this density-based correction allows us to obtain CCSD(T) atomization and correlation energies near the CBS limit for the G2 set of molecules with compact Gaussian basis sets., Comment: 7 pages, 2 figures and 1 table (supporting information available)

We investigate the performance of the range-separated hybrid (RSH) scheme, which combines long-range Hartree-Fock (HF) and a short-range density-functional approximation (DFA), for calculating photoexcitation/photoionization spectra of the H and He atoms, using a B-spline basis set in order to correctly describe the continuum part of the spectra. The study of these simple systems allows us to quantify the influence on the spectra of the errors coming from the short-range exchange-correlation DFA and from the missing long-range correlation in the RSH scheme. We study the differences between using the long-range HF exchange (nonlocal) potential and the long-range exact exchange (local) potential. Contrary to the former, the latter supports a series of Rydberg states and gives reasonable photoexcitation/photoionization spectra, even without applying linear-response theory. The most accurate spectra are obtained with the linear-response time-dependent range-separated hybrid (TDRSH) scheme. In particular, for the He atom at the optimal value of the range-separation parameter, TDRSH gives slightly more accurate photoexcitation and photoioniza-tion spectra than standard linear-response time-dependent HF. More generally, the present work shows the potential of range-separated density-functional theory for calculating linear and nonlinear optical properties involving continuum states.

\textsc{Quantum Package} is an open-source programming environment for quantum chemistry specially designed for wave function methods. Its main goal is the development of determinant-driven selected configuration interaction (sCI) methods and multi-reference second-order perturbation theory (PT2). The determinant-driven framework allows the programmer to include any arbitrary set of determinants in the reference space, hence providing greater methodological freedoms. The sCI method implemented in \textsc{Quantum Package} is based on the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) algorithm which complements the variational sCI energy with a PT2 correction. Additional external plugins have been recently added to perform calculations with multireference coupled cluster theory and range-separated density-functional theory. All the programs are developed with the IRPF90 code generator, which simplifies collaborative work and the development of new features. \textsc{Quantum Package} strives to allow easy implementation and experimentation of new methods, while making parallel computation as simple and efficient as possible on modern supercomputer architectures. Currently, the code enables, routinely, to realize runs on roughly 2\,000 CPU cores, with tens of millions of determinants in the reference space. Moreover, we have been able to push up to 12\,288 cores in order to test its parallel efficiency. In the present manuscript, we also introduce some key new developments: i) a renormalized second-order perturbative correction for efficient extrapolation to the full CI limit, and ii) a stochastic version of the CIPSI selection performed simultaneously to the PT2 calculation at no extra cost., Comment: 21 pages, 8 figures

We introduce an approximation to the short-range correlation energy functional with multide-terminantal reference involved in a variant of range-separated density-functional theory. This approximation is a local functional of the density, the density gradient, and the on-top pair density, which locally interpolates between the standard Perdew-Burke-Ernzerhof correlation functional at vanishing range-separation parameter and the known exact asymptotic expansion at large range-separation parameter. When combined with (selected) configuration-interaction calculations for the long-range wave function, this approximation gives accurate dissociation energy curves of the H2, Li2, and Be2 molecules, and thus appears as a promising way to accurately account for static correlation in range-separated density-functional theory., Comment: The Journal of Chemical Physics, In press