Ruddlesden-Popper Oxyfluorides La2Ni1-xCuxO3F2 (0 ≤ x ≤ 1): Impact of the Ni/Cu Ratio on the Thermal Stability and Magnetic Properties.

Ruddlesden-Popper oxyfluorides of the substitution series La2Ni1-xCuxO3F2 (0 ≤ x ≤ 1) were obtained by topochemical fluorination with polyvinylidene fluoride (PVDF) of oxide precursors La2Ni1-xCuxO4. The thermal stability and the temperature-dependent unit cell evolution of the oxyfluorides were investigated by high-temperature XRD measurements. The oxyfluoride with x = 0.6 shows the highest decomposition temperature of θdec ∼ 520 °C, which is significantly higher than the ones found for the end members La2NiO3F2 (x = 0) θdec ∼ 460 °C and La2CuO3F2 (x = 1) θdec ∼ 430 °C. The magnetic properties of all La2Ni1-xCuxO3F2 oxyfluorides were characterized by field- and temperature-dependent measurements as well as DFT calculations of the magnetic ground state. An antiferromagnetic ordering was derived for all substitution levels. For the Néel temperature (TN), a nonlinear dependence on the copper content was found, and comparably high values of TN in the region of 200-250 K were observed in the broad composition range of 0.3 ≤ x ≤ 0.8.

Ruddlesden-Popper Oxyfluorides La2Ni1-xCuxO3F2 (0 ≤ x ≤ 1): Impact of the Ni/Cu Ratio on the Structure.

Ruddlesden-Popper oxyfluorides La2Ni1-xCuxO3F2 (0 ≤ x ≤ 1) were obtained by topochemical reaction of oxide precursors La2Ni1-xCuxO4, prepared by citrate-based soft chemistry synthesis, with polyvinylidene fluoride (PVDF) as the fluorine source. Systematic changes of the crystal structure in the oxide as well as the oxyfluoride substitution series were investigated. For 0.2 ≤ x ≤ 0.9, the oxyfluorides adopt the monoclinic (C2/c) structural distortion previously solved for the x = 0.8 compound based on neutron powder diffraction data, whereas the sample with a lower Cu content of x = 0.1 crystallizes in the orthorhombic (Cccm) structure variant of La2NiO3F2. The orthorhombic-to-monoclinic structural transition was found to be the result of an additional tilt component of the Jahn-Teller elongated CuO4F2 octahedra. The structural transitions were additionally studied by DFT calculations, confirming the monoclinic space group symmetry. The "channel-like" anionic ordering of the endmembers La2NiO3F2 and La2CuO3F2 was checked by 19F MAS NMR experiments and was found to persist throughout the entire substitution series.

Reproducibility of density functional approximations: How new functionals should be reported.

Density functional theory is the workhorse of chemistry and materials science, and novel density functional approximations are published every year. To become available in program packages, the novel density functional approximations (DFAs) need to be (re)implemented. However, according to our experience as developers of Libxc [Lehtola et al., SoftwareX 7, 1 (2018)], a constant problem in this task is verification due to the lack of reliable reference data. As we discuss in this work, this lack has led to several non-equivalent implementations of functionals such as Becke-Perdew 1986, Perdew-Wang 1991, Perdew-Burke-Ernzerhof, and Becke's three-parameter hybrid functional with Lee-Yang-Parr correlation across various program packages, yielding different total energies. Through careful verification, we have also found many issues with incorrect functional forms in recent DFAs. The goal of this work is to ensure the reproducibility of DFAs. DFAs must be verifiable in order to prevent the reappearance of the above-mentioned errors and incompatibilities. A common framework for verification and testing is, therefore, needed. We suggest several ways in which reference energies can be produced with free and open source software, either with non-self-consistent calculations with tabulated atomic densities or via self-consistent calculations with various program packages. The employed numerical parameters-especially the quadrature grid-need to be converged to guarantee a ≲0.1 µEh precision in the total energy, which is nowadays routinely achievable in fully numerical calculations. Moreover, as such sub-µEh level agreement can only be achieved when fully equivalent implementations of the DFA are used, the source code of the reference implementation should also be made available in any publication describing a new DFA.

Many recent density functionals are numerically ill-behaved.

Most computational studies in chemistry and materials science are based on the use of density functional theory. Although the exact density functional is unknown, several density functional approximations (DFAs) offer a good balance of affordable computational cost and semi-quantitative accuracy for applications. The development of DFAs still continues on many fronts, and several new DFAs aiming for improved accuracy are published every year. However, the numerical behavior of these DFAs is an often-overlooked problem. In this work, we look at all 592 DFAs for three-dimensional systems available in Libxc 5.2.2 and examine the convergence of the density functional total energy based on tabulated atomic Hartree-Fock wave functions. We show that several recent DFAs, including the celebrated SCAN family of functionals, show impractically slow convergence with typically used numerical quadrature schemes, making these functionals unsuitable both for routine applications and high-precision studies, as thousands of radial quadrature points may be required to achieve sub-µEh accurate total energies for these functionals, while standard quadrature grids like the SG-3 grid only contain O(100) radial quadrature points. These results are both a warning to users to always check the sufficiency of the quadrature grid when adopting novel functionals, as well as a guideline to the theory community to develop better-behaved density functionals.

Efficient and improved prediction of the band offsets at semiconductor heterojunctions from meta-GGA density functionals: A benchmark study.

Accurate theoretical prediction of the band offsets at interfaces of semiconductor heterostructures can often be quite challenging. Although density functional theory has been reasonably successful to carry out such calculations, efficient, accurate semilocal functionals are desirable to reduce the computational cost. In general, the semilocal functionals based on the generalized gradient approximation (GGA) significantly underestimate the bulk bandgaps. This, in turn, results in inaccurate estimates of the band offsets at the heterointerfaces. In this paper, we investigate the performance of several advanced meta-GGA functionals in the computational prediction of band offsets at semiconductor heterojunctions. In particular, we investigate the performance of r2SCAN (two times revised strongly constrained and appropriately normed functional), rMGGAC (revised semilocal functional based on cuspless hydrogen model and Pauli kinetic energy density functional), mTASK (modified Aschebrock and Kümmel meta-GGA functional), and local modified Becke-Johnson exchange-correlation functionals. Our results strongly suggest that these meta-GGA functionals for supercell calculations perform quite well, especially, when compared to computationally more demanding GW calculations. We also present band offsets calculated using ionization potentials and electron affinities, as well as band alignment via the branch point energies. Overall, our study shows that the aforementioned meta-GGA functionals can be used within the density functional theory framework to estimate the band offsets in semiconductor heterostructures with predictive accuracy.

Structure, Magnetism, and Thermal Stability of La2NiO2.5F3: A Ruddlesden-Popper Oxyfluoride Crystallizing in Space Group P42/nnm.

We report on the new Ruddlesden-Popper (RP) oxyfluoride La2NiO2.5F3 containing an unprecedented high amount of fluorine and Ni2+. This oxyfluoride was prepared by topochemical low-temperature fluorination of La2NiO4, which was obtained by a soft chemistry synthesis, with poly(vinylidene difluoride) (PVDF) as fluorinating agent. La2NiO2.5F3 is the first n = 1 RP compound crystallizing in the tetragonal space group P42/nnm (a = 5.7297(6) Å and c = 13.0106(2) Å). The crystal structure shows a unique tilting scheme of the NiO4F2 octahedra that has so far been only theoretically predicted. Combined neutron and X-ray powder diffraction experiments together with bond-valence-sum and DFT+U calculations reveal an unusual anion ordering with fluoride being located on the apical anion sites of the NiO4F2 octahedra. Excess fluorine ions were found to populate two of the four interstitial anion sites in an ordered fashion. A third interstitial anion position is occupied by oxygen ions while the fourth site remains unoccupied. This hitherto unobserved ordering scenario in RP oxyfluorides promotes a strong layerwise alternating tilting of the NiO4F2 octahedra. Magnetic measurements show strong antiferromagnetic interactions with a high Néel temperature of about 225 K and a pronounced ZFC/FC splitting most likely as the result of a small ferromagnetic moment arising from spin canting. The electronic structure was characterized by DFT and UV-vis spectroscopy, and a strong increase of Eg was found compared to La2NiO4 (3.4 eV vs 1.3 eV).

Bishop's hat silicene: a planar square silicon bilayer decorated with adatoms.

We investigate a family of free-standing quasi-two-dimensional silicon structures based on a planar square bilayer with adatom decorations. When attached to the bilayer, these adatoms form local reconstructions which resemble either a bishop's hat or elongated square bipyramids. We systematically constructed members of this family via exhaustive enumeration and then studied them using tight-binding and density-functional theory. We find that this geometry contributes significantly to the stability of the resulting structures, with some squared bilayers energetically more stable than the honeycomb bilayers. The most interesting phases were then characterized in more detail, and they all turned out metallic. Finally, we propose the [100] surface of ZrO2 as the most suitable substrate for the synthesis of these two-dimensional phases.

Halogen molecular modifications at high pressure: the case of iodine.

Metallization and dissociation are key transformations in diatomic molecules at high densities particularly significant for modeling giant planets. Using X-ray absorption spectroscopy and atomistic modeling, we demonstrate that in halogens, the formation of a connected molecular structure takes place at pressures well below metallization. Here we show that the iodine diatomic molecule first elongates by ∼0.007 Å up to a critical pressure of Pc ∼ 7 GPa, developing bonds between molecules. Then its length continuously decreases with pressure up to 15-20 GPa. Universal trends in halogens are shown and allow us to predict for chlorine a pressure of 42 ± 8 GPa for molecular bond-length reversal. Our findings contribute to tackling the molecule invariability paradigm in diatomic molecular phases at high pressures and may be generalized to other abundant diatomic molecules in the universe, including hydrogen.

First-Principles Identification of Single Photon Emitters Based on Carbon Clusters in Hexagonal Boron Nitride.

A recent study associates carbon with single photon emitters (SPEs) in hexagonal boron nitride (h-BN). This observation, together with the high mobility of carbon in h-BN, suggests the existence of SPEs based on carbon clusters. Here, by means of density functional theory calculations, we studied clusters of substitutional carbon atoms up to tetramers in h-BN. Two different conformations of neutral carbon trimers have zero-point line energies and shifts of the phonon sideband compatible with typical photoluminescence spectra. Moreover, some conformations of two small C clusters next to each other result in photoluminescence spectra similar to those found in the experiments. We also showed that vacancies are unable to reproduce the typical features of the phonon sideband observed in most measurements because of the large spectral weight of low-energy breathing modes, ubiquitous in such defects.

Bandgap of two-dimensional materials: Thorough assessment of modern exchange-correlation functionals.

The density-functional theory (DFT) approximations that are the most accurate for the calculation of bandgap of bulk materials are hybrid functionals, such as HSE06, the modified Becke-Johnson (MBJ) potential, and the GLLB-SC potential. More recently, generalized gradient approximations (GGAs), such as HLE16, or meta-GGAs, such as (m)TASK, have also proven to be quite accurate for the bandgap. Here, the focus is on two-dimensional (2D) materials and the goal is to provide a broad overview of the performance of DFT functionals by considering a large test set of 298 2D systems. The present work is an extension of our recent studies [T. Rauch, M. A. L. Marques, and S. Botti, Phys. Rev. B 101, 245163 (2020); Patra et al., J. Phys. Chem. C 125, 11206 (2021)]. Due to the lack of experimental results for the bandgap of 2D systems, G0W0 results were taken as reference. It is shown that the GLLB-SC potential and mTASK functional provide the bandgaps that are the closest to G0W0. Following closely, the local MBJ potential has a pretty good accuracy that is similar to the accuracy of the more expensive hybrid functional HSE06.

Noninvasive neuromodulation of the prefrontal cortex in young women with obesity: a randomized clinical trial.

BACKGROUND/OBJECTIVES: Obesity is associated with reduced neurocognitive performance. Individuals with obesity show decreased activation in the left dorsolateral prefrontal cortex (DLPFC), a key brain region relevant to the regulation of eating behavior. Transcranial direct current stimulation (tDCS) has emerged as a potential technique to correct these abnormalities. However, there is limited information to date, particularly in clinical settings and regarding long-term effects of tDCS. This study aimed to investigate the effects of DLPFC-targeted tDCS in young women with obesity. SUBJECT/METHODS: Randomized, double-blind, sham-controlled parallel-design clinical trial conducted in 38 women, aged 20-40 years, with BMI 30-35 kg/m2. STUDY DESIGN: Phase I: target engagement (immediate effects of tDCS on working memory performance), Phase II: tDCS only (ten sessions, 2 weeks), Phase III: tDCS + hypocaloric diet (six sessions, 30% energy intake reduction, 2 weeks, inpatient), Phase IV: follow-up at 1, 3, and 6 months. PRIMARY OUTCOME: change in body weight. SECONDARY OUTCOMES: change in eating behavior and appetite. Additional analyses: effect of Catechol-O-methyl transferase (COMT) gene variability. Data were analyzed as linear mixed models. RESULTS: There was no group difference in change in body weight during the tDCS intervention. At follow-up, the active group lost less weight than the sham group. In addition, the active group regained weight at 6-month follow-up, compared with sham. Genetic analysis indicated that COMT Met noncarriers were the subgroup that accounted for this paradoxical response in the active group. CONCLUSION: Our results suggest that in young women with class I obesity, tDCS targeted to the DLPFC does not facilitate weight loss. Indeed, we found indications that tDCS could have a paradoxical effect in this population, possibly connected with individual differences in dopamine availability. Future studies are needed to confirm these findings.

Reduced Density Matrix Functional Theory for Bosons.

Based on a generalization of Hohenberg-Kohn's theorem, we propose a ground state theory for bosonic quantum systems. Since it involves the one-particle reduced density matrix γ as a variable but still recovers quantum correlations in an exact way it is particularly well suited for the accurate description of Bose-Einstein condensates. As a proof of principle we study the building block of optical lattices. The solution of the underlying v-representability problem is found and its peculiar form identifies the constrained search formalism as the ideal starting point for constructing accurate functional approximations: The exact functionals F[γ] for this N-boson Hubbard dimer and general Bogoliubov-approximated systems are determined. For Bose-Einstein condensates with N_{BEC}≈N condensed bosons, the respective gradient forces are found to diverge, ∇_{γ}F∝1/sqrt[1-N_{BEC}/N], providing a comprehensive explanation for the absence of complete condensation in nature.

Curcuma supplementation in high-fat-fed C57BL/6 mice: no beneficial effect on lipid and glucose profile or prevention of weight gain.

PURPOSE: This experimental study investigated the effects of curcuma supplementation on weight gain, Body Adiposity Index, glucose and lipid profile, and liver and pancreas histology in C57BL/6 mice fed with a high-fat diet. METHODS: 40 animals were separated into four groups: standard diet (SD), standard diet plus curcuma (SD + C), high-fat diet (HFD), and high-fat diet plus curcuma (HFD + C). Curcuma dose was 8 mg/animal/day. Histological and biochemical analyses were performed at the end of the experimental period. RESULTS: Curcuma prevented weight gain, despite a higher food intake, and increased brown adipose tissue weight only in mice receiving standard diet. However, these changes were not observed in HFD + C group. The groups that received curcuma (SD + C and HFD + C) showed a pancreas with diffuse macro- and microgoticular steatosis. CONCLUSIONS: Curcuma supplementation did not prevent weight gain or improved glucose and lipid profile in mice receiving high-fat diet. Furthermore, there was evidence of possible curcuma toxicity in the pancreas of C57BL/6 mice. The implications of these findings on humans still need to be investigated.

Novel two-dimensional silicon-carbon binaries by crystal structure prediction.

The semimetallic bandstructure of graphene and silicene limit their use in functional devices. Mixing silicon and carbon offers a rather unexplored pathway to build semiconducting sheets compatible with current Si-based electronics. We present here a complete theoretical study of the phase diagram of two-dimensional silicon-carbon binaries. To scan the composition range, we employ an ab initio global structural prediction method, complemented by exhaustive enumeration of two-dimensional structure prototypes. We find a wealth of two-dimensional low-energy structures, from standard honeycomb single- and double-layers, passing by dumbbell geometries, to carbon nanosheets bridged by Si atoms. Many of these phases depart from planarity, either through buckling, or by germinating three-dimensional networks with a mixture of sp2 and sp3 bonds. We further characterize the most interesting crystal structures, unveiling a large variety of electronic properties, that could be exploited to develop high-performance electronic devices at the nanoscale.

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems.

Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light-matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).

The CECAM electronic structure library and the modular software development paradigm.

First-principles electronic structure calculations are now accessible to a very large community of users across many disciplines, thanks to many successful software packages, some of which are described in this special issue. The traditional coding paradigm for such packages is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, essentially from the compiler up, possibly with the exception of linear-algebra and message-passing libraries. This model has endured and been quite successful for decades. The successful evolution of the electronic structure methodology itself, however, has resulted in an increasing complexity and an ever longer list of features expected within all software packages, which implies a growing amount of replication between different packages, not only in the initial coding but, more importantly, every time a code needs to be re-engineered to adapt to the evolution of computer hardware architecture. The Electronic Structure Library (ESL) was initiated by CECAM (the European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure codes and redesign them as open-source libraries available to everybody. Such libraries include "heavy-duty" ones that have the potential for a high degree of parallelization and adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by, e.g., physicists and chemists when implementing new ideas. We envisage that this modular paradigm will improve overall coding efficiency and enable specialists (whether they be computer scientists or computational scientists) to use their skills more effectively and will lead to a more dynamic evolution of software in the community as well as lower barriers to entry for new developers. The model comes with new challenges, though. The building and compilation of a code based on many interdependent libraries (and their versions) is a much more complex task than that of a code delivered in a single self-contained package. Here, we describe the state of the ESL, the different libraries it now contains, the short- and mid-term plans for further libraries, and the way the new challenges are faced. The ESL is a community initiative into which several pre-existing codes and their developers have contributed with their software and efforts, from which several codes are already benefiting, and which remains open to the community.

Computational acceleration of prospective dopant discovery in cuprous iodide.

The zinc blende (γ) phase of copper iodide holds the record hole conductivity for intrinsic transparent p-type semiconductors. In this work, we employ a high-throughput approach to systematically explore strategies for enhancing γ-CuI further by impurity incorporation. Our objectives are not only to find a practical approach to increase the hole conductivity in CuI thin films, but also to explore the possibility for ambivalent doping. In total 64 chemical elements were investigated as possible substitutionals on either the copper or the iodine site. All chalcogen elements were found to display acceptor character when substituting iodine, with sulfur and selenium significantly enhancing carrier concentrations produced by the native VCu defects under conditions most favorable for impurity incorporation. Furthermore, eight impurities suitable for n-type doping were discovered. Unfortunately, our work also reveals that donor doping is hindered by compensating native defects, making ambipolar doping unlikely. Finally, we investigated how the presence of impurities influences the optical properties. In the majority of the interesting cases, we found no deep states in the band-gap, showing that CuI remains transparent upon doping.

Neural network force fields for simple metals and semiconductors: construction and application to the calculation of phonons and melting temperatures.

We present a practical procedure to obtain reliable and unbiased neural network based force fields for solids. Training and test sets are efficiently generated from global structural prediction runs, at the same time assuring the structural variety and importance of sampling the relevant regions of phase space. The neural networks are trained to yield not only good formation energies, but also accurate forces and stresses, which are the quantities of interest for molecular dynamics simulations. Finally, we construct, as an example, several force fields for both semiconducting and metallic elements, and prove their accuracy for a variety of structural and dynamical properties. These are then used to study the melting of bulk copper and gold.

On the time evolution of fermionic occupation numbers.

We derive an approximate equation for the time evolution of the natural occupation numbers for fermionic systems. The evolution of such numbers is connected with the symmetry-adapted generalized Pauli exclusion principle, as well as with the evolution of the natural orbitals and a set of many-body relative phases. We then relate the evolution of these phases to a geometrical and a dynamical term attached to some of the Slater determinants appearing in the configuration-interaction expansion of the wave function. Our approach becomes exact for highly symmetric systems whenever the wave function possesses as many Slater determinants as independent occupation numbers.

On the calculation of the bandgap of periodic solids with MGGA functionals using the total energy.

During the last few years, it has become more and more clear that functionals of the meta generalized gradient approximation (MGGA) are more accurate than GGA functionals for the geometry and energetics of electronic systems. However, MGGA functionals are also potentially more interesting for the electronic structure, in particular, when the potential is nonmultiplicative (i.e., when MGGAs are implemented in the generalized Kohn-Sham framework), which may help to get more accurate bandgaps. Here, we show that the calculation of bandgap of solids with MGGA functionals can also be done very accurately in a non-self-consistent manner. This scheme uses only the total energy and can, therefore, be very useful when the self-consistent implementation of a particular MGGA functional is not available. Since self-consistent MGGA calculations may be difficult to converge, the non-self-consistent scheme may also help to speed up the calculations. Furthermore, it can be applied to any other types of functionals, for which the implementation of the corresponding potential is not trivial.