D. Van Neck

Read-Green resonances in a topological superconductor coupled to a bath

P. Claeys, S. De Baerdemacker, D. Van Neck
Physical Review B
93 (22) 220503(R)
2016
A1

Abstract 

We study a topological superconductor capable of exchanging particles with an environment. This additional interaction breaks particle-number symmetry and can be modeled by means of an integrable Hamiltonian, building on the class of Richardson-Gaudin pairing models. The isolated system supports zero-energy modes at a topological phase transition, which disappear when allowing for particle exchange with an environment. However, it is shown from the exact solution that these still play an important role in system-environment particle exchanges, which can be observed through resonances in low-energy and low-momentum level occupations. These fluctuations signal topologically protected Read-Green points and cannot be observed within traditional mean-field theory.

DMRG-CASPT2 study of the longitudinal static second hyperpolarizability of all-trans polyenes

S. Wouters, V. Van Speybroeck, D. Van Neck
Journal of Chemical Physics
145 (5), 054120
2016
A1

Abstract 

We have implemented internally contracted complete active space second order perturbation theory (CASPT2) with the density matrix renormalization group (DMRG) as active space solver [Y. Kurashige and T. Yanai, J. Chem. Phys. 135, 094104 (2011)]. Internally contracted CASPT2 requires to contract the generalized Fock matrix with the 4-particle reduced density matrix (4-RDM) of the reference wavefunction. The required 4-RDM elements can be obtained from 3-particle reduced density matrices (3-RDM) of different wavefunctions, formed by symmetry-conserving single-particle excitations op top of the reference wavefunction. In our spin-adapted DMRG code chemps2 [https://github.com/sebwouters/chemps2], we decompose these excited wavefunctions as spin-adapted matrix product states, and calculate their 3-RDM in order to obtain the required contraction of the generalized Fock matrix with the 4-RDM of the reference wavefunction. In this work, we study the longitudinal static second hyperpolarizability of all-trans polyenes C_{2n}H_{2n+2} [n = 4 - 12] in the cc-pVDZ basis set. DMRG-SCF and DMRG-CASPT2 yield substantially lower values and scaling with system size compared to RHF and MP2, respectively.

Maximum probability domains for Hubbard models

G. Acke, S. De Baerdemacker, P. Claeys, M. Van Raemdonck, W. Poelmans, D. Van Neck, P. Bultinck
Molecular Physics
114 (7-8), 1392-1405
2016
A1

Abstract 

The theory of maximum probability domains (MPDs) is formulated for the Hubbard model in terms of projection operators and generating functions for both exact eigenstates as well as Slater determinants. A fast MPD analysis procedure is proposed, which is subsequently used to analyse numerical results for the Hubbard model. It is shown that the essential physics behind the considered Hubbard models can be exposed using MPDs. Furthermore, the MPDs appear to be in line with what is expected from Valence Bond (VB) Theory-based knowledge.

Open Access version available at UGent repository

Glide of dislocations in < 1 1 1 >{321} slip system: an atomistic study

D. Terentyev, A. Bakaev, D. Van Neck, E. Zhurkin
Philosophical Magazine
96 (1), 71-83
2016
A1

Abstract 

Atomistic calculations are performed to investigate plastic slip in the {321} system in body-centred cubic iron. Several modern interatomic potentials, developed over the last decade, are applied to compute the stacking fault -line energy in the {321} plane and the results are compared with the ab initio prediction. The applied potentials have shown strong deviations, but several potentials acquired good qualitative agreement with the ab initio data. Depending on the applied potential, the lowest value of the Peierls stress for the edge dislocation (ED) is 50MPa (Ackland and Bacon from 1997) and the highest is 550MPa (Dudarev and Derlet from 2005), while for the screw dislocation it is much higher, in the range 1-2GPa. At finite temperature, however, the flow stress of the ED is found to decrease exponentially reaching a negligible value at about 200K, irrespective of the applied potential. On the basis of the data obtained using Ackland-Mendelev potential from 2004, we conclude that the slip resistance of the {321} system is in between the resistance of the {110} and {112} slip systems.

Eigenvalue-based determinants for scalar products and form factors in Richardson–Gaudin integrable models coupled to a bosonic mode

P. Claeys, S. De Baerdemacker, M. Van Raemdonck, D. Van Neck
Journal of Physics A: Mathematical and Theoretical
48 (42), 425201
2015
A1

Abstract 

Starting from integrable su(2) (quasi-)spin Richardson–Gaudin (RG) XXZ models we derive several properties of integrable spin models coupled to a bosonic mode. We focus on the Dicke–Jaynes–Cummings–Gaudin models and the two-channel (p + ip)-wave pairing Hamiltonian. The pseudo-deformation of the underlying su(2) algebra is here introduced as a way to obtain these models in the contraction limit of different RG models. This allows for the construction of the full set of conserved charges, the Bethe ansatz state, and the resulting RG equations. For these models an alternative and simpler set of quadratic equations can be found in terms of the eigenvalues of the conserved charges. Furthermore, the recently proposed eigenvalue-based determinant expressions for the overlaps and form factors of local operators are extended to these models, linking the results previously presented for the Dicke–Jaynes–Cummings–Gaudin models with the general results for RG XXZ models.

This paper has been selected as editor's pick

Polynomial scaling approximations and Dynamic Correlation Corrections to Doubly Occupied Configuration Interaction wave functions

M. Van Raemdonck, D. Alcoba, W. Poelmans, S. De Baerdemacker, A. Torre, L. Lain, G. Massaccesi, D. Van Neck, P. Bultinck
Journal of Chemical Physics
143 (10), 104106
2015
A1

Abstract 

A class of polynomial scaling methods that approximate Doubly Occupied Configuration Interaction (DOCI) wave functions and improve the description of dynamic correlation is introduced. The accuracy of the resulting wave functions is analysed by comparing energies and studying the overlap between the newly developed methods and full configuration interaction wave functions, showing that a low energy does not necessarily entail a good approximation of the exact wave function. Due to the dependence of DOCI wave functions on the single-particle basis chosen, several orbital optimisation algorithms are introduced. An energy-based algorithm using the simulated annealing method is used as a benchmark. As a computationally more affordable alternative, a seniority number minimising algorithm is developed and compared to the energy based one revealing that the seniority minimising orbital set performs well. Given a well-chosen orbital basis, it is shown that the newly developed DOCI based wave functions are especially suitable for the computationally efficient description of static correlation and to lesser extent dynamic correlation.

Open Access version available at UGent repository

Variational optimization of the second order density matrix corresponding to a seniority-zero configuration interaction wave function

W. Poelmans, M. Van Raemdonck, B. Verstichel, S. De Baerdemacker, A. Torre, L. Lain, G. Massaccesi, D. Alcoba, P. Bultinck, D. Van Neck
Journal of Chemical Theory and Computation (JCTC)
11 (9), 4064–4076
2015
A1

Abstract 

We perform a direct variational determination of the second-order (two-particle) density matrix corresponding to a many-electron system, under a restricted set of the two-index $N$-representability $\mathcal{P}$-, $\mathcal{Q}$-, and $\mathcal{G}$-conditions. In addition, we impose a set of necessary constraints that the two-particle density matrix must be derivable from a doubly-occupied many-electron wave function, i.e.\ a singlet wave function for which the Slater determinant decomposition only contains determinants in which spatial orbitals are doubly occupied. We rederive the two-index $N$-representability conditions first found by Weinhold and Wilson and apply them to various benchmark systems (linear hydrogen chains, He, $\text{N}_2$ and $\text{CN}^-$). This work is motivated by the fact that a doubly-occupied many-electron wave function captures in many cases the bulk of the static correlation. Compared to the general case, the structure of doubly-occupied two-particle density matrices causes the associate semidefinite program to have a very favorable scaling as $L^3$, where $L$ is the number of spatial orbitals. Since the doubly-occupied Hilbert space depends on the choice of the orbitals, variational calculation steps of the two-particle density matrix are interspersed with orbital-optimization steps (based on Jacobi rotations in the space of the spatial orbitals). We also point to the importance of symmetry breaking of the orbitals when performing calculations in a doubly-occupied framework.

The Dicke model as the contraction limit of a pseudo-deformed Richardson-Gaudin model

P. Claeys, S. De Baerdemacker, M. Van Raemdonck, D. Van Neck
Journal of Physics Conference Series
597, UNSP 012025
2015
P1

Abstract 

The Dicke model is derived in the contraction limit of a pseudo-deformation of the quasispin algebra in the su(2)-based Richardson-Gaudin models. Likewise, the integrability of the Dicke model is established by constructing the full set of conserved charges, the form of the Bethe Ansatz state, and the associated Richardson-Gaudin equations. Thanks to the formulation in terms of the pseudo-deformation, the connection from the su(2)-based Richardson-Gaudin model towards the Dicke model can be performed adiabatically.

Open Access version available at UGent repository

Normal mode analysis of macromolecular systems with the mobile block Hessian method

A. Ghysels, V. Van Speybroeck, D. Van Neck, B.R. Brooks, M. Waroquier
AIP Conference Proceedings
1642 (2015), 559
2015
P1

Abstract 

Until recently, normal mode analysis (NMA) was limited to small proteins, not only because the required energy minimization is a computationally exhausting task, but also because NMA requires the expensive diagonalization of a 3N(a) x 3N(a) matrix with N-a the number of atoms. A series of simplified models has been proposed, in particular the Rotation-Translation Blocks (RTB) method by Tama et al. for the simulation of proteins. It makes use of the concept that a peptide chain or protein can be seen as a subsequent set of rigid components, i.e. the peptide units. A peptide chain is thus divided into rigid blocks with six degrees of freedom each.

Recently we developed the Mobile Block Hessian (MBH) method, which in a sense has similar features as the RTB method. The main difference is that MBH was developed to deal with partially optimized systems. The position/orientation of each block is optimized while the internal geometry is kept fixed at a plausible - but not necessarily optimized - geometry. This reduces the computational cost of the energy minimization. Applying the standard NMA on a partially optimized structure however results in spurious imaginary frequencies and unwanted coordinate dependence. The MBH avoids these unphysical effects by taking into account energy gradient corrections. Moreover the number of variables is reduced, which facilitates the diagonalization of the Hessian.

In the original implementation of MBH, atoms could only be part of one rigid block. The MBH is now extended to the case where atoms can be part of two or more blocks. Two basic linkages can be realized: (1) blocks connected by one link atom, or (2) by two link atoms, where the latter is referred to as the hinge type connection. In this work we present the MBH concept and illustrate its performance with the crambin protein as an example.

PPV Polymerization via the Gilch Route: Diradical Character of Monomers

J.D. Nikolic, S. Wouters, J. Romanova, A. Shimizu, B. Champagne, T. Junkers, D. Vanderzande, D. Van Neck, M. Waroquier, V. Van Speybroeck, S. Catak
Chemistry - A European Journal
21, 19176-19185
2015
A1

Abstract 

Despite various studies on the polymerization of poly(p-phenylene vinylene) (PPV) through different precursor routes, detailed mechanistic knowledge on the individual reaction steps and intermediates is still incomplete. The present study aims to gain more insight into the radical polymerization of PPV through the Gilch route. The initial steps of the polymerization involve the formation of a p-quinodimethane intermediate, which spontaneously self-initiates through a dimerization process leading to the formation of diradical species; chain propagation ensues on both sides of the diradical or chain termination occurs by the formation of side products, such as [2.2]paracyclophanes. Furthermore, different p-quinodimethane systems were assessed with respect to the size of their aromatic core as well as the presence of heteroatoms in/on the conjugated system. The nature of the aromatic core and the specific substituents alter the electronic structure of the p-quinodimethane monomers, affecting the mechanism of polymerization. The diradical character of the monomers has been investigated with several advanced methodologies, such as spin-projected UHF, CASSCF, CASPT2, and DMRG calculations. It was shown that larger aromatic cores led to a higher diradical character in the monomers, which in turn is proposed to cause rapid initiation.

Pages

Subscribe to RSS - D. Van Neck