Abstract: The first step in the construction of a regression model or a data-driven analysis, aiming to predict or elucidate the relationship between the atomic scale structure of matter and its properties, involves transforming the Cartesian coordinates of the atoms into a suitable representation. The development of atomic-scale representations has played, and continues to play, a central role in the success of machine-learning methods for chemistry and materials science. This review summarizes the current understanding of the nature and characteristics of the most commonly used structural and chemical descriptions of atomistic structures, highlighting the deep underlying connections between different frameworks, and the ideas that lead to computationally efficient and universally applicable models. It emphasizes the link between properties, structures, their physical chemistry and their mathematical description, provides examples of recent applications to a diverse set of chemical and materials science problems, and outlines the open questions and the most promising research directions in the field.

Abstract: Machine-learned force fields (ML-FFs) combine the accuracy of ab initio methods with the efficiency of conventional force fields. However, current ML-FFs typically ignore electronic degrees of freedom, such as the total charge or spin state, and assume chemical locality, which is problematic when molecules have inconsistent electronic states, or when nonlocal effects play a significant role. This work introduces SpookyNet, a deep neural network for constructing ML-FFs with explicit treatment of electronic degrees of freedom and quantum nonlocality. Chemically meaningful inductive biases and analytical corrections built into the network architecture allow it to properly model physical limits. SpookyNet improves upon the current state-of-the-art (or achieves similar performance) on popular quantum chemistry data sets. Notably, it is able to generalize across chemical and conformational space and can leverage the learned chemical insights, e.g. by predicting unknown spin states, thus helping to close a further important remaining gap for today's machine learning models in quantum chemistry.

Abstract: Implicit solvation is an effective, highly coarse-grained approach in atomic-scale simulations to account for a surrounding liquid electrolyte on the level of a continuous polarizable medium. Originating in molecular chemistry with finite solutes, implicit solvation techniques are now increasingly used in the context of first-principles modeling of electrochemistry and electrocatalysis at extended (often metallic) electrodes. The prevalent ansatz to model the latter electrodes and the reactive surface chemistry at them through slabs in periodic boundary condition supercells brings its specific challenges. Foremost this concerns the diffculty to describe the entire double layer forming at the electrified solid-liquid interface (SLI) within supercell sizes tractable by commonly employed density-functional theory (DFT). We review liquid solvation methodology from this specific application angle, highlighting in particular its use in the widespread {\em ab initio} thermodynamics approach to surface catalysis. Notably, implicit solvation can be employed to mimic a polarization of the electrode's electronic density under the applied potential and the concomitant capacitive charging of the entire double layer beyond the limitations of the employed DFT supercell. Most critical for continuing advances of this effective methodology for the SLI context is the lack of pertinent (experimental or high-level theoretical) reference data needed for parametrization.

Abstract: The photoinduced excited-state dynamics of the keto and enol forms of cytosine is investigated using ab initio surface hopping in order to understand the outcome of molecular beam femtosecond pump-probe photoionization spectroscopy experiments. Both singlet and triplet states are included in the dynamics. The results show that triplet states play a significant role in the relaxation of the keto tautomer, while they are less important in the enol tautomer. In both forms, the T$_1$ state minimum is found too low in energy to be detected in standard photoionization spectroscopy experiments and therefore experimental decay times should arise from a simultaneous relaxation to the ground state and additional intersystem crossing followed by internal conversion to the T$_1$ state. In agreement with available experimental lifetimes, we observe three decay constants of 7 fs, 270 fs and 1900 fs - the first two coming from the keto tautomer and the longer one from the enol tautomer. Deactivation of the enol form is due to internal conversion to the ground state via two S$_1$/S$_0$ conical intersections of ethylenic type.

Abstract: Machine learning (ML) methods are being used in almost every conceivable area of electronic structure theory and molecular simulation. In particular, ML has become firmly established in the construction of high-dimensional interatomic potentials. Not a day goes by without another proof of principle being published on how ML methods can represent and predict quantum mechanical properties - be they observable, such as molecular polarizabilities, or not, such as atomic charges. As ML is becoming pervasive in electronic structure theory and molecular simulation, we provide an overview of how atomistic computational modeling is being transformed by the incorporation of ML approaches. From the perspective of the practitioner in the field, we assess how common workflows to predict structure, dynamics, and spectroscopy are affected by ML. Finally, we discuss how a tighter and lasting integration of ML methods with computational chemistry and materials science can be achieved and what it will mean for research practice, software development, and postgraduate training.

Abstract: Performing electronic structure calculations is a canonical many-body problem that has recently emerged as a challenging new paradigm for neural network quantum states (NNQS). Here, we parameterise the electronic wavefunction with a novel autoregressive neural network (ARN) that permits highly efficient and scalable sampling, whilst also embedding physical priors that reflect the structure of molecular systems without sacrificing expressibility. This allows us to perform electronic structure calculations on molecules with up to 30 spin-orbitals - which consider multiple orders of magnitude more Slater determinants than previous applications of conventional NNQS - and we find that our ansatz can outperform the de-facto gold-standard coupled cluster methods even in the presence of strong quantum correlations. With a highly expressive neural network for which sampling is no longer a computational bottleneck, we conclude that the barriers to further scaling are not associated with the wavefunction ansatz itself, but rather are inherent to any variational Monte Carlo approach.

Abstract: Light-induced chemical processes are ubiquitous in nature and have widespread technological applications. For example, the photoisomerization of azobenzene allows a drug with an azo scaffold to be activated with light. In principle, photoswitches with useful reactive properties, such as high isomerization yields, can be identified through virtual screening with reactive simulations. In practice these simulations are rarely used for screening, since they require hundreds of trajectories and expensive quantum chemical methods to account for non-adiabatic excited state effects. Here we introduce a neural network potential to accelerate such simulations for azobenzene derivatives. The model, which is based on diabatic states, is called the \textit{diabatic artificial neural network} (DANN). The network is six orders of magnitude faster than the quantum chemistry method used for training. DANN is transferable to molecules outside the training set, predicting quantum yields for unseen species that are correlated with experiment. We use the model to virtually screen 3,100 hypothetical molecules, and identify several species with extremely high quantum yields. Our results pave the way for fast and accurate virtual screening of photoactive compounds.

Abstract: There has been recent interest in the deployment of ab initio density matrix renormalization group computations on high performance computing platforms. Here, we introduce a reformulation of the conventional distributed memory ab initio DMRG algorithm that connects it to the conceptually simpler and advantageous sum of sub-Hamiltonians approach. Starting from this framework, we further explore a hierarchy of parallelism strategies, that includes (i) parallelism over the sum of sub-Hamiltonians, (ii) parallelism over sites, (iii) parallelism over normal and complementary operators, (iv) parallelism over symmetry sectors, and (v) parallelism over dense matrix multiplications. We describe how to reduce processor load imbalance and the communication cost of the algorithm to achieve higher efficiencies. We illustrate the performance of our new open-source implementation on a recent benchmark ground-state calculation of benzene in an orbital space of 108 orbitals and 30 electrons, with a bond dimension of up to 6000, and a model of the FeMo cofactor with 76 orbitals and 113 electrons. The observed parallel scaling from 448 to 2800 CPU cores is nearly ideal.

Abstract: Retrosynthesis—the process of identifying a set of reactants to synthesize a target molecule—is of vital importance to material design and drug discovery. Existing machine learning approaches based on language models and graph neural networks have achieved encouraging results. However, the inner connections of these models are rarely discussed, and rigorous evaluations of these models are largely in need. In this paper, we propose a framework that unifies sequence- and graph-based methods as energy-based models (EBMs) with different energy functions. This unified point of view establishes connections between different models and identifies the differences between them, thereby promoting the understanding of model design. We also provide a comprehensive assessment of performance to the community. Moreover, we present a novel “dual” variant within the framework that performs consistent training over Bayesian forward- and backward-prediction by constraining the agreement between the two directions. This model improves the state of the art for template-free approaches where the reaction type is unknown and known.

Abstract: The conversion of photon energy into other energetic forms in molecules is accompanied by charge moving on ultrafast timescales. We directly observe the charge motion at a specific site in an electronically excited molecule using time-resolved x-ray photoelectron spectroscopy (TR-XPS). We extend the concept of static chemical shift from conventional XPS by the excited-state chemical shift (ESCS), which is connected to the charge in the framework of a potential model. This allows us to invert TR-XPS spectra to the dynamic charge at a specific atom. We demonstrate the power of TR-XPS by using sulphur 2p-core-electron-emission probing to study the UV-excited dynamics of 2-thiouracil. The new method allows us to discover that a major part of the population relaxes to the molecular ground state within 220-250 fs. In addition, a 250-fs oscillation, visible in the kinetic energy of the TR-XPS, reveals a coherent exchange of population among electronic states.

Abstract: Machine learning (ML)-accelerated discovery requires large amounts of high-fidelity data to reveal predictive structure-property relationships. For many properties of interest in materials discovery, the challenging nature and high cost of data generation has resulted in a data landscape that is both scarcely populated and of dubious quality. Data-driven techniques starting to overcome these limitations include the use of consensus across functionals in density functional theory, the development of new functionals or accelerated electronic structure theories, and the detection of where computationally demanding methods are most necessary. When properties cannot be reliably simulated, large experimental data sets can be used to train ML models. In the absence of manual curation, increasingly sophisticated natural language processing and automated image analysis are making it possible to learn structure-property relationships from the literature. Models trained on these data sets will improve as they incorporate community feedback.

Abstract: In search of the potential cathode materials for sodium-ion batteries and to understand the diffusion kinetics, we report the detailed analysis of electrochemical investigation of honeycomb structured Na$_{2}$Ni$_{2}$TeO$_{6}$ material using cyclic voltammetry (CV), electrochemical impedance spectroscopy (EIS), galvanostatic charge-discharge (GCD) and galvanostatic intermittent titration technique (GITT). We found the discharge capacities of 82 and 77 mAhg$^{-1}$ at 0.05~C and 0.1~C current rates, respectively, and the mid-working potential of $\approx$3.9~V at 1~C and high capacity retention of 80\% after 500 cycles at 0.5~C as well as excellent rate capability. The analysis of CV data at different scan rates reveals the pseudo-capacitive mechanism of sodium-ion storage. Interestingly, the {\it in-situ} EIS measurements show a systematic change in the charge-transfer resistance at different charge/discharge stages as well as after different number of cycles. The diffusion coefficient extracted using CV, EIS and GITT lies mainly in the range of 10$^{-10}$ to 10$^{-12}$ cm$^{2}$s$^{-1}$ and the de-insertion/insertion of Na$^+$-ion concentration during electrochemical cycling is consistent with the ratio of Ni$^{3+}$/Ni$^{2+}$ valence state determined by photoemission study. Moreover, the post-cyclic results of retrieved active material show very stable structure and morphology even after various charge-discharge cycles. Our detailed electrochemical investigation and diffusion kinetics studies establish the material as a high working potential and long life electrode for sodium-ion batteries.

Abstract: We present an extension of Alchemical Transfer Method (ATM) for the estimation of relative binding free energies of molecular complexes applicable to conventional as well as scaffold-hopping alchemical transformations. The method, named ATM-RBFE, implemented in the free and open-source OpenMM molecular simulation package, aims to provide a simpler and more generally applicable route to the calculation of relative binding free energies than is currently available. The method is based on sound statistical mechanics theory and a novel coordinate perturbation scheme designed to swap the positions of a pair of ligands such that one is transferred from the bulk solvent to the receptor binding site while the other moves simultaneously in the opposite direction. The calculation is conducted directly using a single solvent box prepared using conventional setup tools, without splitting of electrostatic and non-electrostatic transformations, and without pairwise soft-core potentials. ATM-RBFE is validated here against the absolute binding free energies of the SAMPL8 GDCC host-guest benchmark set and against a benchmark set of estrogen receptor $\alpha$ complexes. In each case, the method yields self-consistent and converged relative binding free energy estimates in agreement with absolute binding free energies, reference literature values as well as experimental measurements.

Abstract: The solid-state diffusion coefficient of the electrode active material is one of the key parameters in lithium-ion battery modelling. Conventionally, this diffusion coefficient is estimated through the galvanostatic intermittent titration technique (GITT). In this work, the validity of GITT and a faster alternative technique, intermittent current interruption (ICI), are investigated regarding their effectiveness through a black-box testing approach. A Doyle-Fuller-Newman model with parameters for a LiNi$_{0.8}$Mn$_{0.1}$Co$_{0.1}$O$_2$ electrode is used as a fairly faithful representation as a real battery system, and the GITT and ICI experiments are simulated to extract the diffusion coefficient. With the parameters used in this work, the results show that both the GITT and ICI methods can identify the solid-state diffusion coefficient very well compared to the value used as input into the simulation model. The ICI method allows more frequent measurement but the experiment time is 85% less than what takes to perform a GITT test. Different fitting approaches and fitting length affected the estimation accuracy, however not significantly. Moreover, a thinner electrode, a higher C-rate and a greater electrolyte diffusion coefficient will lead to an estimation of a higher solid-state diffusion coefficient, generally closer to the target value.

Abstract: Spiky/hollow metal nanoparticles have applications across a broad range of fields. However, current bottom-up methods to produce spiky/hollow metal nanoparticles rely heavily on the use of strongly adsorbing surfactant molecules, which is undesirable since these passivate the product particle surfaces. Here we report a high-yield surfactant-free synthesis of spiky hollow Au-Ag nanostars (SHAANs). Each SHAAN is composed of more than 50 spikes attached to a hollow ca. 150 nm diameter cubic core, which makes SHAANs highly plasmonically and catalytically active. Moreover, the surfaces of SHAANs are chemically exposed which gives them significantly enhanced functionality compared to their surfactant-capped counterparts, as demonstrated in surface-enhanced Raman spectroscopy (SERS) and catalysis. The chemical accessibility of the pristine SHAANs also allows the use of hydroxyethyl cellulose as a weakly-bound stabilizing agent. This produces colloidal SHAANs which remain stable for more than 1 month while retaining the functionalities of the pristine particles and allow even single-particle SERS to be realized.

Abstract: Hyperpolarized nuclear magnetic resonance and lab-on-a-chip microfluidics are two dynamic, but until recently quite distinct, fields of research. Recent developments in both areas increased their synergistic overlap. By microfluidic integration, many complex experimental steps can be brought together onto a single platform. Microfluidic devices are therefore increasingly finding applications in medical diagnostics, forensic analysis, and biomedical research. In particular, they provide novel and powerful ways to culture cells, cell aggregates, and even functional models of entire organs. Nuclear magnetic resonance is a non-invasive, high-resolution spectroscopic technique which allows real-time process monitoring with chemical specificity. It is ideally suited for observing metabolic and other biological and chemical processes in microfluidic systems. However, its intrinsically low sensitivity has limited its application. Recent advances in nuclear hyperpolarization techniques may change this: under special circumstances, it is possible to enhance NMR signals by up to 5 orders of magnitude, which dramatically extends the utility of NMR in the context of microfluidic systems. At the same time, hyperpolarization requires complex chemical and/or physical manipulations, which in turn may benefit from microfluidic implementation. In fact, many hyperpolarization methodologies rely on processes that are more efficient at the micro-scale, such as molecular diffusion, penetration electromagnetic radiation into the sample, or restricted molecular mobility on a surface. In this review we examine the confluence between the fields of hyperpolarization-enhanced NMR and microfluidics, and assess how these areas of research have mutually benefited one another, and will continue to do so.

Abstract: Accurate force fields are necessary for predictive molecular simulations. However, developing force fields that accurately reproduce experimental properties is challenging. Here, we present a machine learning directed, multiobjective optimization workflow for force field parameterization that evaluates millions of prospective force field parameter sets while requiring only a small fraction of them to be tested with molecular simulations. We demonstrate the generality of the approach and identify multiple low-error parameter sets for two distinct test cases: simulations of hydrofluorocarbon (HFC) vapor-liquid equilibrium (VLE) and an ammonium perchlorate (AP) crystal phase. We discuss the challenges and implications of our force field optimization workflow.

Abstract: Prediction of a molecule's 3D conformer ensemble from the molecular graph holds a key role in areas of cheminformatics and drug discovery. Existing generative models have several drawbacks including lack of modeling important molecular geometry elements (e.g. torsion angles), separate optimization stages prone to error accumulation, and the need for structure fine-tuning based on approximate classical force-fields or computationally expensive methods such as metadynamics with approximate quantum mechanics calculations at each geometry. We propose GeoMol--an end-to-end, non-autoregressive and SE(3)-invariant machine learning approach to generate distributions of low-energy molecular 3D conformers. Leveraging the power of message passing neural networks (MPNNs) to capture local and global graph information, we predict local atomic 3D structures and torsion angles, avoiding unnecessary over-parameterization of the geometric degrees of freedom (e.g. one angle per non-terminal bond). Such local predictions suffice both for the training loss computation, as well as for the full deterministic conformer assembly (at test time). We devise a non-adversarial optimal transport based loss function to promote diverse conformer generation. GeoMol predominantly outperforms popular open-source, commercial, or state-of-the-art machine learning (ML) models, while achieving significant speed-ups. We expect such differentiable 3D structure generators to significantly impact molecular modeling and related applications.

Abstract: We develop a quartic-scaling implementation of coupled-cluster singles and doubles based on low-rank tensor hypercontraction (THC) factorizations of both the electron repulsion integrals (ERIs) and the doubles amplitudes. This extends our rank-reduced coupled-cluster method to incorporate higher-order tensor factorizations. The THC factorization of the doubles amplitudes accounts for most of the gain in computational efficiency as it is sufficient, in conjunction with a Cholesky decomposition of the ERIs, to reduce the computational complexity of most contributions to the CCSD amplitude equations. Further THC factorization of the ERIs reduces the complexity of certain terms arising from nested commutators between the doubles excitation operator and the two-electron operator. We implement this new algorithm using graphical processing units (GPUs) and demonstrate that it enables CCSD calculations for molecules with 250 atoms and 2500 basis functions using a single computer node. Further, we show that the new method computes correlation energies with comparable accuracy to the underlying RR-CCSD method.

Abstract: Physically-motivated and mathematically robust atom-centred representations of molecular structures are key to the success of modern atomistic machine learning (ML) methods They lie at the foundation of a wide range of methods to predict the properties of both materials and molecules as well as to explore and visualize the chemical compound and configuration space Recently, it has become clear that many of the most effective representations share a fundamental formal connection: that they can all be expressed as a discretization of N-body correlation functions of the local atom density, suggesting the opportunity of standardizing and, more importantly, optimizing the calculation of such representations We present an implementation, named librascal, whose modular design lends itself both to developing refinements to the density-based formalism and to rapid prototyping for new developments of rotationally equivariant atomistic representations As an example, we discuss SOAP features, perhaps the most widely used member of this family of representations, to show how the expansion of the local density can be optimized for any choice of radial basis set We discuss the representation in the context of a kernel ridge regression model, commonly used with SOAP features, and analyze how the computational effort scales for each of the individual steps of the calculation By applying data reduction techniques in feature space, we show how to further reduce the total computational cost by at up to a factor of 4 or 5 without affecting the model's symmetry properties and without significantly impacting its accuracy

Abstract: Signal Amplification By Reversible Exchange (SABRE) creates hyperpolarization (large spin magnetization) using a transition metal catalyst and parahydrogen, addressing the sensitivity limitations of magnetic resonance. SABRE and its heteronuclear variant X-SABRE are simple, fast and general, but to date have not produced polarization levels as large as more established methods. We show here that inaccuracies in the theoretical framework for these applications, which focuses on avoided crossings (also called level anti-crossings or LACs) steer current SABRE and X-SABRE experiments away from optimal solutions. Accurate simulations show astonishingly rich and unexpected dynamics in SABRE/X-SABRE, which we explain with a combination of perturbation theory and average Hamiltonian approaches. This theoretical picture predicts simple pulse sequences with field values far from LACs (both instantaneously and on average), using different terms in the effective Hamiltonian to strategically control evolution and improve polarization transfer. Significant signal enhancements under such highly non-intuitive conditions are verified experimentally.

Abstract: Non-resonant second harmonic generation phase and amplitude measurements obtained from the silica:water interface at varying pH and 05 M ionic strength point to the existence of a nonlinear susceptibility term, which we call chi(3)X, that is associated with a 90 deg phase shift Including this contribution in a model for the total effective second-order nonlinear susceptibility produces reasonable point estimates for interfacial potentials and second-order nonlinear susceptibilities when chi(3)Xis about 15 times chi(3)water A model without this term and containing only traditional chi(2) and chi(3) terms cannot recapitulate the experimental data The new model also provides a demonstrated utility for distinguishing apparent differences in the second-order nonlinear susceptibility when the electrolyte is NaCl vs MgSO4, pointing to the possibility of using HD-SHG to investigate ion-specificity in interfacial processes

Abstract: The exact factorization of the time-dependent electron-nuclear wavefunction has been employed successfully in the field of quantum molecular dynamics simulations for interpreting and simulating light-induced ultrafast processes. In this work, we summarize the major developments leading to the formulation of a trajectory-based approach, derived from the exact factorization equations, capable of dealing with nonadiabatic electronic processes, and including spin-orbit coupling and the non-perturbative effect of an external time-dependent field. This trajectory-based quantum-classical approach has been dubbed coupled-trajectory mixed quantum-classical (CT-MQC) algorithm, whose performance is tested here to study the photo-dissociation dynamics of IBr.

Abstract: A variety of homologous carbon chains (HCnH, HCnN, CnS, CnO, and OCnO) are found to exhibit an appealing even-odd effect. Chains containing a number of carbon atoms of a certain parity possess singlet ground states, while members of opposite parity have triplet ground states. From a general perspective, it is important that this even-odd effect confounds straightforward chemical intuition. Whether the most stable form is a triplet or a singlet is neither simply related to the fact that the species in question is a normal (closed-shell, nonradical) molecule nor a (di)radical or to the (e.g., cumulene-type) C-C bond succession across the chain. From a computational perspective, the present results are important also because they demonstrate that electron correlations in carbon-based chains are extremely strong. Whether the gold-standard CCSD(T) (coupled-cluster expansions with single and double excitations and triple excitations corrections) framework suffices to describe such strongly correlated systems remains an open question that calls for further clarification. Most importantly for astrochemistry, the present results may explain why certain members are not astronomically observed although larger members of the same homologous series are detected; the missing species are exactly those for which the present calculations predict triplet ground states.

Abstract: We report a new deep learning message passing network that takes inspiration from Newton's equations of motion to learn interatomic potentials and forces. With the advantage of directional information from trainable latent force vectors, and physics-infused operators that are inspired by the Newtonian physics, the entire model remains rotationally equivariant, and many-body interactions are inferred by more interpretable physical features. We test NewtonNet on the prediction of several reactive and non-reactive high quality ab initio data sets including single small molecule dynamics, a large set of chemically diverse molecules, and methane and hydrogen combustion reactions, achieving state-of-the-art test performance on energies and forces with far greater data and computational efficiency than other deep learning models.

Abstract: Recent experimental advances in strongly coupled light-matter systems has sparked the development of general ab-initio methods capable of describing interacting light-matter systems from first principles. One of these methods, quantum-electrodynamical density-functional theory (QEDFT), promises computationally efficient calculations for large correlated light-matter systems with the quality of the calculation depending on the underlying approximation for the exchange-correlation functional. So far no true density-functional approximation has been introduced limiting the efficient application of the theory. In this paper, we introduce the first gradient-based density functional for the QEDFT exchange-correlation energy derived from the adiabatic-connection fluctuation-dissipation theorem. We benchmark this simple-to-implement approximation on small systems in optical cavities and demonstrate its relatively low computational costs for fullerene molecules up to C$_{180}$ coupled to 400,000 photon modes in a dissipative optical cavity. This work now makes first principle calculations of much larger systems possible within the QEDFT framework effectively combining quantum optics with large-scale electronic structure theory.

Abstract: We present an efficient implementation of periodic Gaussian density fitting (GDF) using the Coulomb metric The three-center integrals are divided into two parts by range-separating the Coulomb kernel, with the short-range part evaluated in real space and the long-range part in reciprocal space With a few algorithmic optimizations, we show that this new method -- which we call range-separated GDF (RSGDF) -- scales sublinearly to linearly with the number of $k$-points for small to medium-sized $k$-point meshes that are commonly used in periodic calculations with electron correlation Numerical results on a few three-dimensional solids show about $10$-fold speedups over the previously developed GDF with little precision loss The error introduced by RSGDF is about $10^{-5}~E_{\textrm{h}}$ in the converged Hartree-Fock energy with default auxiliary basis sets and can be systematically reduced by increasing the size of the auxiliary basis with little extra work [The article has been accepted by The Journal of Chemical Physics]

Abstract: Quantum machine learning algorithms have emerged to be a promising alternative to their classical counterparts as they leverage the power of quantum computers. Such algorithms have been developed to solve problems like electronic structure calculations of molecular systems and spin models in magnetic systems. However the discussion in all these recipes focus specifically on targeting the ground state. Herein we demonstrate a quantum algorithm that can filter any energy eigenstate of the system based on either symmetry properties or on a predefined choice of the user. The work horse of our technique is a shallow neural network encoding the desired state of the system with the amplitude computed by sampling the Gibbs- Boltzmann distribution using a quantum circuit and the phase information obtained classically from the non-linear activation of a separate set of neurons. We show that the resource requirements of our algorithm is strictly quadratic. To demonstrate its efficacy, we use state-filtration in monolayer transition metal-dichalcogenides which are hitherto unexplored in any flavor of quantum simulations. We implement our algorithm not only on quantum simulators but also on actual IBM-Q quantum devices and show good agreement with the results procured from conventional electronic structure calculations. We thus expect our protocol to provide a new alternative in exploring band-structures of exquisite materials to usual electronic structure methods or machine learning techniques which are implementable solely on a classical computer

Abstract: Following our recent work on the benzene molecule [\href{this https URL}{J.~Chem.~Phys.~\textbf{153}, 176101 (2020)}], itself motivated by the blind challenge of Eriksen \textit{et al.} [\href{this https URL}{J.~Phys.~Chem.~Lett.~\textbf{11}, 8922 (2020)}] on the same system, we report accurate full configuration interaction (FCI) frozen-core correlation energy estimates for twelve five- and six-membered ring molecules in the standard correlation-consistent double-$\zeta$ Dunning basis set (cc-pVDZ). Our FCI correlation energy estimates, with estimated error smaller than 1 millihartree, are based on energetically optimized-orbital selected configuration interaction (SCI) calculations performed with the \textit{Configuration Interaction using a Perturbative Selection made Iteratively} (CIPSI) algorithm. Having at our disposal these accurate reference energies, the respective performance and convergence properties of several popular and widely-used families of single-reference quantum chemistry methods are investigated. In particular, we study the convergence properties of i) the Moller-Plesset perturbation series up to fifth-order (MP2, MP3, MP4, and MP5), ii) the iterative approximate coupled-cluster series CC2, CC3, and CC4, and iii) the coupled-cluster series CCSD, CCSDT, and CCSDTQ. The performance of the ground-state gold standard CCSD(T) as well as the completely renormalized CC model, CR-CC(2,3), are also investigated.