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

Michał Lesiuk

Michał Lesiuk contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phphysics.atom-phcond-mat.str-elphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
4topics
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

9 published item(s)

preprint2022arXiv

A systematic construction of Gaussian basis sets for the description of laser field ionization and high-harmonic generation

A precise understanding of mechanisms governing the dynamics of electrons in atoms and molecules subjected to intense laser fields has a key importance for the description of attosecond processes such as the high-harmonic generation and ionization. From the theoretical point of view, this is still a challenging task, as new approaches to solve the time-dependent Schrödinger equation with both good accuracy and efficiency are still emerging. Until recently, the purely numerical methods of real-time propagation of the wavefunction using finite grids have been frequently and successfully used to capture the electron dynamics in small one- or two-electron systems. However, as the main focus of attoscience shifts toward many-electron systems, such techniques are no longer effective and need to be replaced by more approximate but computationally efficient ones. In this paper, we explore the increasingly popular method of expanding the wavefunction of the examined system into a linear combination of atomic orbitals and present a novel systematic scheme for constructing an optimal Gaussian basis set suitable for the description of excited and continuum atomic or molecular states. We analyze the performance of the proposed basis sets by carrying out a series of time-dependent configuration interaction calculations for the hydrogen atom in fields of intensity varying from $5 \times 10^{13}\:\rm W/cm^2$ to $5 \times 10^{14}\:\rm W/cm^2$. We also compare the results with the data obtained using Gaussian basis sets proposed previously by other authors.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Explicitly correlated electronic structure calculations with transcorrelated matrix product operators

In this work, we present the first implementation of the transcorrelated electronic Hamiltonian in an optimization procedure for matrix product states by the density matrix renormalization group (DMRG) algorithm. In the transcorrelation ansatz, the electronic Hamiltonian is similarity-transformed with a Jastrow factor to describe the cusp in the wave function at electron-electron coalescence. As a result, the wave function is easier to approximate accurately with the conventional expansion in terms of one-particle basis functions and Slater determinants. The transcorrelated Hamiltonian in first quantization comprises up to three-body interactions, which we deal with in the standard way by applying robust density fitting to two- and three-body integrals entering the second-quantized representation of this Hamiltonian. The lack of hermiticity of the transcorrelated Hamiltonian is taken care of along the lines of the first work on transcorrelated DMRG [J. Chem. Phys. 153, 164115 (2020)] by encoding it as a matrix product operator and optimizing the corresponding ground state wave function with imaginary-time time-dependent DMRG. We demonstrate our quantum chemical transcorrelated DMRG approach at the example of several atoms and first-row diatomic molecules. We show that transcorrelation improves the convergence rate to the complete basis set limit in comparison to conventional DMRG. Moreover, we study extensions of our approach that aim at reducing the cost of handling the matrix product operator representation of the transcorrelated Hamiltonian.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Near-exact CCSDT energetics from rank-reduced formalism supplemented by non-iterative corrections

We introduce a non-iterative energy correction, added on top of the rank-reduced coupled-cluster method with single, double, and triple substitutions, that accounts for excitations excluded from the parent triple excitation subspace. The formula for the correction is derived by employing the coupled-cluster Lagrangian formalism with an additional assumption that the parent excitation subspace is closed under the action of the Fock operator. Owning to the rank-reduced form of the triple excitation amplitudes tensor, the computational cost of evaluating the correction scales as $N^7$ with the system size, $N$. The accuracy and computational efficiency of the proposed method is assessed both for total and relative correlation energies. We show that the non-iterative correction can fulfill two separate roles. If an accuracy level of a fraction of kJ/mol is sufficient for a given system the correction significantly reduces the dimension of the parent triple excitation subspace needed in the iterative part of the calculations. Simultaneously, it enables to reproduce the exact CCSDT results to an accuracy level below 0.1 kJ/mol with a larger, yet still reasonable, dimension of the parent excitation subspace. This typically can be achieved at a computational cost only several times larger than required for the CCSD(T) method. The proposed method retains black-box features of the single-reference coupled-cluster theory; the dimension of the parent excitation subspace remains the only additional parameter that has to be specified.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quintic-scaling rank-reduced coupled cluster theory with single and double excitations

We consider the rank-reduced coupled-cluster theory with single and double excitations (RR-CCSD) introduced recently [Parrish \emph{et al.}, J. Chem. Phys. {\bf 150}, 164118 (2019)]. The main feature of this method is the decomposed form of the doubly-excited amplitudes which are expanded in the basis of largest magnitude eigenvectors of the MP2 or MP3 amplitudes. This approach enables a substantial compression of the amplitudes with only minor loss of accuracy. However, the formal scaling of the computational costs with the system size ($N$) is unaffected in comparison with the conventional CCSD theory ($\propto N^6$) due to presence of some terms quadratic in the amplitudes. We show how to solve this problem, exploiting the fact that their effective rank increases only linearly with the system size and reduce the scaling of the RR-CCSD iterations down to the level of $N^5$. This is combined with an iterative method of finding dominant eigenpairs of the MP2 or MP3 amplitudes which eliminates the necessity to perform the complete diagonalization. Next, we consider the evaluation of the perturbative corrections to the CCSD energies resulting from triply excited configurations. The triply-excited amplitudes present in the CCSD(T) method are decomposed to the Tucker-3 format using the higher-order orthogonal iteration (HOOI) procedure. This enables to compute the energy correction due to triple excitations non-iteratively with $N^6$ cost. The accuracy of the resulting rank-reduced CCSD(T) method is studied both for total and relative correlation energies of a diverse set of molecules. Accuracy levels better than 99.9\% can be achieved with a substantial reduction of the computational costs. Concerning the computational timings, break-even point between the rank-reduced and conventional CCSD implementations occurs for systems with about $30-40$ active electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Potential energy curve for the $a^3Σ_u^+$ state of lithium dimer with Slater-type orbitals

We report state-of-the-art ab initio calculations of the potential energy curve for the $a^3Σ_u^+$ state of the lithium dimer conducted to achieve spectroscopic accuracy ($<$1cm$^{-1}$) without any prior adjustment to fit the corresponding experimental data. The nonrelativistic clamped-nuclei component of the interaction energy is calculated with a composite method involving six-electron coupled cluster and full configuration interaction theories combined with basis sets of Slater-type orbitals ranging in quality from double- to sextuple-zeta. To go beyond the nonrelativistic Born-Oppenheimer picture we include both the leading-order relativistic and adiabatic corrections, and find both of these effects to be non-negligible within the present accuracy standards. The potential energy curve developed by us allowed to calculate molecular parameters ($D_e$, $D_0$, $ω_e$ etc.) for this system, as well as the corresponding vibrational energy levels, with an error of only a few tenths of a wavenumber ($0.2-0.4\,$cm$^{-1}$). We also report an ab initio value for the scattering length of two $^2S$ lithium atoms which determines the stability of the related Bose-Einstein condensate.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Spin-orbit coupling matrix elements from the explicitly connected expressions of the response functions within the coupled-cluster theory

In this work we present a coupled cluster based approach to the computation of the spin orbit coupling matrix elements. The working expressions are derived from the quadratic response function with the coupled cluster parametrization, using the auxiliary excitation operator S. Systematic approximations are proposed with the CCSD and CC3 levels of theory. The new method is tested by computing lifetimes of several electronic states of Ca, Sr and Ba atoms, with Gaussian and Slater basis sets. The results are compared with available theoretical and experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Theoretical determination of polarizability and magnetic susceptibility of neon

We report theoretical determination of the dipole polarizability of the neon atom, including its frequency dependence. Corrections for the relativistic, quantum electrodynamics, finite nuclear mass, and finite nuclear size effects are taken into account. We obtain the value $α_0=2.66080(36)$ for the static polarizability, and $α_2=2.850(7)$ and $α_4=4.932(14)$ for the first two polarizability dispersion coefficients (Cauchy moments); all values are in atomic units (a.u.). In the case of static polarizability, our result agrees with the best experimental determination [C. Gaiser and B. Fellmuth, Phys. Rev. Lett. 120, 123203 (2018)], but our estimated uncertainty is significantly larger. For the dispersion coefficients, the results obtained in this work appear to be the most accurate to date overall compared to published theoretical and experimental data. We also calculated the static magnetic susceptibility of the neon atom, needed to obtain the refractive index of gaseous neon. Our result, $χ_0 = -8.484(19) \cdot 10^{-5}$ a.u., is about 9% larger in absolute value than the recommended experimental value [CRC Handbook of Chemistry and Physics, CRC Press, 2019, p. 4-145].

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A straightforward a posteriori method for reduction of density-fitting error in coupled-cluster calculations

We present a simple method for \emph{a posteriori} removal of a significant fraction of the density-fitting error from the calculated total coupled-cluster energies. The method treats the difference between the exact and density-fitted integrals as a perturbation, and simplified response-like equations allow to calculate improved amplitudes and the corresponding energy correction. The proposed method is tested at the coupled-cluster singles and doubles level of theory for a diverse set of moderately-sized molecules. On average, error reductions by a factor of approximately ten and twenty are observed in double-zeta and triple-zeta basis sets, respectively. Similar reductions are observed in calculations of interaction energies of several model complexes. The computational cost of the procedure is small in comparison with the preceding coupled-cluster iterations. The applicability of the method is not limited to the density-fitting approximation; in principle, it can be used in conjunction with an arbitrary decomposition scheme of the electron repulsion integrals.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Implementation of the full CCSDT electronic structure model with tensor decompositions

We report a complete implementation of the coupled-cluster method with single, double, and triple excitations (CCSDT) where tensor decompositions are used to reduce its scaling and overall computational costs. For the decomposition of the electron repulsion integrals the standard density fitting (or Cholesky decomposition) format is used. The coupled-cluster single and double amplitudes are treated conventionally, and for the triple amplitudes tensor we employ the Tucker-3 compression formula, $t_{ijk}^{abc} \approx t_{XYZ} \,U^X_{ai}\,U^Y_{bj} \,U^Z_{ck}$. The auxiliary quantities $U^X_{ai}$ come from singular value decomposition (SVD) of an approximate triple amplitudes tensor based on perturbation theory. The efficiency of the proposed method relies on an observation that the dimension of the &#34;compressed&#34; tensor $t_{XYZ}$ sufficient to deliver a constant relative accuracy of the correlation energy grows only linearly with the size of the system, $N$. This fact, combined with proper factorization of the coupled-cluster equations, leads to practically $N^6$ scaling of the computational costs of the proposed method, as illustrated numerically for linear alkanes with increasing chain length. This constitutes a considerable improvement over the $N^8$ scaling of the conventional (uncompressed) CCSDT theory. The accuracy of the proposed method is verified by benchmark calculations of total and relative energies for several small molecular systems and comparison with the exact CCSDT method. The accuracy levels of 1 kJ/mol are easily achievable with reasonable SVD subspace size, and even more demanding levels of accuracy can be reached with a considerable reduction of the computational costs. Extensions of the proposed method to include higher excitations are briefly discussed, along with possible strategies of reducing other residual errors.

0 citations0 reviewsNo code linkedOpen access