BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Giulia Galli

Giulia Galli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-sciquant-phphysics.chem-phphysics.comp-phcond-mat.mes-hallcond-mat.softcond-mat.str-elcond-matcond-mat.stat-mechEmerging Technologiesmtrl-thNeural and Evolutionary Computingphysics.app-ph

Catalog footprint

What is connected

28works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2025arXiv

Probing aqueous interfaces with spin defects

Understanding the physical and chemical properties of aqueous interfaces is important in diverse fields of science, ranging from biology and chemistry to materials science. In spite of crucial progress in surface sensitive spectroscopic techniques over the past decades, the microscopic properties of aqueous interfaces remain difficult to measure. Here we explore the use of noise spectroscopy to characterize interfacial properties, specifically of quantum sensors hosted in two-dimensional materials in contact with water. We combine molecular dynamics simulations of water/graphene interfaces and the calculations of the spin dynamics of an NV-like color center, and we investigate the impact of interfacial water and simple ions on the decoherence time of the defect. We show that the Hahn echo coherence time of the NV center is sensitive to motional narrowing and to the hydrogen bonding arrangement and the dynamical properties of water and ions at the interface. We present results as a function of the liquid temperature, strength of the water-surface interaction, and for varied mono-valent and di-valent ions, highlighting the broad applicability of near-surface qubits to gain insight into the properties of aqueous interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Green's function formulation of quantum defect embedding theory

We present a Green's function formulation of the quantum defect embedding theory (QDET) where a double counting scheme is rigorously derived within the $G_0 W_0$ approximation. We then show the robustness of our methodology by applying the theory with the newly derived scheme to several defects in diamond. Additionally, we discuss a strategy to obtain converged results as a function of the size and composition of the active space. Our results show that QDET is a promising approach to investigate strongly correlated states of defects in solids.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum Embedding Theories to Simulate Condensed Systems on Quantum Computers

Quantum computers hold promise to improve the efficiency of quantum simulations of materials and to enable the investigation of systems and properties more complex than tractable at present on classical architectures. Here, we discuss computational frameworks to carry out electronic structure calculations of solids on noisy intermediate scale quantum computers using embedding theories, and we give examples for a specific class of materials, i.e., spin defects in solids. These are promising systems to build future quantum technologies, e.g., computers, sensors and devices for quantum communications. Although quantum simulations on quantum architectures are in their infancy, promising results for realistic systems appear to be within reach.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum materials for energy-efficient neuromorphic computing

Neuromorphic computing approaches become increasingly important as we address future needs for efficiently processing massive amounts of data. The unique attributes of quantum materials can help address these needs by enabling new energy-efficient device concepts that implement neuromorphic ideas at the hardware level. In particular, strong correlations give rise to highly non-linear responses, such as conductive phase transitions that can be harnessed for short and long-term plasticity. Similarly, magnetization dynamics are strongly non-linear and can be utilized for data classification. This paper discusses select examples of these approaches, and provides a perspective for the current opportunities and challenges for assembling quantum-material-based devices for neuromorphic functionalities into larger emergent complex network systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Vibrationally resolved optical excitations of the nitrogen-vacancy center in diamond

A comprehensive description of the optical cycle of spin defects in solids requires the understanding of the electronic and atomistic structure of states with different spin multiplicity, including singlet states which are particularly challenging from a theoretical standpoint. We present a general framework, based on spin-flip time-dependent density function theory, to determine the excited state potential energy surfaces of the many-body singlet states of spin defects; we then predict the vibrationally resolved absorption spectrum between singlet shelving states of a prototypical defect, the nitrogen-vacancy center in diamond. Our results, which are in excellent agreement with experiments, provide an interpretation of the measured spectra and reveal the key role of specific phonons in determining absorption processes, and the notable influence of non-adiabatic interactions. The insights gained from our calculations may be useful in defining strategies to improve infrared-absorption-based magnetometry and optical pumping schemes. The theoretical framework developed here is general and applicable to a variety of other spin defects and materials.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Generalized scaling of spin qubit coherence in over 12,000 host materials

Spin defect centers with long quantum coherence times ($T_2$) are key solid-state platforms for a variety of quantum applications. Recently, cluster correlation expansion (CCE) techniques have emerged as a powerful tool to simulate the $T_2$ of defect electron spins in these solid-state systems with good accuracy. Here, based on CCE, we uncover an algebraic expression for $T_2$ generalized for host compounds with dilute nuclear spin baths, which enables a quantitative and comprehensive materials exploration with a near instantaneous estimate of the coherence. We investigate more than 12,000 host compounds at natural isotopic abundance, and find that silicon carbide (SiC), a prominent widegap semiconductor for quantum applications, possesses the longest coherence times among widegap non-chalcogenides. In addition, more than 700 chalcogenides are shown to possess a longer $T_2$ than SiC. We suggest new potential host compounds with promisingly long $T_2$ up to 47 ms, and pave the way to explore unprecedented functional materials for quantum applications.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Machine Learning Dielectric Screening for the Simulation of Excited State Properties of Molecules and Materials

Accurate and efficient calculations of absorption spectra of molecules and materials are essential for the understanding and rational design of broad classes of systems. Solving the Bethe-Salpeter equation (BSE) for electron-hole pairs usually yields accurate predictions of absorption spectra, but it is computationally expensive, especially if thermal averages of spectra computed for multiple configurations are required. We present a method based on machine learning to evaluate a key quantity entering the definition of absorption spectra: the dielectric screening. We show that our approach yields a model for the screening that is transferable between multiple configurations sampled during first principles molecular dynamics simulations; hence it leads to a substantial improvement in the efficiency of calculations of finite temperature spectra. We obtained computational gains of one to two orders of magnitude for systems with 50 to 500 atoms, including liquids, solids, nanostructures, and solid/liquid interfaces. Importantly, the models of dielectric screening derived here may be used not only in the solution of the BSE but also in developing functionals for time-dependent density functional theory (TDDFT) calculations of homogeneous and heterogeneous systems. Overall, our work provides a strategy to combine machine learning with electronic structure calculations to accelerate first principles simulations of excited-state properties.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Spin-spin interactions in solids from mixed all-electron and pseudopotential calculations $-$ a path to screening materials for spin qubits

Understanding the quantum dynamics of spin defects and their coherence properties requires accurate modeling of spin-spin interaction in solids and molecules, for example by using spin Hamiltonians with parameters obtained from first-principles calculations. We present a real-space approach based on density functional theory for the calculation of spin-Hamiltonian parameters, where only selected atoms are treated at the all-electron level, while the rest of the system is described with the pseudopotential approximation. Our approach permits calculations for systems containing more than 1000 atoms, as demonstrated for defects in diamond and silicon carbide. We show that only a small number of atoms surrounding the defect needs to be treated at the all-electron level, in order to obtain an overall all-electron accuracy for hyperfine and zero-field splitting tensors. We also present results for coherence times, computed with the cluster correlation expansion method, highlighting the importance of accurate spin-Hamiltonian parameters for quantitative predictions of spin dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Stability and molecular pathways to the formation of spin defects in silicon carbide

Spin defects in wide-bandgap semiconductors provide a promising platform to create qubits for quantum technologies. Their synthesis, however, presents considerable challenges, and the mechanisms responsible for their generation or annihilation are poorly understood. Here, we elucidate spin defect formation processes in a binary crystal for a key qubit candidate--the divacancy complex (VV) in silicon carbide (SiC). Using atomistic models, enhanced sampling simulations, and density functional theory calculations, we find that VV formation is a thermally activated process that competes with the conversion of silicon ($V_{Si}$) to carbon monovacancies ($V_{C}$), and that VV reorientation can occur without dissociation. We also find that increasing the concentration of $V_{Si}$ relative to $V_{C}$ favors the formation of divacancies. Moreover, we identify pathways to create spin defects consisting of antisite-double vacancy complexes and determine their electronic properties. The detailed view of the mechanisms that underpin the formation and dynamics of spin defects presented here may facilitate the realization of qubits in an industrially relevant material.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Probing the coherence of solid-state qubits at avoided crossings

Optically addressable paramagnetic defects in wide-band-gap semiconductors are promising platforms for quantum communications and sensing. The presence of avoided crossings between the electronic levels of these defects can substantially alter their quantum dynamics and be both detrimental and beneficial for quantum information applications. Avoided crossings give rise to clock transitions, which can significantly improve protection from magnetic noise and favorably increase coherence time. However, the reduced coupling between electronic and nuclear spins at an avoided crossing may be detrimental to applications where nuclear spins act as quantum memories. Here we present a combined theoretical and experimental study of the quantum dynamics of paramagnetic defects interacting with a nuclear spin bath at avoided crossings. We develop a computational approach based on a generalization of the cluster expansion technique, which can account for processes beyond pure dephasing and describe the dynamics of any solid-state spin-qubits near avoided crossings. Using this approach and experimental validation, we determine the change in nature and source of noise at avoided crossings for divacancies in SiC. We find that we can condition the clock transition of the divacancies in SiC on multiple adjacent nuclear spins states. In our experiments, we demonstrate that one can suppress the effects of fluctuating charge impurities with depletion techniques, leading to an increased coherence time at clock transition, limited purely by magnetic noise. Combined with ab-initio predictions of spin Hamiltonian parameters, the proposed theoretical approach paves the way to designing the coherence properties of spin qubits from first principles.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

PyCDFT: A Python package for constrained density functional theory

We present PyCDFT, a Python package to compute diabatic states using constrained density functional theory (CDFT). PyCDFT provides an object-oriented, customizable implementation of CDFT, and allows for both single-point self-consistent-field calculations and geometry optimizations. PyCDFT is designed to interface with existing density functional theory (DFT) codes to perform CDFT calculations where constraint potentials are added to the Kohn-Sham Hamiltonian. Here we demonstrate the use of PyCDFT by performing calculations with a massively parallel first-principles molecular dynamics code, Qbox, and we benchmark its accuracy by computing the electronic coupling between diabatic states for a set of organic molecules. We show that PyCDFT yields results in agreement with existing implementations and is a robust and flexible package for performing CDFT calculations. The program is available at https://github.com/hema-ted/pycdft/.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantum simulations of materials on near-term quantum computers

Quantum computers hold promise to enable efficient simulations of the properties of molecules and materials; however, at present they only permit ab initio calculations of a few atoms, due to a limited number of qubits. In order to harness the power of near-term quantum computers for simulations of larger systems, it is desirable to develop hybrid quantum-classical methods where the quantum computation is restricted to a small portion of the system. This is of particular relevance for molecules and solids where an active region requires a higher level of theoretical accuracy than its environment. Here we present a quantum embedding theory for the calculation of strongly-correlated electronic states of active regions, with the rest of the system described within density functional theory. We demonstrate the accuracy and effectiveness of the approach by investigating several defect quantum bits in semiconductors that are of great interest for quantum information technologies. We perform calculations on quantum computers and show that they yield results in agreement with those obtained with exact diagonalization on classical architectures, paving the way to simulations of realistic materials on near-term quantum computers.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Improving the efficiency of $G_0W_0$ calculations with approximate spectral decompositions of dielectric matrices

Recently it was shown that the calculation of quasiparticle energies using the $G_0W_0$ approximation can be performed without computing explicitly any virtual electronic states, by expanding the Green function and screened Coulomb interaction in terms of the eigenstates of the static dielectric matrix. Avoiding the evaluation of virtual electronic states leads to improved efficiency and ease of convergence of $G_0W_0$ calculations. Here we propose a further improvement of the efficiency of these calculations, based on an approximation of density-density response functions of molecules and solids. The approximation relies on the calculation of a subset of eigenvectors of the dielectric matrix using the kinetic operator instead of the full Hamiltonian, and it does not lead to any substantial loss of accuracy for the quasiparticle energies. The computational savings introduced by this approximation depend on the system, and they become more substantial as the number of electrons increases.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

MatD3: A Database and Online Presentation Package for Research Data Supporting Materials Discovery, Design, and Dissemination

The discovery of new materials as well as the determination of a vast set of materials properties for science and technology is a fast growing field of research, with contributions from many groups worldwide. Materials data from individual research groups is traditionally disseminated by means of loosely interconnected, peer-reviewed publications. MatD3 is an open-source, dedicated database and web application framework designed to store, curate and disseminate experimental and theoretical materials data generated by individual research groups or research consortia. A research group can set up its own instance of MatD3 and publish scientific results or simply use an existing online MatD3 instance. Disseminating research data in this form enables broader access, reproducibility, and repurposing of scientific products. MatD3 is a general purpose database that does not focus on any specific level of theory or experimental method. Instead, the focus is on storing and making accessible the data and making it straightforward to curate them.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Design of defect spins in piezoelectric aluminum nitride for solid-state hybrid quantum technologies

Spin defects in wide-band gap semiconductors are promising systems for the realization of quantum bits, or qubits, in solid-state environments. To date, defect qubits have only been realized in materials with strong covalent bonds. Here, we introduce a strain-driven scheme to rationally design defect spins in functional ionic crystals, which may operate as potential qubits. In particular, using a combination of state-of-the-art ab-initio calculations based on hybrid density functional and many-body perturbation theory, we predicted that the negatively charged nitrogen vacancy center in piezoelectric aluminum nitride exhibits spin-triplet ground states under realistic uni- and bi-axial strain conditions; such states may be harnessed for the realization of qubits. The strain-driven strategy adopted here can be readily extended to a wide range of point defects in other wide-band gap semiconductors, paving the way to controlling the spin properties of defects in ionic systems for potential spintronic technologies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonempirical Range-separated Hybrid Functionals for Solids and Molecules

Dielectric-dependent hybrid (DDH) functionals were recently shown to yield accurate energy gaps and dielectric constants for a wide variety of solids, at a computational cost considerably less than that of GW calculations. The fraction of exact exchange included in the definition of DDH functionals depends (self-consistently) on the dielectric constant of the material. Here we introduce a range-separated (RS) version of DDH functionals where short and long-range components are matched using system dependent, non-empirical parameters. We show that RS DDHs yield accurate electronic properties of inorganic and organic solids, including energy gaps and absolute ionization potentials. Furthermore we show that these functionals may be generalized to finite systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quantum decoherence dynamics of divacancy spins in silicon carbide

Long coherence times are key to the performance of quantum bits (qubits). Here, we experimentally and theoretically show that the Hahn-echo coherence time (T2) of electron spins associated with divacancy defects in 4H-SiC reaches 1.3 ms, one of the longest T2 times of an electron spin in a naturally isotopic crystal. Using a first-principles microscopic quantum-bath model, we find that two factors determine the unusually robust coherence. First, in the presence of moderate magnetic fields (300 G and above), the 29Si and 13C paramagnetic nuclear spin baths are decoupled. In addition, because SiC is a binary crystal, homo-nuclear spin pairs are both diluted and forbidden from forming strongly coupled, nearest-neighbor spin pairs. Longer neighbor distances result in fewer nuclear spin flip-flops, a less fluctuating intra-crystalline magnetic environment, and thus a longer T2 time. Our results point to polyatomic crystals as promising hosts for coherent qubits in the solid state.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Large scale GW calculations

We present GW calculations of molecules, ordered and disordered solids and interfaces, which employ an efficient contour deformation technique for frequency integration, and do not require the explicit evaluation of virtual electronic states, nor the inversion of dielectric matrices. We also present a parallel implementation of the algorithm which takes advantage of separable expressions of both the single particle Green's function and the screened Coulomb interaction. The method can be used starting from density functional theory calculations performed with semi-local or hybrid functionals. We applied the newly developed technique to GW calculations of systems of unprecedented size, including water/semiconductor interfaces with thousands of electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Self-consistent hybrid functional for condensed systems

A self-consistent scheme for determining the optimal fraction of exact exchange for full-range hybrid functionals is presented and applied to the calculation of band gaps and dielectric constants of solids. The exchange-correlation functional is defined in a similar manner to the PBE0 functional, but the mixing parameter is set equal to the inverse macroscopic dielectric function and it is determined self-consistently by computing the optimal dielectric screening. We found excellent agreement with experiments for the properties of a broad class of systems, with band gaps ranging between 0.7 and 21.7 eV and dielectric constants within 1.23 and 15.9. We propose that the eigenvalues and eigenfunctions obtained with the present self-consistent hybrid scheme may be excellent inputs for G$_0$W$_0$ calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tetrahedrally coordinated carbonates in Earth's lower mantle

Carbonates are the main species that bring carbon deep into our planet through subduction. They are an important rock-forming mineral group, fundamentally distinct from silicates in Earth's crust in that carbon binds to three oxygen atoms, while silicon is bonded to four oxygens. Here, we present experimental evidence that under the sufficiently high pressures and high temperatures existing in the lower mantle, ferromagnesian carbonates transform to a phase with tetrahedrally coordinated carbons. Above 80 GPa, in situ synchrotron infrared experiments show the unequivocal spectroscopic signature of the high-pressure phase of (Mg,Fe)CO$_3$. Using ab-initio calculations, we assign the new IR signature to C-O bands associated with tetrahedrally coordinated carbon with asymmetric C-O bonds. Tetrahedrally coordinated carbonates are expected to exhibit substantially different reactivity than low pressure three-fold coordinated carbonates, as well as different chemical properties in the liquid state. Hence this may have significant implications on carbon reservoirs and fluxes and the global geodynamic carbon cycle.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The refractive index and electronic gap of water and ice increase with increasing pressure

Determining the electronic and dielectric properties of water at high pressure and temperature is an essential prerequisite to understand the physical and chemical properties of aqueous environments under supercritical conditions, e.g. in the Earth interior. However optical measurements of compressed ice and water remain challenging and it has been common practice to assume that their band gap is inversely correlated to the measured refractive index, consistent with observations reported for hundreds of materials. Here we report ab initio molecular dynamics and electronic structure calculations showing that both the refractive index and the electronic gap of water and ice increase with pressure, at least up to 30 GPa. Subtle electronic effects, related to the nature of interband transitions and band edge localization under pressure, are responsible for this apparently anomalous behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ab initio investigation of the melting line of nitrogen at high pressure

Understanding the behavior of molecular systems under pressure is a fundamental problem in condensed matter physics. In the case of nitrogen, the determination of the phase diagram and in particular of the melting line, are largely open problems. Two independent experiments have reported the presence of a maximum in the nitrogen melting curve, below 90 GPa, however the position and the interpretation of the origin of such maximum differ. By means of ab initio molecular dynamics simulations based on density functional theory and thermodynamic integration techniques, we have determined the phase diagram of nitrogen in the range between 20 and 100 GPa. We find a maximum in the melting line, related to a transformation in the liquid, from molecular N_2 to polymeric nitrogen accompanied by an insulator-to-metal transition.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

First-Principles Study of Electronic and Vibrational Properties of BaHfN$_2$

The transition metal nitride BaHfN$_2$, which consists of weakly bonded neutral slabs of closed shell ions, has structural and chemical similarities to other layered nitrides which have impressive superconducting T$_c$ when electron doped: A$_x$HfNCl, A$_x$ZrNCl, A$_x$TiNCl, with $T_c= 25.5$, $15.2$ and $16.5$ K, respectively for appropriate donor (A) concentrations $x$. These similarities suggest the possibility of BaHfN$_2$ being another relatively high T$_c$ nitride upon doping, with effects of structure and the role of specific transition metal ions yet to be understood. We report first-principles electronic structure calculations for stoichiometric BaHfN$_2$ using density functional theory with plane-wave basis sets and separable dual-space Gaussian pseudopotentials. An indirect band gap of 0.8 eV was obtained and the lowest conduction band is primarily of Hf 5$d_{xy}$ character, similar to $β$-ZrNCl and $α$-TiNCl. The two N sites, one in the Hf layer and another one in the Ba layer, were found to have very anisotropic Born effective charges (BEC):deviations from the formal charge (-3) are opposite for the two sites, and opposite for the two orientations (in-plane, out of plane). LO-TO splittings and comparison of BECs and dielectric constant tensors to those of related compounds are discussed, and the effect of electron doping on the zone-center phonons is reported.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The nature and strength of inter-layer binding in graphite

We computed the inter-layer bonding properties of graphite using an ab-initio many body theory. We carried out variational and diffusion quantum Monte Carlo calculations and found an equilibrium inter-layer binding energy in good agreement with most recent experiments. We also analyzed the behavior of the total energy as a function of interlayer separation at large distances comparing the results with the predictions of the random phase approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Theoretical investigation of methane under pressure

We present computer simulations of liquid and solid phases of condensed methane at pressures below 25 GPa, between 150 and 300 K, where no appreciable molecular dissociation occurs. We used molecular dynamics (MD) and metadynamics techniques, and empirical potentials in the rigid molecule approximation, whose validity was confirmed a posteriori by carrying out it ab initio MD simulations for selected pressure and temperature conditions. Our results for the melting line are in satisfactory agreement with existing measurements. We find that the fcc crystal transforms into a hcp structure with 4 molecules per unit cell (B phase) at about 10 GPa and 150 K, and that the B phase transforms into a monoclinic high pressure phase above 20 GPa. Our results for solid/solid phase transitions are consistent with those of Raman studies but the phase boundaries estimated in our calculations are at higher pressure than those inferred from spectroscopic data.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

The electronic structure of liquid water within density functional theory

In the last decade, computational studies of liquid water have mostly concentrated on ground state properties. However recent spectroscopic measurements have been used to infer the structure of water, and the interpretation of optical and x-ray spectra requires accurate theoretical models of excited electronic states, not only of the ground state. To this end, we investigate the electronic properties of water at ambient conditions using ab initio density functional theory within the generalized gradient approximation (DFT/GGA), focussing on the unoccupied subspace of Kohn-Sham eigenstates. We generate long (250 ps) classical trajectories for large supercells, up to 256 molecules, from which uncorrelated configurations of water molecules are extracted for use in DFT/GGA calculations of the electronic structure. We find that the density of occupied states of this molecular liquid is well described with 32 molecule supercells using a single k-point (k = 0) to approximate integration over the first Brillouin zone. However, the description of the density of unoccupied states (u-EDOS) is sensitive to finite size effects. Small, 32 molecule supercell calculations, using Gamma-the point approximation, yield a spuriously isolated state above the Fermi level. Nevertheless, the more accurate u-EDOS of large, 256 molecule supercells may be reproduced using smaller supercells and increased k-point sampling. This indicates that the electronic structure of molecular liquids like water is relatively insensitive to the long-range disorder in the molecular structure. These results have important implications for efficiently increasing the accuracy of spectral calculations for water and other molecular liquids.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Large scale quantum simulations: C_60 impacts on a semiconducting surface

We present tight binding molecular dynamics simulations of C_60 collisions on the reconstructed diamond(111) surface, carried out with an O(N) method and with cells containing 1140 atoms. The results of our simulations are in very good agreement with experiments performed under the same impact conditions. Furthermore our calculations provide a detailed characterization of the microscopic processes occuring during the collision, and allow the identification of three impact regimes, as a function of the fullerene incident energy. Finally, the study of the reactivity between the cluster and the surface gives insight into the deposition mechanisms of C_60 on semiconducting substrates.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Total energy global optimizations using non orthogonal localized orbitals

An energy functional for orbital based $O(N)$ calculations is proposed, which depends on a number of non orthogonal, localized orbitals larger than the number of occupied states in the system, and on a parameter, the electronic chemical potential, determining the number of electrons. We show that the minimization of the functional with respect to overlapping localized orbitals can be performed so as to attain directly the ground state energy, without being trapped at local minima. The present approach overcomes the multiple minima problem present within the original formulation of orbital based $O(N)$ methods; it therefore makes it possible to perform $O(N)$ calculations for an arbitrary system, without including any information about the system bonding properties in the construction of the input wavefunctions. Furthermore, while retaining the same computational cost as the original approach, our formulation allows one to improve the variational estimate of the ground state energy, and the energy conservation during a molecular dynamics run. Several numerical examples for surfaces, bulk systems and clusters are presented and discussed.

0 citations0 reviewsNo code linkedOpen access