BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Ali Alavi

Ali Alavi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.comp-phcond-mat.mtrl-scicond-mat.str-elphysics.chem-phcond-mat.supr-conquant-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2022arXiv

Ab initio wavefunction analysis of electron removal quasi-particle state of NdNiO$_2$ with fully correlated quantum chemical methods

The discovery of superconductivity in hole-doped infinite-layer NdNiO$_2$ -- a transition metal (TM) oxide that is both isostructural and isoelectronic to cuprate superconductors -- has lead to renewed enthusiasm in the hope of understanding the origin of unconventional superconductivity. Here, we investigate the electron-removal states in infinite-layered Ni$^{1+}$ oxide, NdNiO$_2$, which mimics hole-doping, with the state-of-the-art many-body multireference quantum chemistry methods. From the analysis of the many-body wavefunction, we find that the hole-doped $d^8$ ground state of NdNiO$_2$ is very different from the $d^8$ ground state in isostructural cuprate analog CaCuO$_2$, although the parent $d^9$ ground states are for the most part identical. We show that the doped hole in NdNiO$_2$ mainly localizes on the Ni $3d_{x^2-y^2}$ orbital to form a closed-shell singlet, and this singlet configuration contributes to $\sim$40% of the wavefunction. In contrast, in CaCuO$_2$ the Zhang-Rice singlet configurations contribute to $\sim$65% of the wavefunction. With the help of the quantum information concept of entanglement entropy, we quantify the different types of electronic correlations in the nickelate and cuprate compounds and find that the dynamic radial-type correlations within the Ni $d$ manifold are persistent in hole-doped NdNiO$_2$. As a result, the $d^8$ multiplet effects are stronger and the additional hole foot-print is more three-dimensional in NdNiO$_2$. Our analysis shows that the most commonly used three-band Hubbard model employed to express the doped scenario in cuprates represents $\sim$90% of the $d^8$ wavefunction for CaCuO$_2$, but such a model grossly approximates the $d^8$ wavefunction for NdNiO$_2$ as it only stands for $\sim$60% of the wavefunction.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

General analytical nuclear force and molecular potential energy surface from full configuration interaction quantum Monte Carlo

Full configuration interaction quantum Monte Carlo (FCIQMC) is a state-of-the-art stochastic electronic structure method, providing a methodology to compute FCI-level state energies of molecular systems within a quantum chemical basis. However, especially to probe {\em dynamics} at the FCIQMC level, it is necessary to devise more efficient schemes to produce nuclear forces and potential energy surfaces (PES) from FCIQMC. In this work, we derive the general formula for nuclear force from FCIQMC, and clarify different contributions of the total force. This method to obtain FCIQMC forces eliminates previous restrictions, and can be used with frozen core approximation and free selection of orbitals, making it promising for more efficient nuclear force calculations. After numerical check of this procedure on the binding curve of N$_2$ molecule, we use the FCIQMC energy and force to obtain the full-dimensional ground state PES of water molecule via Gaussian processes regression. The new water FCIQMC PES can be used as the basis for H$_2$O ground state nuclear dynamics, structure optimization, and rotation-vibrational spectrum calculation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

General embedded cluster protocol for accurate modeling of oxygen vacancies in metal-oxides

The O vacancy (Ov) formation energy, $E_\textrm{Ov}$, is an important property of a metal-oxide, governing its performance in applications such as fuel cells or heterogeneous catalysis. These defects are routinely studied with density functional theory (DFT). However, it is well-recognized that standard DFT formulations (e.g. the generalized gradient approximation) are insufficient for modeling the Ov, requiring higher levels of theory. The embedded cluster method offers a promising approach to compute $E_\textrm{Ov}$ accurately, giving access to all electronic structure methods. Central to this approach is the construction of quantum(-mechanically treated) clusters placed within suitable embedding environments. Unfortunately, current approaches to constructing the quantum clusters either require large system sizes, preventing application of high-level methods, or require significant manual input, preventing investigations of multiple systems simultaneously. In this work, we present a systematic and general quantum cluster design protocol that can determine small converged quantum clusters for studying the Ov in metal-oxides with accurate methods such as local coupled cluster with single, double and perturbative triple excitations [CCSD(T)]. We apply this protocol to study the Ov in the bulk and surface planes of rutile TiO2 and rocksalt MgO, producing the first accurate and well-converged determinations of $E_\textrm{Ov}$ with this method. These reference values are used to benchmark exchange-correlation functionals in DFT and we find that all studied functionals underestimate $E_\textrm{Ov}$, with the average error decreasing along the rungs of Jacob's ladder. This protocol is automatable for high-throughput calculations and can be generalized to study other point defects or adsorbates.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Combined unitary and symmetric group approach applied to low-dimensional spin systems

A novel combined unitary and symmetric group approach is used to study the spin-$\frac{1}{2}$ Heisenberg model and related Fermionic systems in a spin-adapted representation, using a linearly-parameterised Ansatz for the many-body wave function. We show that a more compact ground state wave function representation is obtained when combining the symmetric group, $\mathcal{S}_n$, in the form of permutations of the underlying lattice site ordering, with the cumulative spin-coupling based on the unitary group, $\mathrm{U}(n)$. In one-dimensional systems the observed compression of the wave function is reminiscent of block-spin renormalization group approaches, and allows us to study larger lattices (here taken up to 80 sites) with the spin-adapted full configuration interaction quantum Monte Carlo method, which benefits from the sparsity of the Hamiltonian matrix and the corresponding sampled eigenstates that emerge from the reordering. We find that in an optimal lattice ordering the configuration state function with highest weight already captures with high accuracy the spin-spin correlation function of the exact ground state wave function. This feature is found for more general lattice models, such as the Hubbard model, and ab initio quantum chemical models, in this work exemplified by a one-dimensional hydrogen chain. We also provide numerical evidence that the optimal lattice ordering for the unitary group approach is not generally equivalent to the optimal ordering obtained for methods based on matrix-product states, such as the density-matrix renormalization group approach.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A comparative study using state-of-the-art electronic structure theories on solid hydrogen phases under high pressures

Identifying the atomic structure and properties of solid hydrogen under high pressures is a long-standing problem of high-pressure physics with far-reaching significance in planetary and materials science. Determining the pressure-temperature phase diagram of hydrogen is challenging for experiment and theory due to the extreme conditions and the required accuracy in the quantum mechanical treatment of the constituent electrons and nuclei, respectively. Here, we demonstrate explicitly that coupled cluster theory can serve as a computationally efficient theoretical tool to predict solid hydrogen phases with high accuracy. We present a first principles study of solid hydrogen phases at pressures ranging from 100 to 450 GPa. The computed static lattice enthalpies are compared to state-of-the-art diffusion Monte Carlo results and density functional theory calculations. Our coupled cluster theory results for the most stable phases including C2/c-24 and P21/c-24 are in good agreement with those obtained using diffusion Monte Carlo, with the exception of Cmca-4, which is predicted to be significantly less stable. We discuss the scope of the employed methods and how they can contribute as efficient and complementary theoretical tools to solve the long-standing puzzle of understanding solid hydrogen phases at high pressures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

NECI: N-Electron Configuration Interaction with emphasis on state-of-the-art stochastic methods

We present NECI, a state-of-the-art implementation of the Full Configuration Interaction Quantum Monte Carlo algorithm, a method based on a stochastic application of the Hamiltonian matrix on a sparse sampling of the wave function. The program utilizes a very powerful parallelization and scales efficiently to more than 24000 CPU cores. In this paper, we describe the core functionalities of NECI and recent developments. This includes the capabilities to calculate ground and excited state energies, properties via the one- and two-body reduced density matrices, as well as spectral and Green's functions for ab initio and model systems. A number of enhancements of the bare FCIQMC algorithm are available within NECI, allowing to use a partially deterministic formulation of the algorithm, working in a spin-adapted basis or supporting transcorrelated Hamiltonians. NECI supports the FCIDUMP file format for integrals, supplying a convenient interface to numerous quantum chemistry programs and it is licensed under GPL-3.0.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structural and electronic properties of solid molecular hydrogen from many-electron theories

We study the structural and electronic properties of phase III of solid hydrogen using accurate many-electron theories and compare to state-of-the-art experimental findings. The atomic structures of phase III modelled by C2/c-24 crystals are fully optimized on the level of second-order perturbation theory, demonstrating that previously employed structures optimized on the level of approximate density functionals exhibit errors in the H$_2$ bond lengths that cause significant discrepancies in the computed quasi particle band gaps and vibrational frequencies compared to experiment. Using the newly optimized atomic structures, we study the band gap closure and change in vibrational frequencies as a function of pressure. Our findings are in good agreement with recent experimental observations and may prove useful in resolving long-standing discrepancies between experimental estimates of metallization pressures possibly caused by disagreeing pressure calibrations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Adaptive Shift Method in Full Configuration Interaction Quantum Monte Carlo: Development and Applications

In a recent paper, we proposed the adaptive shift method for correcting the undersampling bias of the initiator-FCIQMC. The method allows faster convergence with the number of walkers to the FCI limit than the normal initiator method, particularly for large systems. In its application to strongly correlated molecules, however, the method is prone to overshooting the FCI energy at intermediate walker numbers, with convergence to the FCI limit from below. In this paper, we present a solution to the overshooting problem in strongly correlated molecules, as well as further accelerating convergence to the FCI energy. This is achieved by offsetting the reference energy to a value typically below the Hartree-Fock energy but above the exact energy. This offsetting procedure does not change the exactness property of the algorithm, namely convergence to the exact FCI solution in the large-walker limit, but at its optimal value greatly accelerates convergence. There is no overhead cost associated with this offsetting procedure, and is therefore a pure and substantial computational gain. We illustrate the behavior of this offset adaptive shift method by applying it to the N$_2$ molecule, the ozone molecule at three different geometries (equilibrium open minimum, a hypothetical ring minimum, and a transition state) in three basis sets (cc-pV$X$Z, $X$=D,T,Q), and the chromium dimer in cc-pVDZ basis set, correlating 28 electrons in 76 orbitals. We show that in most cases the offset adaptive shift method converges much faster than both the normal initiator method and the original adaptive shift method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The color center singlet state of oxygen vacancies in TiO$_2$

Oxygen vacancies are ubiquitous in TiO$_2$ and play key roles in catalysis and magnetism applications.Despite being extensively investigated, the electronic structure of oxygen vacancies in TiO$_2$ remains controversial both experimentally and theoretically.Here we report a study of a neutral oxygen vacancy in TiO$_2$ using state-of-the-art quantum chemical electronic structure methods.We find that the ground state is a color center singlet state in both the rutile and the anatase phase of TiO$_2$. Specifically, embedded CCSD(T) calculations find, for an oxygen vacancy in rutile, that the lowest triplet state energy is 0.6 eV above the singlet state, and in anatase the triplet state energy is higher by 1.4 eV. Our study provides fresh insights on the electronic structure of the oxygen vacancy in TiO$_2$, clarifying earlier controversies and potentially inspiring future studies of defects with correlated wave function theories.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Small polarons and the Janus nature of $\text{TiO}_\text{2}(110)$

Polarons are ubiquitous in many semiconductors and have been linked with conductivity and optical response of materials for photovoltaics and heterogeneous catalysis, yet how surface polarons influence adsorption remains unclear. Here, by modelling the surface of rutile titania using density functional theory, we reveal the effect of small surface polarons on water adsorption, dissociation, and hydrogen bonding. On the one hand the presence of such polarons significantly suppresses dissociation of water molecules that are bonded directly to polaronic sites. On the other hand, polarons facilitate water dissociation at certain non-polaronic sites. Furthermore, polarons strengthen hydrogen bonds, which in turn affects water dissociation in hydrogen bonded overlayer structures. This study reveals that polarons at the rutile surface have complex, multi-faceted, effects on water adsorption, dissociation and hydrogen bonding, highlighting the importance of polarons on water structure and dynamics on such surfaces. We expect that many of the physical properties of surface polarons identified here will apply more generally to surfaces and interfaces that can host small polarons, beyond titania.

0 citations0 reviewsNo code linkedOpen access