-
Emergent Mesoscale Correlations in Active Solids with Noisy Chiral Dynamics
Authors:
Amir Shee,
Silke Henkes,
Cristián Huepe
Abstract:
We present the linear response theory for an elastic solid composed of active Brownian particles with intrinsic individual chirality, deriving both a normal mode formulation and a continuum elastic formulation. Using this theory, we compute analytically the velocity correlations and energy spectra under different conditions, showing an excellent agreement with simulations. We generate the correspo…
▽ More
We present the linear response theory for an elastic solid composed of active Brownian particles with intrinsic individual chirality, deriving both a normal mode formulation and a continuum elastic formulation. Using this theory, we compute analytically the velocity correlations and energy spectra under different conditions, showing an excellent agreement with simulations. We generate the corresponding phase diagram, identifying chiral and achiral disordered regimes (for high chirality or noise levels), as well as chiral and achiral states with mesoscopic-range order (for low chirality and noise). The chiral ordered states display mesoscopic spatial correlations and oscillating time correlations, but no wave propagation. In the high chirality regime, we find a peak in the elastic energy spectrum that leads to a non-monotonic behavior with increasing noise strength that is consistent with the emergence of the `hammering' state recently identified in chiral glasses. Finally, we show numerically that our theory, despite its linear response nature, can be applied beyond the idealized homogeneous solid assumed in our derivations. Indeed, by increasing the level of activity, we show that it remains a good approximation of the system dynamics until just below the melting transition. In addition, we show that there is still an excellent agreement between our analytical results and simulations when we extend our results to heterogeneous solids composed of mixtures of active particles with different intrinsic chirality and noise levels. Our work provides a thorough analytical and numerical description of the emergent states in a densely packed system of chiral self-propelled Brownian disks, thus allowing a detailed understanding of the phases and dynamics identified in a minimal chiral active systems.
△ Less
Submitted 6 July, 2024;
originally announced July 2024.
-
Performance of wave function and Green's functions based methods for non equilibrium many-body dynamics
Authors:
Cian C. Reeves,
Gaurav Harsha,
Avijit Shee,
Yuanran Zhu,
Chao Yang,
K Birgitta Whaley,
Dominika Zgid,
Vojtech Vlcek
Abstract:
Theoretical descriptions of non equilibrium dynamics of quantum many-body systems essentially employ either (i) explicit treatments, relying on truncation of the expansion of the many-body wave function, (ii) compressed representations of the many-body wave function, or (iii) evolution of an effective (downfolded) representation through Green's functions. In this work, we select representative cas…
▽ More
Theoretical descriptions of non equilibrium dynamics of quantum many-body systems essentially employ either (i) explicit treatments, relying on truncation of the expansion of the many-body wave function, (ii) compressed representations of the many-body wave function, or (iii) evolution of an effective (downfolded) representation through Green's functions. In this work, we select representative cases of each of the methods and address how these complementary approaches capture the dynamics driven by intense field perturbations to non equilibrium states. Under strong driving, the systems are characterized by strong entanglement of the single particle density matrix and natural populations approaching those of a strongly interacting equilibrium system. We generate a representative set of results that are numerically exact and form a basis for critical comparison of the distinct families of methods. We demonstrate that the compressed formulation based on similarity transformed Hamiltonians (coupled cluster approach) is practically exact in weak fields and, hence, weakly or moderately correlated systems. Coupled cluster, however, struggles for strong driving fields, under which the system exhibits strongly correlated behavior, as measured by the von Neumann entropy of the single particle density matrix. The dynamics predicted by Green's functions in the (widely popular) GW approximation are less accurate by improve significantly upon the mean-field results in the strongly driven regime.
△ Less
Submitted 14 May, 2024;
originally announced May 2024.
-
Impact of chirality on active Brownian particle: Exact moments in two and three dimensions
Authors:
Anweshika Pattanayak,
Amir Shee,
Debasish Chaudhuri,
Abhishek Chaudhuri
Abstract:
In this work, we investigate the effects of chirality, accounting for translational diffusion, on active Brownian particles in two and three dimensions. Despite the inherent complexity in solving the Fokker-Planck equation, we demonstrate a Laplace transform method for precisely calculating the temporal evolution of various dynamic moments. Our analysis yields explicit expressions for multiple mom…
▽ More
In this work, we investigate the effects of chirality, accounting for translational diffusion, on active Brownian particles in two and three dimensions. Despite the inherent complexity in solving the Fokker-Planck equation, we demonstrate a Laplace transform method for precisely calculating the temporal evolution of various dynamic moments. Our analysis yields explicit expressions for multiple moments, such as the second and fourth moments of displacement, revealing the impact of persistence and chirality. These moments exhibit oscillatory behaviour, and excess kurtosis indicates deviations from the Gaussian distribution during intermediate time intervals.
△ Less
Submitted 15 April, 2024;
originally announced April 2024.
-
A static quantum embedding scheme based on coupled cluster theory
Authors:
Avijit Shee,
Fabian M. Faulstich,
Birgitta Whaley,
Lin Lin,
Martin Head-Gordon
Abstract:
We develop a static quantum embedding scheme, utilizing projection equations to solve coupled cluster (CC) amplitudes. To reduce the computational cost (for example, of a large basis set calculation), we solve the local fragment problem using a high-level coupled cluster method and address the environment problem with a lower-level Møller-Plesset (MP) perturbative method. This embedding approach i…
▽ More
We develop a static quantum embedding scheme, utilizing projection equations to solve coupled cluster (CC) amplitudes. To reduce the computational cost (for example, of a large basis set calculation), we solve the local fragment problem using a high-level coupled cluster method and address the environment problem with a lower-level Møller-Plesset (MP) perturbative method. This embedding approach is consistently formulated within the coupled cluster framework and will be called MP-CC. We demonstrate the effectiveness of our method through several prototypical molecular examples by analyzing a global quantity, that is, the total correlation energy of the system in the study of potential energy curves (PEC) and thermochemical reaction energies. We have shown that our method can achieve comparable accuracy both with a small and large basis set when the fragment Hilbert space size remains the same. Additionally, our results indicate that increasing the fragment size can systematically enhance the accuracy of observables, approaching the precision of the full coupled cluster solver.
△ Less
Submitted 13 April, 2024;
originally announced April 2024.
-
Self-Aligning Polar Active Matter
Authors:
Paul Baconnier,
Olivier Dauchot,
Vincent Démery,
Gustavo Düring,
Silke Henkes,
Cristián Huepe,
Amir Shee
Abstract:
Self-alignment describes the property of a polar active unit to align or anti-align its orientation towards its velocity. In contrast to mutual alignment, where the headings of multiple active units tend to directly align to each other -- as in the celebrated Vicsek model --, self-alignment impacts the dynamics at the individual level by coupling the rotation and displacements of each active unit.…
▽ More
Self-alignment describes the property of a polar active unit to align or anti-align its orientation towards its velocity. In contrast to mutual alignment, where the headings of multiple active units tend to directly align to each other -- as in the celebrated Vicsek model --, self-alignment impacts the dynamics at the individual level by coupling the rotation and displacements of each active unit. This enriches the dynamics even without interactions or external forces, and allows, for example, a single self-propelled particle to orbit in a harmonic potential. At the collective level, self-alignment modifies the nature of the transition to collective motion already in the mean field description, and it can also lead to other forms of self-organization such as collective actuation in dense or solid elastic assemblies of active units. This has significant implications for the study of dense biological systems, metamaterials, and swarm robotics. Here, we review a number of models that were introduced independently to describe the previously overlooked property of self-alignment and identify some of its experimental realizations. Our aim is three-fold: (i)~underline the importance of self-alignment in active systems, especially in the context of dense populations of active units and active solids; (ii)~provide a unified mathematical and conceptual framework for the description of self-aligning systems; (iii)~discuss the common features and specific differences of the existing models of self-alignment. We conclude by discussing promising research avenues in which the concept of self-alignment could play a significant role.
△ Less
Submitted 15 March, 2024;
originally announced March 2024.
-
Comparing self-consistent GW and vertex corrected G0W0 (G0W0Γ) accuracy for molecular ionization potentials
Authors:
Ming Wen,
Vibin Abraham,
Gaurav Harsha,
Avijit Shee,
K. Birgitta Whaley,
Dominika Zgid
Abstract:
We test the performance of self-consistent GW and several representative implementations of vertex corrected G0W0 (G0W0Γ). These approaches are tested on benchmark data sets covering full valence spectra (first ionization potentials and some inner valence shell excitations). For small molecules, when comparing against state of the art wave function techniques, our results show that performing full…
▽ More
We test the performance of self-consistent GW and several representative implementations of vertex corrected G0W0 (G0W0Γ). These approaches are tested on benchmark data sets covering full valence spectra (first ionization potentials and some inner valence shell excitations). For small molecules, when comparing against state of the art wave function techniques, our results show that performing full self-consistency in the GW scheme either systematically outperforms vertex corrected G0W0 or gives results of at least the same quality. Moreover, the G0W0Γ results in additional computational cost when compared to G0W0 or self-consistent GW and the G0W0Γ dependency on the starting mean-filed solution is frequently larger than the magnitude of the vertex correction. Consequently, for molecular systems self-consistent GW performed on imaginary axis and then followed by modern analytical continuation techniques offers a more reliable approach to make predictions of IP spectra.
△ Less
Submitted 5 April, 2024; v1 submitted 20 November, 2023;
originally announced November 2023.
-
Noise-induced quenched disorder in dense active systems
Authors:
Guozheng Lin,
Zhangang Han,
Amir Shee,
Cristián Huepe
Abstract:
We report and characterize the emergence of a noise-induced state of quenched disorder in a generic model describing a dense sheet of active polar disks with non-isotropic rotational and translational dynamics. In this state, randomly oriented self-propelled disks become jammed, only displaying small fluctuations about their mean positions and headings. The quenched disorder phase appears at inter…
▽ More
We report and characterize the emergence of a noise-induced state of quenched disorder in a generic model describing a dense sheet of active polar disks with non-isotropic rotational and translational dynamics. In this state, randomly oriented self-propelled disks become jammed, only displaying small fluctuations about their mean positions and headings. The quenched disorder phase appears at intermediate noise levels, between the two states that typically define the flocking transition (a standard disordered state that displays continuously changing headings due to rotational diffusion and a polar order state of collective motion). We find that the angular fluctuations in this dense system follow an Ornstein-Uhlenbeck process leading to retrograde forces that oppose self-propulsion, and determine its properties. Using this result, we explain the mechanism behind the emergence of the quenched state and compute analytically its critical noise, showing that it matches our numerical simulations. We argue that this novel type of state could be observed in a broad range of natural and artificial dense active systems with repulsive interactions.
△ Less
Submitted 7 April, 2023;
originally announced April 2023.
-
Triple excitations in Green's function coupled cluster solver for studies of strongly correlated systems in the framework of self-energy embedding theory
Authors:
Avijit Shee,
Chia-Nan Yeh,
Bo Peng,
Karol Kowalski,
Dominika Zgid
Abstract:
Embedding theories became important approaches used for accurate calculations of both molecules and solids. In these theories, a small chosen subset of orbitals is treated with an accurate method, called an impurity solver, capable of describing higher correlation effects. Ideally, such a chosen fragment should contain multiple orbitals responsible for the chemical and physical behavior of the com…
▽ More
Embedding theories became important approaches used for accurate calculations of both molecules and solids. In these theories, a small chosen subset of orbitals is treated with an accurate method, called an impurity solver, capable of describing higher correlation effects. Ideally, such a chosen fragment should contain multiple orbitals responsible for the chemical and physical behavior of the compound. Handing a large number of chosen orbitals presents a very significant challenge for the current generation of solvers used in the physics and chemistry community. Here, we develop a Green's function coupled cluster singles doubles and triples (GFCCSDT) solver that can be used for a quantitative description in both molecules and solids. This solver allows us to treat orbital spaces that are inaccessible to other accurate solvers. At the same time, GFCCSDT maintains high accuracy of the resulting self-energy. Moreover, in conjunction with the GFCCSD solver, it allows us to test the systematic convergence of computational studies. Developing the CC family of solvers paves the road to fully systematic Green's function embedding calculations in solids. In this paper, we focus on the investigation of GFCCSDT self-energies for a strongly correlated problem of SrMnO$_3$ solid. Subsequently, we apply this solver to solid MnO showing that an approximate variant of GFCCSDT is capable of yielding a high accuracy orbital resolved spectral function.
△ Less
Submitted 22 November, 2022;
originally announced November 2022.
-
Relativistic Self-Consistent $GW$: Exact Two-Component Formalism with One-Electron Approximation for Solids
Authors:
Chia-Nan Yeh,
Avijit Shee,
Qiming Sun,
Emanuel Gull,
Dominika Zgid
Abstract:
We present a formulation of relativistic self-consistent $GW$ for solids based on the exact two-component formalism with one-electron approximation (X2C1e) and non-relativistic Coulomb interactions. Our theory allows us to study scalar relativistic effects, spin-orbit coupling, and the interplay of relativistic effects with electron correlation without adjustable parameters. Our all-electron imple…
▽ More
We present a formulation of relativistic self-consistent $GW$ for solids based on the exact two-component formalism with one-electron approximation (X2C1e) and non-relativistic Coulomb interactions. Our theory allows us to study scalar relativistic effects, spin-orbit coupling, and the interplay of relativistic effects with electron correlation without adjustable parameters. Our all-electron implementation is fully $ab$ $initio$ and does not require a pseudopotential constructed from atomic calculations. We examine the effect of the X2C1e approximation by comparison to the established four-component formalism and reach excellent agreement. The simplicity of X2C1e enables the construction of higher order theories, such as embedding theories, on top of perturbative calculations.
△ Less
Submitted 19 May, 2022; v1 submitted 4 February, 2022;
originally announced February 2022.
-
Self-propulsion with speed and orientation fluctuation: exact computation of moments and dynamical bistabilities in displacement
Authors:
Amir Shee,
Debasish Chaudhuri
Abstract:
We consider the influence of active speed fluctuations on the dynamics of a $d$-dimensional active Brownian particle performing a persistent stochastic motion. We use the Laplace transform of the Fokker-Planck equation to obtain exact expressions for time-dependent dynamical moments. Our results agree with direct numerical simulations and show several dynamical crossovers determined by the activit…
▽ More
We consider the influence of active speed fluctuations on the dynamics of a $d$-dimensional active Brownian particle performing a persistent stochastic motion. We use the Laplace transform of the Fokker-Planck equation to obtain exact expressions for time-dependent dynamical moments. Our results agree with direct numerical simulations and show several dynamical crossovers determined by the activity, persistence, and speed fluctuation.The persistence in the motion leads to anisotropy, with the parallel component of displacement fluctuation showing sub-diffusive behavior and non-monotonic variation. The kurtosis remains positive at short times determined by the speed fluctuation, crossing over to a negative minimum at intermediate times governed by the persistence before vanishing asymptotically. The probability distribution of particle displacement obtained from numerical simulations in two-dimension shows two crossovers between contracted and expanded trajectories via two bimodal distributions at intervening times. While the speed fluctuation dominates the first crossover, the second crossover is controlled by persistence like in the worm-like chain model of semiflexible polymers.
△ Less
Submitted 26 December, 2021;
originally announced December 2021.
-
Active Brownian motion with speed fluctuations in arbitrary dimensions: exact calculation of moments and dynamical crossovers
Authors:
Amir Shee,
Debasish Chaudhuri
Abstract:
We consider the motion of an active Brownian particle with speed fluctuations in d-dimensions in the presence of both translational and orientational diffusion. We use an Ornstein-Uhlenbeck process for active speed generation. Using a Laplace transform approach, we describe and use a Fokker-Planck equation-based method to evaluate the exact time dependence of all relevant dynamical moments. We pre…
▽ More
We consider the motion of an active Brownian particle with speed fluctuations in d-dimensions in the presence of both translational and orientational diffusion. We use an Ornstein-Uhlenbeck process for active speed generation. Using a Laplace transform approach, we describe and use a Fokker-Planck equation-based method to evaluate the exact time dependence of all relevant dynamical moments. We present explicit calculations of such moments and compare our analytical predictions against numerical simulations to demonstrate and analyze several dynamical crossovers. The kurtosis of displacement shows positive or negative deviations from a Gaussian behavior at intermediate times depending on the dominance of speed or orientational fluctuations.
△ Less
Submitted 27 August, 2021;
originally announced August 2021.
-
Exploring Coupled Cluster Green's function as a method for treating system and environment in Green's function embedding methods
Authors:
Avijit Shee,
Chia-Nan Yeh,
Dominika Zgid
Abstract:
Within the self-energy embedding theory (SEET) framework, we study coupled cluster Green's function (GFCC) method in two different contexts: as a method to treat either the system or environment present in the embedding construction. Our study reveals that when GFCC is used to treat the environment we do not see improvement in total energies in comparison to the coupled cluster method itself. To r…
▽ More
Within the self-energy embedding theory (SEET) framework, we study coupled cluster Green's function (GFCC) method in two different contexts: as a method to treat either the system or environment present in the embedding construction. Our study reveals that when GFCC is used to treat the environment we do not see improvement in total energies in comparison to the coupled cluster method itself. To rationalize this puzzling result, we analyze the performance of GFCC as an impurity solver with a series of transition metal oxides. These studies shed light on strength and weaknesses of such a solver and demonstrate that such a solver gives very accurate results when the size of the impurity is small. We investigate if it is possible to achieve a systematic accuracy of the embedding solution when we increase the size of the impurity problem. We found that in such a case, the performance of the solver worsens, both in terms of finding the ground state solution of the impurity problem as well as the self-energies produced. We concluded that increasing the rank of GFCC solver is necessary to be able to enlarge impurity problems and achieve a reliable accuracy. We also have shown that natural orbitals from weakly correlated perturbative methods are better suited than symmetrized atomic orbitals (SAO) when the total energy of the system is the target quantity.
△ Less
Submitted 16 July, 2021;
originally announced July 2021.
-
Testing the GFCCSD impurity solver on real materials within the self-energy embedding theory framework
Authors:
Chia-Nan Yeh,
Avijit Shee,
Sergei Iskakov,
Dominika Zgid
Abstract:
We apply the Green's function coupled cluster singles and doubles (GFCCSD) impurity solver to realistic impurity problems arising for strongly correlated solids within the self-energy embedding theory (SEET) framework. We describe the details of our GFCC solver implementation, investigate its performance, and highlight potential advantages and problems on examples of impurities created during the…
▽ More
We apply the Green's function coupled cluster singles and doubles (GFCCSD) impurity solver to realistic impurity problems arising for strongly correlated solids within the self-energy embedding theory (SEET) framework. We describe the details of our GFCC solver implementation, investigate its performance, and highlight potential advantages and problems on examples of impurities created during the self-consistent SEET for antiferromagnetic MnO and paramagnetic SrMnO$_{3}$. GFCCSD provides satisfactory descriptions for weakly and moderately correlated impurities with sizes that are intractable by existing accurate impurity solvers such as exact diagonalization (ED). However, our data also shows that when correlations become strong, the singles and doubles approximation used in GFCC could lead to instabilities in searching for the particle number present in impurity problems. These instabilities appears especially severe when the impurity size gets larger and multiple degenerate orbitals with strong correlations are present. We conclude that to fully check the reliability of GFCCSD results and use them in fully {\em ab initio} calculations in the absence of experiments, a verification from a GFCC solver with higher order excitations is necessary.
△ Less
Submitted 31 December, 2020;
originally announced December 2020.
-
Filament-motor protein system under loading: instability and limit cycle oscillations
Authors:
Amir Shee,
Subhadip Ghosh,
Debasish Chaudhuri
Abstract:
We consider the dynamics of a rigid filament in a motor protein assay under external loading. The motor proteins are modeled as active harmonic linkers with tail ends immobilized on a substrate. Their heads attach to the filament stochastically to extend along it, resulting in a force on the filament, before detaching. The rate of extension and detachment are load dependent. Here we formulate and…
▽ More
We consider the dynamics of a rigid filament in a motor protein assay under external loading. The motor proteins are modeled as active harmonic linkers with tail ends immobilized on a substrate. Their heads attach to the filament stochastically to extend along it, resulting in a force on the filament, before detaching. The rate of extension and detachment are load dependent. Here we formulate and characterize the governing dynamics in the mean field approximation using linear stability analysis, and direct numerical simulations of the motor proteins and filament. Under constant loading, the system shows transition from a stable configuration to instability towards detachment of the filament from motor proteins. Under elastic loading, we find emergence of stable limit cycle oscillations via a supercritical Hopf bifurcation with change in activity and the number of motor proteins. Numerical simulations of the system for large number of motor proteins show good agreement with the mean field predictions.
△ Less
Submitted 17 December, 2020;
originally announced December 2020.
-
Relativistic EOM-CCSD for core-excited and core-ionized state energies based on the 4-component Dirac-Coulomb(-Gaunt) Hamiltonian
Authors:
Loïc Halbert,
Marta Lopez Vidal,
Avijit Shee,
Sonia Coriani,
André Severo Pereira Gomes
Abstract:
We report an implementation of the core-valence separation approach to the 4-component relativistic Hamiltonian based equation-of-motion coupled-cluster with singles and doubles theory (CVS-EOM-CCSD), for the calculation of relativistic core-ionization potentials and core-excitation energies. With this implementation, which is capable of exploiting double group symmetry, we investigate the effects…
▽ More
We report an implementation of the core-valence separation approach to the 4-component relativistic Hamiltonian based equation-of-motion coupled-cluster with singles and doubles theory (CVS-EOM-CCSD), for the calculation of relativistic core-ionization potentials and core-excitation energies. With this implementation, which is capable of exploiting double group symmetry, we investigate the effects of the different CVS-EOM-CCSD variants, and the use of different Hamiltonians based on the exact 2-component (X2C) framework, on the energies of different core ionized and excited states in halogen (CH$_3$I, HX and X$^-$, X = Cl-At) and xenon containing (Xe, XeF$_2$) species. Our results show that the X2C molecular mean-field approach [Sikkema et al., J. Chem. Phys. 2009, 131, 124116], based on 4-component Dirac-Coulomb mean-field calculations ($^2$DC$^M$) is capable of providing core excitations and ionization energies that are nearly indistinguishable from the reference 4-component energies for up to and including fifth-row elements. We observe that two-electron integrals over the small-component basis sets yield non-negligible contributions to core binding energies for the K and L edges for atoms such as iodine or astatine, and that the approach based on Dirac-Coulomb-Gaunt mean-field calculations ($^2$DCG$^M$) are significantly more accurate than X2C calculations for which screened two-electron spin-orbit interactions are included via atomic mean-field integrals.
△ Less
Submitted 30 March, 2021; v1 submitted 17 November, 2020;
originally announced November 2020.
-
Semiflexible polymer in a gliding assay: reentrant transition, role of turnover and activity
Authors:
Amir Shee,
Nisha Gupta,
Abhishek Chaudhuri,
Debasish Chaudhuri
Abstract:
We consider a model of an extensible semiflexible filament moving in two dimensions on a motility assay of motor proteins represented explicitly as active harmonic linkers. Their heads bind stochastically to polymer segments within a capture radius, and extend along the filament in a directed fashion before detaching. Both the extension and detachment rates are load-dependent and generate an activ…
▽ More
We consider a model of an extensible semiflexible filament moving in two dimensions on a motility assay of motor proteins represented explicitly as active harmonic linkers. Their heads bind stochastically to polymer segments within a capture radius, and extend along the filament in a directed fashion before detaching. Both the extension and detachment rates are load-dependent and generate an active drive on the filament. The filament undergoes a first order phase transition from open chain to spiral conformations and shows a reentrant behavior in both the active extension and the turnover, defined as the ratio of attachment-detachment rates. Associated with the phase transition, the size and shape of the polymer changes non-monotonically, and the relevant autocorrelation functions display double-exponential decay. The corresponding correlation times show a maximum signifying the dominance of spirals. The orientational dynamics captures the rotation of spirals, and its correlation time decays with activity as a power law.
△ Less
Submitted 17 December, 2020; v1 submitted 20 June, 2020;
originally announced June 2020.
-
The DIRAC code for relativistic molecular calculations
Authors:
Trond Saue,
Radovan Bast,
Andre Severo Pereira Gomes,
Hans Jørgen Aagaard Jensen,
Lucas Visscher,
Ignacio Agustın Aucar,
Roberto Di Remigio,
Kenneth G. Dyall,
Ephraim Eliav,
Elke Faßhauer,
Timo Fleig,
Loıc Halbert,
Erik Donovan Hedegård,
Benjamin Helmich-Paris,
Miroslav Iliaš,
Christoph R. Jacob,
Stefan Knecht,
Jon K Laerdahl,
Marta L. Vidal,
Malaya K Nayak,
Małgorzata Olejniczak,
Jógvan Magnus Haugaard Olsen,
Markus Pernpointner,
Bruno Senjean,
Avijit Shee
, et al. (2 additional authors not shown)
Abstract:
DIRAC is a freely distributed general-purpose program system for 1-, 2- and 4-component relativistic molecular calculations at the level of Hartree--Fock, Kohn--Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, coupled cluster and electron propagator theory. At the self-consistent-field level a highly original scheme, bas…
▽ More
DIRAC is a freely distributed general-purpose program system for 1-, 2- and 4-component relativistic molecular calculations at the level of Hartree--Fock, Kohn--Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, coupled cluster and electron propagator theory. At the self-consistent-field level a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding, and frozen density embedding models. DIRAC was one of the earliest codes for relativistic molecular calculations and remains a reference in its field.
△ Less
Submitted 14 February, 2020;
originally announced February 2020.
-
Active Brownian particles: mapping to equilibrium polymers and exact computation of moments
Authors:
Amir Shee,
Abhishek Dhar,
Debasish Chaudhuri
Abstract:
It is well known that path probabilities of Brownian motion correspond to the equilibrium configurational probabilities of flexible Gaussian polymers, while those of active Brownian motion correspond to in-extensible semiflexible polymers. Here we investigate the properties of the equilibrium polymer that corresponds to the trajectories of particles acted on simultaneously by both Brownian as well…
▽ More
It is well known that path probabilities of Brownian motion correspond to the equilibrium configurational probabilities of flexible Gaussian polymers, while those of active Brownian motion correspond to in-extensible semiflexible polymers. Here we investigate the properties of the equilibrium polymer that corresponds to the trajectories of particles acted on simultaneously by both Brownian as well as active noise. Through this mapping we can see interesting crossovers in mechanical properties of the polymer with changing contour length. The polymer end-to-end distribution exhibits Gaussian behaviour for short lengths, which changes to the form of semiflexible filaments at intermediate lengths, to finally go back to a Gaussian form for long contour lengths. By performing a Laplace transform of the governing Fokker-Planck equation of the active Brownian particle, we discuss a direct method to derive exact expressions for all the moments of the relevant dynamical variables, in arbitrary dimensions. These are verified via numerical simulations and used to describe interesting qualitative features such as, for example, dynamical crossovers. Finally we discuss the kurtosis of the ABP's position which we compute exactly and show that it can be used to differentiate between active Brownian particles and active Ornstein-Uhlenbeck process.
△ Less
Submitted 11 December, 2020; v1 submitted 5 February, 2020;
originally announced February 2020.
-
Coupled Cluster as an impurity solver for Green's function embedding methods
Authors:
Avijit Shee,
Dominika Zgid
Abstract:
We investigate the performance of Green's function coupled cluster singles and doubles (CCSD) method as a solver for Green's function embedding methods. To develop an efficient CC solver, we construct the one-particle Green's function from the coupled cluster (CC) wave function based on a non-hermitian Lanczos algorithm. The major advantage of this method is that its scaling does not depend on the…
▽ More
We investigate the performance of Green's function coupled cluster singles and doubles (CCSD) method as a solver for Green's function embedding methods. To develop an efficient CC solver, we construct the one-particle Green's function from the coupled cluster (CC) wave function based on a non-hermitian Lanczos algorithm. The major advantage of this method is that its scaling does not depend on the number of frequency points. We have tested the applicability of the CC Green's function solver in the weakly to strongly correlated regimes by employing it for a half-filled 1D Hubbard model projected onto a single site impurity problem and a half-filled 2D Hubbard model projected onto a 4-site impurity problem. For the 1D Hubbard model, for all interaction strengths, we observe an excellent agreement with the full configuration interaction (FCI) technique, both for the self-energy and spectral function. For the 2D Hubbard, we have employed an open-shell version of the current implementation and observed some discrepancies from FCI in the strongly correlated regime. Finally, on an example of a small ammonia cluster, we analyze the performance of the Green's function CCSD solver within the self-energy embedding theory (SEET) with Hartee-Fock (HF) and Green's function second order (GF2) for the treatment of the environment.
△ Less
Submitted 10 June, 2019;
originally announced June 2019.
-
Towards Standardization of Data Licenses: The Montreal Data License
Authors:
Misha Benjamin,
Paul Gagnon,
Negar Rostamzadeh,
Chris Pal,
Yoshua Bengio,
Alex Shee
Abstract:
This paper provides a taxonomy for the licensing of data in the fields of artificial intelligence and machine learning. The paper's goal is to build towards a common framework for data licensing akin to the licensing of open source software. Increased transparency and resolving conceptual ambiguities in existing licensing language are two noted benefits of the approach proposed in the paper. In pa…
▽ More
This paper provides a taxonomy for the licensing of data in the fields of artificial intelligence and machine learning. The paper's goal is to build towards a common framework for data licensing akin to the licensing of open source software. Increased transparency and resolving conceptual ambiguities in existing licensing language are two noted benefits of the approach proposed in the paper. In parallel, such benefits may help foster fairer and more efficient markets for data through bringing about clearer tools and concepts that better define how data can be used in the fields of AI and ML. The paper's approach is summarized in a new family of data license language - \textit{the Montreal Data License (MDL)}. Alongside this new license, the authors and their collaborators have developed a web-based tool to generate license language espousing the taxonomies articulated in this paper.
△ Less
Submitted 20 March, 2019;
originally announced March 2019.
-
Predictive simulations of ionization energies of solvated halide ions with relativistic embedded Equation of Motion Coupled-Cluster Theory
Authors:
Yassine Bouchafra,
Avijit Shee,
Florent Real,
Valerie Vallet,
Andre Severo Pereira Gomes
Abstract:
A subsystem approach for obtaining electron binding energies in the valence region and apply it to the case of halide ions (X$^-$, X = F-At) in water. This approach is based on electronic structure calculations combining the relativistic equation of motion coupled-cluster method for electron detachment (EOM-IP-CCSD) and density functional theory via the frozen density embedding (FDE) approach, usi…
▽ More
A subsystem approach for obtaining electron binding energies in the valence region and apply it to the case of halide ions (X$^-$, X = F-At) in water. This approach is based on electronic structure calculations combining the relativistic equation of motion coupled-cluster method for electron detachment (EOM-IP-CCSD) and density functional theory via the frozen density embedding (FDE) approach, using structures from classical molecular dynamics with polarizable force fields for discrete systems (in the present study, droplets containing the anion and 50 water molecules). Our results indicate one can accurately capture both the large solvent effect observed for the halides as well as the splitting of their ionization signals due to the increasingly large spin-orbit coupling of the p$_{3/2}$-p$_{1/2}$ manifold across the series, at an affordable computational cost. Furthermore, due to the quantum mechanical treatment of both solute and solvent, electron binding energies of semi-quantitative quality are also obtained also obtained for (bulk) water as by-products of the calculations for the halogens (in droplets).
△ Less
Submitted 4 December, 2018; v1 submitted 5 November, 2018;
originally announced November 2018.
-
Equation-of-Motion Coupled-Cluster Theory based on the 4-component Dirac-Coulomb(-Gaunt) Hamiltonian. Energies for single electron detachment, attachment and electronically excited states
Authors:
Avijit Shee,
Trond Saue,
Lucas Visscher,
Andre Severo Pereira Gomes
Abstract:
We report in this paper an implementation of 4-component relativistic Hamiltonian based Equation-of-Motion Coupled-Cluster with singles and doubles (EOM-CCSD) theory for the calculation of ionization potential (IP), electron affinity (EA) and excitation energy (EE). In this work we utilize previously developed double group symmetry-based generalized tensor contraction scheme, and also extend it in…
▽ More
We report in this paper an implementation of 4-component relativistic Hamiltonian based Equation-of-Motion Coupled-Cluster with singles and doubles (EOM-CCSD) theory for the calculation of ionization potential (IP), electron affinity (EA) and excitation energy (EE). In this work we utilize previously developed double group symmetry-based generalized tensor contraction scheme, and also extend it in order to carry out tensor contractions involving non-totally symmetric and odd-ranked tensors. Several approximated spin-free and two-component Hamiltonians can also be accessed in this implementation. We have applied this method to the halogen monoxide (XO, X= Cl, Br, I, At, Ts) species, in order to assess the quality of a few other recent EOMCC implementations, where spin-orbit coupling contribution has been approximated in different degree. Besides, we also have studied various excited states of CH$_2$IBr, CH$_2$I$_2$ and I$_2^-$(as well as single electron attachment and detachment electronic states of the same species) where comparison has been made with a closely related multi-reference coupled-cluster method, namely Intermediate Hamiltonian Fock Space Coupled-Cluster singles and doubles (IHFS-CCSD) theory.
△ Less
Submitted 17 October, 2018; v1 submitted 24 August, 2018;
originally announced August 2018.
-
Testing self-energy embedding theory in combination with GW
Authors:
Tran Nguyen Lan,
Avijit Shee,
Jia Li,
Emanuel Gull,
Dominika Zgid
Abstract:
We present a theoretical framework and implementation details for self-energy embedding theory (SEET) with the GW approximation for the treatment of weakly correlated degrees of freedom and configuration interactions solver for handing the strongly correlated degrees. On a series of molecular examples, for which the exact results are known within a given basis, we demonstrate that SEET(CI/GW) is a…
▽ More
We present a theoretical framework and implementation details for self-energy embedding theory (SEET) with the GW approximation for the treatment of weakly correlated degrees of freedom and configuration interactions solver for handing the strongly correlated degrees. On a series of molecular examples, for which the exact results are known within a given basis, we demonstrate that SEET(CI/GW) is a systematically improvable and well controlled method capable of giving accurate results and well behaved causal self-energies and Green's functions. We compare the theoretical framework of SEET(CI/GW) to that of GW+DMFT and comment on differences between these to approaches that aim to treat both strongly and weakly correlated simultaneously.
△ Less
Submitted 18 June, 2017;
originally announced June 2017.