Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT Calculations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Bayesian automated weighting of aggregated DFT, MD, and experimental data for candidate thermodynamic models of aluminum with uncertainty quantification

Atomic-scale modeling methods such as density functional theory (DFT) and molecular dynamics (MD) can predict the thermodynamic properties of materials at a lower cost than experimental measurements. However, their regular usage in thermodynamic model construction is hampered by the lack of quantitative agreement with experimental measurements and the lack of uncertainty estimates on the data. To make regular usage of this atomistic simulation data, it is important to assess whether the atomistic simulation datasets, by themselves or in combination with experimental measurements, result in the same physics-informed models best supported by experimental measurements alone. Here, models of aluminum thermodynamic properties are discussed using three data sources: atomistic calculations (DFT and MD), experiments, and a combination of atomistic calculations and experiments. The study shows that, after ensuring self-consistency in predicting key invariant points, both experimental measurements and atomistic calculations can significantly contribute to an optimal model.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fundamental Studies of the Vibrational, Electronic, and Photophysical Properties of Tetrapyrrolic Architectures

The ability to capture and utilize light in the near-ultraviolet (NUV), visible and near-infrared (NIR-I and NIR-II) spectral regions (i.e., 320–400, 400–700, 700–1000, 1000–1700 nm) is essential for any solar-energy conversion scheme. Nature employs chlorophylls and bacteriochlorophylls in light-harvesting architectures to absorb light in the blue and red/NIR regions. Accessory pigments (carotenoids, bilins) augment absorption of the (bacterio)chlorophylls in the green region. The harvested energy is funneled to a reaction center protein, where charge separation occurs. Subsequent migration of the electron and the hole stabilizes and stores the energy from light via redox chemistry. The long-term objective of the Bocian/Holten&Kirmaier/Lindsey research program under this DOE grant has been to develop tetrapyrrole-based molecular architectures that absorb sunlight, funnel energy and separate charge with high efficiency. Integral to the program has been iterative cycles of design, synthesis and characterization that provided deep insights into the relationships between chemical composition, electronic structure, and key static and dynamic properties (vibrational, redox, photophysical, energy/charge transfer) of tetrapyrrolic systems. Such architectures included monomers, dyads, larger arrays, and complexes with accessory components. The objective was to develop molecular designs and guiding principles to enhance current and future energy-conversion schemes. Molecular arrays targeted to address one or more fundamental questions concerning light harvesting and energy/charge transfer were constructed from analogues of the naturally occurring hemes, chlorophylls and bacteriochlorophylls. Diverse, tunable synthetic building blocks were prepared that spanned the three respective tetrapyrrole families, which are the porphyrins, chlorins and bacteriochlorins. Thus, the research focused on porphyrins as well as synthetic surrogates for chlorophylls (chlorins, 13 1 -oxophorbines and chlorin-imides) and bacteriochlorophylls (bacteriochlorins, bacterio-13 1 -oxophorbines and bacteriochlorin-imides), generically termed hydroporphyrins. Although the three tetrapyrrole classes (porphyrins, chlorins and bacteriochlorins) absorb light strongly in the violet-blue spectral region, the long-wavelength absorption band typically lies in the green-orange, red, and NIR regions, respectively, with increasing intensity. Understanding the spectra, electronic structure, and energy/charge-transfer properties of such tetrapyrrolic macrocycles is of central importance for the rational design of molecular architectures for solar-energy conversion. Our integrated program of molecular design and synthesis coupled with a variety of spectroscopic, electrochemical, and computational studies have probed from first principles how structural and electronic properties of tetrapyrrolic macrocycles dictate spectral properties as well as the rates of ground-state hole/electron transfer and excited-state energy flow in multicomponent architectures. Individual molecules and multicomponent architectures were designed to test ideas of fundamental importance, often requiring the development of new synthetic methodology. The members of the collaborative team had almost daily discussions by phone and/or e-mail concerning design of molecules, flow of compounds between the labs, planning of physical characterization studies, discussing results and analysis and integrating into design of next generation architectures, and the preparation of manuscripts. Furthermore, students and postdocs in the different labs routinely communicated with one another to facilitate the advancement of the research activities. In short, a highly integrated and collaborative research program was well established among the groups. The research effort involved molecular design and synthesis of synthetic molecular architectures by the Lindsey group integrated with physicochemical and photophysical characterization by the Bocian group and the Holten&Kirmaier group (Figure 2). The Bocian group carried out electrochemical, electron paramagnetic resonance (EPR), resonance Raman (RR), and Fourier-transform infrared (FT-IR) studies, as well as density functional theory (DFT) calculations and the time-dependent extension (TDDFT) to gain insight into excited-state properties. The Holten&Kirmaier group carried out static and time-resolved absorption and fluorescence spectroscopy studies and simulated absorption spectra using molecular orbital (MO) energies from DFT as input to the four-orbital model to complement the TDDFT calculations. The combined measurements provided understanding of the vibrational/electronic properties of the individual molecules and the changes that occur upon incorporation into multicomponent architectures. This information underpinned elucidating the mechanisms and timescales of ground-state hole/electron transfer and excited-state energy and charge transfer.

14 SOLAR ENERGY↗
DOE OSTI2022

Low-energy moiré phonons in twisted bilayer van der Waals heterostructures

We develop a low-energy continuum model for phonons in twisted moiré bilayers, based on a configuration-space approach. In this approach, interatomic force constants are obtained from density functional theory (DFT) calculations of untwisted bilayers with various in-plane shifts. This allows for efficient computation of phonon properties for any small twist angle, while maintaining DFT-level accuracy. Based on this framework, we show how the low-energy phonon modes, including interlayer shearing and layer-breathing modes, vary with the twist angle. As the twist angle decreases, the frequencies of the low-energy modes are reordered and the atomic displacement fields corresponding to phonon eigenmodes break translational symmetry, developing periodicity on the moiré length scale. Here, we demonstrate the capabilities of our model by calculating the phonon properties of three specific structures: bilayer graphene, bilayer molybdenum disulfide (MoS 2 ), and molybdenum diselenide-tungsten diselenide (MoSe 2 -WSe 2 ).

2-dimensional systems↗
DOE OSTI2021

Increasing Iridium Oxide Activity for the Oxygen Evolution Reaction with Hafnium Modification

Synthesis and implementation of highly active, stable, and affordable electrocatalysts for the oxygen evolution reaction (OER) is a major challenge in developing energy efficient and economically viable energy conversion devices such as electrolyzers, rechargeable metal-air batteries, and regenerative fuel cells. The current benchmark electrocatalyst for OER is based on iridium oxide (IrO x ) due to its superior performance and excellent stability. However, large scale applications using IrO x are impractical due to its low abundance and high cost. In this work, we report a highly active hafnium-modified iridium oxide (IrHf x O y ) electrocatalyst for OER. The IrHf x O y electrocatalyst demonstrated ten times higher activity in alkaline conditions (pH = 11) and four times higher activity in acid conditions (pH = 1) than a IrO x electrocatalyst. The highest intrinsic mass activity of the IrHf x O y catalyst in acid conditions was calculated as 6950 A gIrO x -1 at an overpotential (η) of 0.3 V. Combined studies utilizing operando surface enhanced Raman spectroscopy (SERS) and DFT calculations revealed that the active sites for OER are the Ir-O species for both IrO x and IrHf x O y catalysts. The presence of Hf sites leads to more negative charge states on nearby O sites, and shortening the bond lengths of Ir-O, and lowering free energies for OER intermediates to accelerate the OER process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Enhancing high-fidelity neural network potentials through low-fidelity sampling

The efficacy of neural network potentials (NNPs) critically depends on the quality of the configurational datasets used for training. Prior research using empirical potentials has shown that well-selected liquid–solid transitional configurations of a metallic system can be translated to other metallic systems. This study demonstrates that such validated configurations can be relabeled using density functional theory (DFT) calculations, thereby enhancing the development of high-fidelity NNPs. Training strategies and sampling approaches are efficiently assessed using empirical potentials and subsequently relabeled via DFT in a highly parallelized fashion for high-fidelity NNP training. Our results reveal that relying solely on energy and force for NNP training is inadequate to prevent overfitting, highlighting the necessity of incorporating stress terms into the loss functions. To optimize training involving force and stress terms, we propose employing transfer learning to fine-tune the weights, ensuring that the potential surface is smooth for these quantities composed of energy derivatives. This approach markedly improves the accuracy of elastic constants derived from simulations in both empirical potential-based NNPs and relabeled DFT-based NNPs. Overall, this study offers significant insights into leveraging empirical potentials to expedite the development of reliable and robust NNPs at the DFT level.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Mechanisms and site requirements for NO and NH 3 oxidation on Cu/SSZ-13

Two series of Cu/SSZ-13 catalysts were synthesized via aqueous solution and solid-state ion exchange using SSZ-13 supports of varying Si/Al ratios. The isolated and multinuclear Cu content of these catalysts were determined by H 2 temperature programmed reduction (H 2 -TPR). Multinuclear Cu in these catalysts, including in situ Cu- dimers formed from ZCu II OH coupling and permanent CuO clusters, are active species for dry NO oxidation. NH 3 oxidation on these catalysts follows an internal SCR (i-SCR) mechanism, i.e., a portion of NH 3 is first oxidized to NO, then NO is selectively reduced by the remaining NH 3 to N 2 . NH 3 oxidation displays distinct kinetic behavior below ~300 °C and above ~400 °C. At low temperature the results indicate that NH 3 -solvated mobile Cu-ions are the active centers. CuO clusters, when present, also contribute to the low temperature activity by catalyzing NH 3 oxidation to NO. At high temperature, in situ Cu-dimers and CuO clusters catalyze NH 3 oxidation to NO, and isolated Cu-ions catalyze SCR to realize the cascade turnovers. For both NO and NH 3 oxidation, Cu-dimers balanced by framework charges of close proximity appear to be more active than Cu-dimers balanced by distant framework charges. However, the former Cu-dimers are less stable than the latter and tend to split into monomers in the presence of vicinal Brønsted acid sites. Via density functional theory (DFT) calculations, the i-SCR mechanism for low temperature NH 3 oxidation, i.e., the energetic favorability for the involvement of the NO intermediate, is justified. Furthermore, the DFT results also agree with experimental data that the formation of Cu-dimers from ZCu II OH dimerization is essential for NH 3 oxidation at high temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Thorium Monosilicide, ThSi: An Experimental and Theoretical Study

The present theoretical and experimental combination study investigates the ThSi molecule in detail. Computationally, we utilized high-level multireference and coupled-cluster levels of theory conjoined with large correlation consistent basis sets to study a series of electronic and spin–orbit states of ThSi. Here, we report potential energy curves (PECs), electron configurations at equilibrium distances, spectroscopic constants, energetics, and spin–orbit coupling effects for 16 electronic states of ThSi. The studied 16 electronic states are arranged tightly within 0.9 eV, highlighting the complexity of the electronic spectrum of ThSi. The ground electronic state of ThSi is a single-reference 1 1 Σ + state that derives from the 1σ 2 2σ 2 1π 4 electronic configuration. The Ω = 0 + spin–orbit ground state of ThSi is composed of 1 1 Σ + (47%) and 13Π (44%) electronic states. Our measured bond energy (D0) of ThSi, obtained using resonant two-photon ionization (R2PI) spectroscopy is 3.146(4) eV, where the assigned error limit is given in parentheses in units of the last quoted digits. The computed D0 of ThSi (Ω = 0 + ) at the CBS-C-CCSD(T)-δT(Q)-δDK-δSO level (3.181 eV) is in good agreement with the experimental value. Our derived enthalpy of formation for ThSi, Δ f H 0K o (ThSi(g)), is 971.8(6.0) kJ/mol. Finally, we have performed density functional theory (DFT) calculations for ThSi(1 1 Σ + ) using 16 exchange correlation functionals that span multiple rungs of “Jacob’s ladder” of density functional approximation (DFA) to assess the DFT errors on D 0 , r e , and ω e of ThSi with respect to experimental and ab initio coupled-cluster values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Solution-Phase Growth of Cu Nanowires with Aspect Ratios Greater Than 1000: Multiscale Theory

Penta-twinned metal nanowires are finding widespread application in existing and emerging technologies. However, little is known about their growth mechanisms. Here, we probe the origins of chloride- and alkylamine-mediated, solution-phase growth of penta-twinned Cu nanowires from first principles using multi-scale theory. Using quantum density functional theory (DFT) calculations, we characterize the binding and surface diffusion of Cu atoms on chlorine-covered Cu(100) and Cu(111) surfaces. We find stronger binding and slower diffusion of Cu atoms on chlorinated Cu(111) than on chlorinated Cu(100) – a reversal of the trend for bare Cu surfaces. We also probe inter-facet diffusion and find this proceeds faster from Cu(100) to Cu(111) than the reverse. Using the DFT rates for hopping between individual sites at Ångstrom scales, we calculate coarse-grained, inter-facet rates for nanowires of various lengths – up to hundreds of μm – and diameters in the 10-nm range. We predict nanowires with aspect ratios around 100 based on surface diffusion, alone. We also account for the influence of a self-assembled alkylamine layer that covers most of the {100} facets, but is absent or thin and disordered on the {111} facets and in an “end zone” near the {100}-{111} boundary. With an end zone, we predict a wide range of nanowire aspect ratios in the experimental ranges. Our work reveals the mechanisms by which a halide – chloride – promotes the growth of high-aspect-ratio nanowires.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Extending GPU-accelerated Gaussian integrals in the TeraChem software package to f type orbitals: Implementation and applications

Here, the increasing availability of graphics processing units (GPUs) for scientific computing has prompted interest in accelerating quantum chemical calculations through their use. However, the complexity of integral kernels for high angular momentum basis functions often limits the utility of GPU implementations with large basis sets or for metal containing systems. In this work, we report the implementation of f function support in the GPU-accelerated TeraChem software package through the development of efficient kernels for the evaluation of Hamiltonian integrals. The high efficiency of the resulting code is demonstrated through density functional theory (DFT) calculations on increasingly large organic molecules and transition metal complexes, as well as coupled cluster singles and doubles calculations on water clusters. Preliminary investigations into Ni(I) catalysis with DFT and the photochemistry of MnH(CH 3 ) with complete active space self-consistent field are also carried out. Overall, our GPU-accelerated software appears to be well-suited for fast simulation of large transition metal containing systems, as well as organic molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Understanding the influence of solvents on the Pt-catalyzed hydrodeoxygenation of guaiacol

Bio-oil derived from fast pyrolysis of lignocellulosic biomass needs to be deoxygenated to become a substitute for petroleum fuels. Here, we study the hydrodeoxygenation mechanism of guaiacol, a bio-oil model compound derived from the lignin fraction of biomass, on Pt(1 1 1) terrace sites in the presence of water, diethyl ether, 1-butanol, and n-hexane as solvent. Using first-principles periodic density functional theory (DFT) calculations and mean-field microkinetic reactor modeling, a detailed reaction mechanism is investigated targeting various products such as catechol, phenol, anisole, benzene, cyclohexanone, and cyclohexanol. Solvent phase DFT outcomes are mostly similar to that of the vapor phase; however, microkinetic modeling results suggest that rate controlling species and transition states differ somewhat in the various reaction environments. Catechol was found to be the major aromatic product across all reaction environments. Over Pt(1 1 1), unsaturated monooxygenate production from catechol is unlikely.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Trends in Siting of Metals in Heterometallic Nd–Yb Metal–Organic Frameworks and Molecular Crystals

Several studies suggest that metal ordering within metal–organic frameworks (MOFs) is important for understanding how MOFs behave in relevant applications; however, these siting trends can be difficult to determine experimentally. To garner insight into the energetic driving forces that may lead to nonrandom ordering within heterometallic MOFs, we employ density functional theory (DFT) calculations on several bimetallic metal–organic crystals composed of Nd and Yb metal atoms. We also investigate the metal siting trends for a newly synthesized MOF. Our DFT-based energy of mixing results suggest that Nd will likely occupy sites with greater access to electronegative atoms and that local homometallic domains within a mixed-metal Nd–Yb system are favored. We also explore the use of less computationally extensive methods such as classical force fields and cluster expansion models to understand their feasibility for large system sizes. Here, this study highlights the impact of metal ordering on the energetic stability of heterometallic MOFs and crystal structures.

energy of mixing↗
DOE OSTI2022

Sulfated Zirconium Metal–Organic Frameworks as Well-Defined Supports for Enhancing Organometallic Catalysis

Understanding heterogeneous catalysts is a challenging pursuit due to surface site nonuniformity and aperiodicity in traditionally-used materials. One example is sulfated metal oxides, which function as highly-active catalysts and as supports for organometallic complexes. These applications are due to traits such as acidity, ability to act as a weakly coordinating ligand, and aptitude for promoting transformations via radical cation intermediates. Research is ongoing about the structural features of sulfated metal oxides that imbue the aforementioned properties, such as sulfate geometry and coordination. To better understand these materials, metal–organic frameworks (MOFs) have been targeted as structurally-defined analogs. Composed of inorganic nodes and organic linkers, MOFs possess features such as high porosity and crystallinity, which make them attractive for mechanistic studies of heterogeneous catalysts. In this work, a Zr 6 -based MOF NU-1000 is sulfated and characterized using atomically-precise techniques such as single crystal X-ray diffraction (SCXRD) in addition to diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS). Here, the dynamic nature of the sulfate binding motif is found to transition from monodentate, to bidentate, to tridentate depending on the degree of hydration, as supported by density functional theory (DFT) calculations. Heightened Brønsted acidity compared to the parent MOF was observed upon sulfation, and probed through trimethylphosphine oxide (TMPO) physisorption, ammonia sorption, in-situ ammonia DRIFTS, and DFT studies. With the support structure benchmarked, an organoiridium complex was chemisorbed onto the sulfated MOF node and the efficacy of this supported catalyst was demonstrated for stoichiometric and catalytic activation of benzene-d6 and toluene with structure-activity relationships derived.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Comparing Classical and Machine Learning Force Fields for Modeling Deformation of Metal–Organic Frameworks Relevant for Direct Air Capture

Deformation of metal–organic frameworks (MOFs) induced by adsorbate molecules can affect adsorption properties such as capacity and selectivity, but most computational studies of MOFs assume framework rigidity to simplify calculations. Although flexible force fields (FFs) for MOFs have been parametrized for specific materials, the generality of FFs for reliably modeling adsorbate-induced deformation to accuracy nearing that of density functional theory (DFT) has not been established. This work confirms using DFT calculations that adsorbate-induced deformation can affect CO 2 and H 2 O adsorption energies in a considerable fraction of MOFs promising for direct air capture (DAC). We then benchmark the efficacy of several general-purpose FFs in describing adsorbate-induced deformation for DAC against DFT. Our results show that current classical FFs are insufficient for describing MOF deformation, especially in cases of interest for DAC where strong interactions exist between adsorbed molecules and MOF frameworks. Some emerging machine learning force fields (MLFFs) we tested, particularly CHGNet, MACE-MP-0, and Equiformer V2, appear to be more promising than the classical FF for emulating the deformation behavior described by DFT. The best performing FF (CHGNet), however, fails to achieve the accuracy required for practical predictions with a mean absolute adsorption energy error of 0.124 eV.

adsorption↗
DOE OSTI2022

Electron Counting and High-Pressure Phase Transformations in Metal Hexaborides

Pressure-induced structural transitions of the alkaline earth hexaborides, CaB 6 , SrB 6 , and BaB 6 , are studied theoretically using electron counting rules and density functional theory calculations. Here, we demonstrate the applicability of gas-phase borane electron counting methods to solid-state metal borides under pressure and validate the assumptions of the rules by density functional theory (DFT) calculations. All three compounds share ambient-pressure and high-pressure structures, but BaB 6 differs from CaB 6 and SrB 6 at intermediate pressures. The unique BaB 6 phase is shown to break electron counting rules, while all other phases obey them. This anomaly is resolved by DFT, which reveals B–Ba covalency and unusual B–B π bonding under pressure. The relationships between structure and bonding can help us to understand the exotic behavior of lanthanide hexaborides and design new borides with desirable properties. Developing electron counting procedures for solids will enhance materials discovery efforts with chemical intuition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Bulk Fermi surfaces of the Dirac type-II semimetallic candidate NiTe 2

Here, we present a study on the Fermi surface of the Dirac type-II semimetallic candidate NiTe 2 via the temperature and angular dependence of the de Haas–van Alphen (dHvA) effect measured in single crystals grown through Te flux. In contrast to its isostructural compounds like PtSe 2 , band-structure calculations predict NiTe 2 to display a tilted Dirac node very close to its Fermi level that is located along the Γ to A high-symmetry direction within its first Brillouin zone. The angular dependence of the dHvA frequencies is found to be in agreement with the first-principles calculations when the electronic bands are slightly shifted with respect to the Fermi level (ε F ), and therefore provide support for the existence of a Dirac type-II node in NiTe 2 . Nevertheless, we observed mild disagreements between experimental observations and density functional theory (DFT) calculations as, for example, nearly isotropic and light experimental effective masses. This indicates that the dispersion of the bands is not well captured by DFT. Despite the coexistence of Dirac-like fermions with topologically trivial carriers, samples of the highest quality display an anomalous and large either linear or sublinear magnetoresistivity. Finally, this suggests that Lorentz invariance breaking Dirac-like quasiparticles dominate the carrier transport in this compound.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Hydride Accessibility and Reactivity in the Configurational and Stoichiometric Space of β-Ga 2 O 3 for CO 2 Hydrogenation

Understanding how surface species evolve under reaction conditions is essential for improving catalyst design for efficient CO 2 hydrogenation. This work combines systematic DFT calculations with grand canonical sampling to investigate the stability and reactivity of Ga–H species on β-Ga 2 O 3 across a range of reaction conditions. Initial DFT studies reveal that when Ga–H species are present, they facilitate formate formation via a low-barrier pathway, largely independent of the surface termination or hydrogen site. However, grand canonical sampling shows that under a broad range of reaction conditions─especially at high oxygen chemical potentials associated with high water content─Ga–H species are thermodynamically inaccessible. Furthermore, adsorbed water molecules can block reactive sites, inhibiting CO 2 activation even when hydrides are present. These findings suggest that the lack of accessible hydride species, rather than their intrinsic reactivity, could contribute to reduced catalytic performance of β-Ga 2 O 3 under more oxidizing, high-conversion conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Kohn–Sham accuracy from orbital-free density functional theory via Δ-machine learning

Here, we present a Δ-machine learning model for obtaining Kohn–Sham accuracy from orbital-free density functional theory (DFT) calculations. In particular, we employ a machine-learned force field (MLFF) scheme based on the kernel method to capture the difference between Kohn–Sham and orbital-free DFT energies/forces. We implement this model in the context of on-the-fly molecular dynamics simulations and study its accuracy, performance, and sensitivity to parameters for representative systems. We find that the formalism not only improves the accuracy of Thomas–Fermi–von Weizsäcker orbital-free energies and forces by more than two orders of magnitude but is also more accurate than MLFFs based solely on Kohn–Sham DFT while being more efficient and less sensitive to model parameters. We apply the framework to study the structure of molten Al 0.88 Si 0.12 , the results suggesting no aggregation of Si atoms, in agreement with a previous Kohn–Sham study performed at an order of magnitude smaller length and time scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗