D. Van Neck

CheMPS2: Improved DMRG-SCF routine and correlation functions

S. Wouters, W. Poelmans, S. De Baerdemacker, P.W. Ayers, D. Van Neck
Computer Physics Communications
191, 235-237
2015
A1

Abstract 

CheMPS2, our spin-adapted implementation of the density matrix renormalization group (DMRG) for ab initio quantum chemistry (Wouters et al., 2014), has several new features. A speed-up of the augmented Hessian Newton–Raphson DMRG self-consistent field (DMRG-SCF) routine is achieved with the direct inversion of the iterative subspace (DIIS). For extended molecules, the active space orbitals can be localized by maximizing the Edmiston–Ruedenberg cost function. These localized orbitals can be ordered according to the topology of the molecule by approximately minimizing the bandwidth of the exchange matrix with the Fiedler vector. The electronic structure can be analyzed by means of the two-orbital mutual information, spin, spin-flip, density, and singlet diradical correlation functions.

Eigenvalue-based method and form-factor determinant representations for integrable XXZ Richardson-Gaudin models

P. Claeys, S. De Baerdemacker, M. Van Raemdonck, D. Van Neck
Physical Review B
91 (15), 155102
2015
A1

Abstract 

We propose an extension of the numerical approach for integrable Richardson-Gaudin models based on a new set of eigenvalue-based variables [A. Faribault et al., Phys. Rev. B 83, 235124 (2011); O. El Araby et al., Phys. Rev. B 85, 115130 (2012)]. Starting solely from the Gaudin algebra, the approach is generalized towards the full class of XXZ Richardson-Gaudin models. This allows for a fast and robust numerical determination of the spectral properties of these models, avoiding the singularities usually arising at the so-called singular points. We also provide different determinant expressions for the normalization of the Bethe ansatz states and form factors of local spin operators, opening up possibilities for the study of larger systems, both integrable and nonintegrable. These expressions can be written in terms of the new set of variables and generalize the results previously obtained for rational Richardson-Gaudin models [A. Faribault and D. Schuricht, J. Phys. A 45, 485202 (2012)] and Dicke-Jaynes-Cummings-Gaudin models [H. Tschirhart and A. Faribault,  J. Phys. A 47, 405204 (2014)]. Remarkably, these results are independent of the explicit parametrization of the Gaudin algebra, exposing a universality in the properties of Richardson-Gaudin integrable systems deeply linked to the underlying algebraic structure.

Open Access version available at UGent repository

A hybrid configuration interaction treatment based on seniority number and excitation schemes

D.R. Alcoba, A. Torre, L. Lain, O.B. Ona, P. Capuzzi, M. Van Raemdonck, P. Bultinck, D. Van Neck
Journal of Chemical Physics
141 (24), 244118
2014
A1

Abstract 

We present a configuration interaction method in which the Hamiltonian of an N-electron system is projected on Slater determinants selected according to the seniority-number criterion along with the traditional excitation-based procedure. This proposed method is especially useful to describe systems which exhibit dynamic (weak) correlation at determined geometric arrangements (where the excitation-based procedure is more suitable) but show static (strong) correlation at other arrangements (where the seniority-number technique is preferred). The hybrid method amends the shortcomings of both individual determinant selection procedures, yielding correct shapes of potential energy curves with results closer to those provided by the full configuration interaction method. (c) 2014 AIP Publishing LLC.

Open Access version available at UGent repository

Probing pairing correlations in Sn isotopes using Richardson-Gaudin integrability

S. De Baerdemacker, V. Hellemans, R. van den Berg, J.-S. Caux, K. Heyde, M. Van Raemdonck, D. Van Neck, P.A. Johnson, A. Buekenhoudt
Journal of Physics: Conference series
533, 012058
2014
P1

Abstract 

Pairing correlations in the even-even A = 102 − 130 Sn isotopes are discussed, based on the Richardson-Gaudin variables in an exact Woods-Saxon plus reduced BCS pairing framework. The integrability of the model sheds light on the pairing correlations, in particular on the previously reported sub-shell structure.

Open Access version available at UGent repository

Non-Variational Orbital Optimization Techniques for the AP1roG Wave Function

K. Boguslawski, P. Tecmer, P.W. Ayers, P. Bultinck, S. De Baerdemacker, D. Van Neck
Journal of Chemical Theory and Computation (JCTC)
10 (11), 4873-4882
2014
A1

Abstract 

We introduce new nonvariational orbital optimization schemes for the antisymmetric product of one-reference orbital geminal (AP1roG) wave function (also known as pair-coupled cluster doubles) that are extensions to our recently proposed projected seniority-two (PS2-AP1roG) orbital optimization method [ J. Chem. Phys. 2014, 140, 214114)]. These approaches represent less stringent approximations to the PS2-AP1roG ansatz and prove to be more robust approximations to the variational orbital optimization scheme than PS2-AP1roG. The performance of the proposed orbital optimization techniques is illustrated for a number of well-known multireference problems: the insertion of Be into H2, the automerization process of cyclobutadiene, the stability of the monocyclic form of pyridyne, and the aromatic stability of benzene.

Efficient description of strongly correlated electrons with mean-field cost

K. Boguslawski, P. Tecmer, P.W. Ayers, P. Bultinck, S. De Baerdemacker, D. Van Neck
Physical Review B
89 (20), 201106
2014
A1

Abstract 

We present an efficient approach to the electron correlation problem that is well suited for strongly interacting many-body systems, but requires only mean-field-like computational cost. The performance of our approach is illustrated for one-dimensional Hubbard rings with different numbers of sites, and for the nonrelativistic quantum-chemical Hamiltonian exploring the symmetric dissociation of the H-50 hydrogen chain.

Open Access version available at UGent repository

Interaction of carbon-vacancy complex with minor alloying elements of ferritic steels

A. Bakaev, D. Terentyev, X. He, E. Zhurkin, D. Van Neck
Journal of Nuclear Materials
451 (1-3), 82-87
2014
A1

Abstract 

Interstitial carbon, dissolved in bcc matrix of ferritic steels, plays an important role in the evolution of radiation-induced microstructure since it exhibits strong interaction with vacancies. Frequent formation and break-up of carbon-vacancy pairs, occurring in the course of irradiation, affect both kinetics of the accumulation of point defect clusters and carbon spatial distribution. The interaction of typical alloying elements (Mn, Ni, Cu, Si, Cr and P) in ferritic steels used as structural materials in nuclear reactors with a carbon-vacancy complex is analyzed using ab initio techniques. It is found that all the considered solutes form stable triple clusters resulting in the increase of the total binding energy by 0.2-0.3 eV. As a result of the formation of energetically favourable solute-carbon-vacancy triplets, the dissociation energy for vacancy/carbon emission is also increased by similar to 0.2-0.3 eV, suggesting that the solutes enhance thermal stability of carbon-vacancy complex. Association of carbon-vacancy pairs with multiple solute clusters is found to be favorable for Ni, Cu and P. The energetic stability of solute(s)-carbon-vacancy complexes was rationalized on the basis of pairwise interaction data and by analyzing the variation of local magnetic moments on atoms constituting the clusters. (C) 2014 Elsevier B.V. All rights reserved.

DOI 

10.1016/j.jnucmat.2014.03.031

Dislocation mechanism of deuterium retention in tungsten under plasma implantation

V. Dubinko, D. Terentyev, P. Grigorev, A. Bakaev, G. Van Oost, F. Gao, D. Van Neck, E. Zhurkin
Journal of Physics: Condensed Matter
26, 395001
2014
A1

Abstract 

We have developed a new theoretical model for deuterium (D) retention in tungsten-based alloys on the basis of its being trapped at dislocations and transported to the surface via the dislocation network with parameters determined by ab initio calculations. The model is used to explain experimentally observed trends of D retention under sub-threshold implantation, which does not produce stable lattice defects to act as traps for D in conventional models. Saturation of D retention with implantation dose and effects due to alloying of tungsten with, e.g. tantalum, are evaluated, and comparison of the model predictions with experimental observations under high-flux plasma implantation conditions is presented.

Energetic stability of solute-carbon-vacancy complexes in bcc iron

A. Bakaev, D. Terentyev, E. Zhurkin, D. Van Neck
Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms
352, 47-50
2015
A1

Abstract 

The strong binding between a vacancy and carbon in bcc iron plays an important role in the evolution of radiation-induced microstructure. Our previous ab initio study points to the fact that the vacancy-carbon (V-C) pair can serve as a nucleus for the solute-rich clusters. Here, we continue the ab initio study by considering the interaction of mixed solute clusters (Mn, Ni and Si) with the V-C pair, and the interaction of typical alloying elements of Fe-based steels (i.e., Mn, Ni, Cu, Si, Cr and P) with di-carbon-vacancy pair (V-C-2). We have identified the sequence of growth of Ni, Si and Mn solute-rich clusters nucleating on the V-C pair. The mixed-solute-V-C configurations are found to be less stable clusters than pure-solute-V-C clusters with the energy difference up to 0.22 eV per four atoms. The V-C-2 pair is found to be as strong nucleation site for the solute-rich clusters as the V-C pair. Only Si solute atom stands out from the trend showing a weaker affinity to the V-C-2 complex by 0.09 eV compared to the attraction to the V-C pair. The overall results point to the importance of taking into account the existence of both V-C and V-C-2 complexes in studying the formation of solute-rich clusters in Fe-based steels for nuclear applications. (C) 2014 Elsevier B.V. All rights reserved.

Review of many-body central force potentials for tungsten

G. Bonny, D. Terentyev, A. Bakaev, P. Grigorev, D. Van Neck
Modelling and Simulation in Materials Science and Engineering
22, 053001
2014
A1

Abstract 

Tungsten and tungsten-based alloys are the primary candidate materials for plasma facing components in fusion reactors. The exposure to high-energy radiation, however, severely degrades the performance and lifetime limits of the in-vessel components. In an effort to better understand the mechanisms driving the materials' degradation at the atomic level, large-scale atomistic simulations are performed to complement experimental investigations. At the core of such simulations lies the interatomic potential, on which all subsequent results hinge. In this work we review 19 central force many-body potentials and benchmark their performance against experiments and density functional theory (DFT) calculations. As basic features we consider the relative lattice stability, elastic constants and point-defect properties. In addition, we also investigate extended lattice defects, namely: free surfaces, symmetric tilt grain boundaries, the 1/2〈1 1 1〉{1 1 0} and 1/2〈1 1 1〉 {1 1 2} stacking fault energy profiles and the 1/2〈1 1 1〉 screw dislocation core. We also provide the Peierls stress for the 1/2〈1 1 1〉 edge and screw dislocations as well as the glide path of the latter at zero Kelvin. The presented results serve as an initial guide and reference list for both the modelling of atomically-driven phenomena in bcc tungsten, and the further development of its potentials.

Pages

Subscribe to RSS - D. Van Neck