BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Hannes Jónsson

Hannes Jónsson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.chem-phcond-mat.mes-hallcond-mat.mtrl-sciphysics.comp-phcond-mat.softBiological Physicscond-mat.stat-mechphysics.atm-clusquant-ph

Catalog footprint

What is connected

19works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

Indirect Mechanism of Au adatom Diffusion on the Si(100) Surface

Calculations of the diffusion of a Au adatom on the dimer reconstructed Si(100)-2x1 surface reveal an interesting mechanism that differs significantly from a direct path between optimal binding sites, which are located in between dimer rows. Instead, the active diffusion mechanism involves promotion of the adatom to higher energy sites on top of a dimer row and then fast migration along the row, visiting ca. a hundred sites at room temperature, before falling back down into an optimal binding site. This top-of-row mechanism becomes more important the lower the temperature is. The calculations are carried out by finding minimum energy paths on the energy surface obtained from density functional theory within the PBEsol functional approximation followed by kinetic Monte Carlo simulations of the diffusion over a range of temperature from 200 K to 900 K. While the activation energy for the direct diffusion mechanism is calculated to be 0.84 eV, the effective activation energy for the indirect mechanism is on average 0.56 eV.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reassignment of magic numbers for icosahedral Au clusters: 310, 564, 928 and 1426

Icosahedral Au clusters with three and four shells of atoms are found to deviate significantly from the commonly assumed Mackay structures. By introducing additional atoms in the surface shell and creating a vacancy in the center of the cluster, the calculated energy per atom can be lowered significantly, according to several different descriptions of the interatomic interaction. Analogous icosahedral structures with five and six shells of atoms are generated using the same structural motifs and are similarly found to be more stable than Mackay icosahedra. The lowest energy per atom is obtained with clusters containing 310, 564, 928 and 1426 atoms, as compared with the commonly assumed magic numbers of 309, 561, 923 and 1415. Some of the vertices in the optimized clusters have a hexagonal ring of atoms, rather than a pentagon, with the vertex atom missing. An inner shell atom in some cases moves outwards by more than an Ångström into the surface shell at the vertex site. This feature, as well as the wide distribution of nearest-neighbor distances in the surface layer, can strongly influence the catalytic properties of icosahedral clusters. The structural optimization is initially carried out using the GOUST method with atomic forces estimated with the EMT empirical potential function, but the atomic coordinates are then refined by minimization using electron density functional theory (DFT) or Gaussian approximation potential (GAP). A single energy barrier is found to separate the Mackay icosahedron from a lower energy structure where a string of atoms moves outwards in a concerted manner from the center so as to create a central vacancy while placing an additional atom in the surface shell.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Simulations of the Electrochemical Oxidation of Pt Nanoparticles of Various Shapes

The activity and stability of a platinum nanoparticle (NP) is not only affected by its size but additionally depends on its shape. To this end, simulations can identify structure-property relationships to make a priori decisions on the most promising structures. While activity is routinely probed by electronic structure calculations on simplified surface models, modeling the stability of NP model systems in electrochemical reactions is challenging due to the long timescale of relevant processes such as oxidation beyond the point of reversibility. In this work, a routine for simulating electrocatalyst stability is presented. The procedure is referred to as GREG after its main ingredients - a grand-canonical simulation approach using reactive force fields to model electrochemical reactions as a function of the galvanic cell potential. The GREG routine is applied to study the oxidation of 3 nm octahedral, cubic, dodecahedral, cuboctahedral, spherical, and tetrahexahedral platinum NPs. The oxidation process is analyzed using adsorption isobars as well as interaction energy heat maps that provide the basis for constructing electrochemical phase diagrams. Onset potentials for surface oxidation increase in the sequence cube ~= dodecahedron <= octahedron <= tetrahexahdron < sphere < cuboctahedron, establishing a relationship between oxidation behavior and surface facet structure. The electrochemical results are rationalized using structural and electronic analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Transferable Potential Function for Flexible H$_2$O Molecules Based on the Single Center Multipole Expansion

A potential function is presented for describing a system of flexible H$_2$O molecules based on the single center multipole expansion (SCME) of the electrostatic interaction. The model, referred to as SCME/f, includes the variation of the molecular quadrupole moment as well as the dipole moment with changes in bond length and angle so as to reproduce results of high level electronic structure calculations. The multipole expansion also includes fixed octupole and hexadecapole moments, as well as anisotropic dipole-dipole, dipole-quadrupole and quadrupole-quadrupole polarizability tensors. The model contains five adjustable parameters related to the repulsive interaction and damping functions in the electrostatic and dispersion interactions. Their values are adjusted to reproduce the lowest energy isomers of small clusters, (H$_2$O)$_n$ with $n=2-6$, as well as measured properties of the ice Ih crystal. Subsequent calculations of the energy difference between the various isomer configurations of the clusters show that SCME/f gives good agreement with results of electronic structure calculations and represents a significant improvement over the previously presented rigid SCME potential function. Analysis of the vibrational frequencies of the clusters and structural properties of ice Ih crystal show the importance of accurately describing the variation of the quadrupole moment with molecular structure.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Assessment of the Accuracy of Density Functionals for Calculating Oxygen Reduction Reaction on Nitrogen Doped Graphene

Experimental studies of the oxygen reduction reaction (ORR) at nitrogen doped graphene electrodes have reported a remarkably low overpotential, on the order of 0.5 V, similar to Pt based electrodes. Theoretical calculations using density functional theory have lent support for this claim. However, other measurements have indicated that transition metal impurities are actually responsible for the ORR activity, thereby raising questions about the reliability of both the experiments and the calculations. In order to assess the accuracy of the theoretical calculations, various generalized gradient approximation (GGA), meta-GGA and hybrid functionals are employed here and calibrated against high-level wave function based coupled cluster calculations (CCSD(T)) of the overpotential as well as self-interaction corrected density functional calculations and published quantum Monte Carlo calculations of O adatom binding to graphene. The PBE0 and HSE06 hybrid functionals are found to give more accurate results than the GGA and meta-GGA functionals, as would be expected, and for low dopant concentration, 3.1%, the overpotential is calculated to be 1.0 V. The GGA and meta-GGA functionals give a lower estimate by as much as 0.4 V. When the dopant concentration is doubled, the overpotential calculated with hybrid functionals drops, while it increases in GGA functional calculations. The opposite trends result from different potential determining steps, the *OOH species being of central importance in the hybrid functional calculations while the reduction of *O determines the overpotential obtained in GGA and meta-GGA calculations. The results presented here are mainly based on calculations of periodic representations of the system, but a comparison is also made with molecular flake models which are found to give erratic results.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Elastic Collision Based Dynamic Partitioning Scheme for Hybrid Simulations

The scattering-adapted flexible inner region ensemble separator (SAFIRES) is a partitioning scheme designed to divide a simulation cell into two regions to be treated with different computational methodologies. SAFIRES prevents particles from crossing between regions and resolves boundary events through elastic collisions of the particles mediated by the boundary, conserving energy and momenta. A multiple-time-step propagation algorithm is introduced where the time step is scaled automatically to identify the moment a collision occurs. If the length of the time step is kept constant, the new propagator reduces to a regular algorithm for Langevin dynamics, and to the velocity Verlet algorithm for classical dynamics if the friction coefficient is set to zero. SAFIRES constitutes the exact limit of the premise behind boundary-based methods such as FIRES, BEST, and BCC which take advantage of the indistinguishability of molecules on opposite sides of the separator. It gives correct average ensemble statistics despite the introduction of an ensemble separator. SAFIRES is tested in simulations where the molecules on the two sides are treated in the same way, for a Lennard-Jones (LJ) liquid and a LJ liquid in contact with a surface, as well as for liquid modelling simulations using the TIP4P force field. Simulations using SAFIRES are shown to reproduce the unconstrained reference simulations without significant deviations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetic Skyrmion Annihilation by Quantum Mechanical Tunneling

Magnetic skyrmions are nano-scale magnetic states that could be used in various spintronics devices. A central issue is the mechanism and rate of various possible annihilation processes and the lifetime of metastable skyrmions. While most studies have focused on classical over-the-barrier mechanism for annihilation, it is also possible that quantum mechanical tunneling through the energy barrier takes place. Calculations of the lifetime of magnetic skyrmions in a two-dimensional lattice are presented and the rate of tunneling compared with the classical annihilation rate. A remarkably strong variation in the crossover temperature and the lifetime of the skyrmion is found as a function of the values of parameters in the extended Heisenberg Hamiltonian, i.e. the out-of-plane anisotropy, Dzyaloshinskii-Moriya interaction (DMI) and applied magnetic field. Materials parameters and conditions are identified where the onset of tunneling could be observed on a laboratory time scale. In particular, it is predicted that skyrmion tunneling could be observed in the PdFe/Ir(111) system when an external magnetic field on the order of 6 T is applied.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Skyrmions in antiferromagnets: thermal stability and the effect of external field and impurities

Calculations of skyrmions in antiferromagnets (AFMs) are presented, and their properties compared with skyrmions in corresponding ferromagnets (FMs). The rates of skyrmion collapse and escape through the boundary of a track, as well as the binding to and collapse at a non-magnetic impurity, are calculated as a function of applied magnetic field. The activation energy for skyrmion annihilation is the same in AFMs and corresponding FMs in the absence of an applied magnetic field. The pre-exponential factor in the Arrhenius rate law is, however, different because skyrmion dynamics is different in the two systems. An applied magnetic field has opposite effects on skyrmions in the two types of materials. In AFMs the rate of collapse of skyrmions as well as the rate of escape through the edge of a magnetic strip decreases slightly with increasing field, while these rates increase strongly for a skyrmion in the corresponding FMs when the field is directed antiparallel to the magnetization in the center of the skyrmion. A non-magnetic impurity is less likely to trap a skyrmion in AFMs especially in the presence of a magnetic field. This, together with the established fact that a spin polarized current moves skyrmions in AFMs in the direction of the current, while in FMs skyrmions move at an angle to the current, demonstrates that skyrmions in AFMs have several advantageous properties over skyrmions in FMs for memory and spintronic devices.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Efficient dynamical correction of the transition state theory rate estimate for a flat energy barrier

The recrossing correction to the transition state theory estimate of a thermal rate can be difficult to calculate when the energy barrier is flat. This problem arises, for example, in polymer escape if the polymer is long enough to stretch between the initial and final state energy wells while the polymer beads undergo diffusive motion back and forth over the barrier. We present an efficient method for evaluating the correction factor by constructing a sequence of hyperplanes starting at the transition state and calculating the probability that the system advances from one hyperplane to another towards the product. This is analogous to what is done in forward flux sampling except that there the hyperplane sequence starts at the initial state. The method is applied to the escape of polymers with up to 64 beads from a potential well. For high temperature, the results are compared with direct Langevin dynamics simulations as well as forward flux sampling and excellent agreement between the three rate estimates is found. The use of a sequence of hyperplanes in the evaluation of the recrossing correction speeds up the calculation by an order of magnitude as compared with the traditional approach. As the temperature is lowered, the direct Langevin dynamics simulations as well as the forward flux simulations become computationally too demanding, while the harmonic transition state theory estimate corrected for recrossings can be calculated without significant increase in the computational effort.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Long-timescale simulations of H$_2$O admolecule diffusion on Ice Ih(0001) surfaces

Long-timescale simulations of the diffusion of a H$_2$O admolecule on the (0001) basal plane of ice Ih were carried out over a temperature range of 100 to 200 K using the adaptive kinetic Monte Carlo method and TIP4P/2005f interaction potential function. The arrangement of dangling H atoms was varied from the proton-disordered surface to the perfectly ordered Fletcher surface. A large variety of sites was found leading to a broad distribution in adsorption energy at both types of surfaces. Up to 4 % of the sites on the proton-disordered surface have an adsorption energy exceeding the cohesive energy of ice Ih. The mean squared displacement of a simulated trajectory at 175 K for the proton-disordered surface gave a diffusion constant of 6$\cdot$10$^{-10}$ cm$^2$/s, consistent with an upper bound previously reported from experimental measurements. During the simulation, dangling H atoms were found to rearrange so as to reduce clustering, thereby approaching a linear Fletcher type arrangement. Diffusion on the perfectly ordered Fletcher surface was estimated to be significantly faster, especially in the direction along the rows of dangling hydrogen atoms. From simulations over the range in temperature, an effective activation energy of diffusion was estimated to be 0.16 eV and 0.22 eV for diffusion parallel and perpendicular to the rows, respectively. Even a slight disruption of the rows of the Fletcher surface made the diffusion isotropic.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Method for finding mechanism and activation energy of magnetic transitions, applied to skyrmion and antivortex annihilation

A method for finding minimum energy paths of transitions in magnetic systems is presented and used to determine the mechanism and estimate the activation energy of skyrmion and antivortex annihiliation in nano-systems. The path is optimized with respect to orientation of the magnetic vectors while their magnitudes are fixed or obtained from separate calculations. The curvature of the configuration space is taken into account by: (1) using geodesics to evaluate distances and displacements of the system during the optimization, and (2) projecting the path tangent and the magnetic force on the tangent space of the manifold defined by all possible orientations of the magnetic vectors. The method, named geodesic nudged elastic band (GNEB), and its implementation are illustrated with calculations of complex transitions involving annihilation and creation of skyrmion and antivortex states. The lifetime of the latter was determined within harmonic transition state theory using a noncollinear extension of the Alexander-Anderson model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Polymer escape from a confining potential

The rate of escape of polymers from a two-dimensionally confining potential well has been evaluated using self-avoiding as well as ideal chain representations of varying length, up to 80 beads. Long timescale Langevin trajectories were calculated using the path integral hyperdynamics method to evaluate the escape rate. A minimum is found in the rate for self-avoiding polymers of intermediate length while the escape rate decreases monotonically with polymer length for ideal polymers. The increase in the rate for long, self-avoiding polymers is ascribed to crowding in the potential well which reduces the free energy escape barrier. An effective potential curve obtained using the centroid as an independent variable was evaluated by thermodynamic averaging and Kramers rate theory then applied to estimate the escape rate. While the qualitative features are well reproduced by this approach, it significantly overestimates the rate, especially for the longer polymers. The reason for this is illustrated by constructing a two-dimensional effective energy surface using the radius of gyration as well as the centroid as controlled variables. This shows that the description of a transition state dividing surface using only the centroid fails to confine the system to the region corresponding to the free energy barrier and this problem becomes more pronounced the longer the polymer is. A proper definition of a transition state for polymer escape needs to take into account the shape as well as the location of the polymer.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Qualitative Insight and Quantitative Analysis of the Effect of Temperature on the Coercivity of a Magnetic System

The temperature dependence of the response of a magnetic system to an applied field can be understood qualitatively by considering variations in the energy surface characterizing the system and estimated quantitatively with rate theory. In the system analysed here, Fe/Sm-Co spring magnet, the width of the hysteresis loop is reduced to a half when temperature is raised from 25~K to 300~K. This narrowing can be explained and reproduced quantitatively without invoking temperature dependence of model parameters as has typically been done in previous data analysis. The applied magnetic field lowers the energy barrier for reorientation of the magnetization but thermal activation brings the system over the barrier. A 2-dimensional representation of the energy surface is developed and used to gain insight into the transition mechanism and to demonstrate how the applied field alters the transition path. Our results show the importance of explicitly including the effect of thermal activation when interpreting experiments involving the manipulation of magnetic systems at finite temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Transition State Theory Approach to Polymer Escape from a One Dimensional Potential Well

The rate of escape of an ideal bead-spring polymer in a symmetric double-well potential is calculated using transition state theory (TST) and the results compared with direct dynamical simulations. The minimum energy path of the transitions becomes flat and the dynamics diffusive for long polymers making the Kramers-Langer estimate poor. However, TST with dynamical corrections based on short time trajectories started at the transition state gives rate constant estimates that agree within a factor of two with the molecular dynamics simulations over a wide range of bead coupling constants and polymer lengths. The computational effort required by the TST approach does not depend on the escape rate and is much smaller than that required by molecular dynamics simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Calculations of magnetic states and minimum energy paths of transitions using a noncollinear extension of the Alexander-Anderson model and a magnetic force theorem

Calculations of stable and metastable magnetic states as well as minimum energy paths for transitions between states are carried out using a noncollinear extension of the multiple-impurity Alexander-Anderson model and a magnetic force theorem which is derived and used to evaluate the total energy gradient with respect to orientation of magnetic moments -- an important tool for efficient navigation on the energy surface. By using this force theorem, the search for stable and metastable magnetic states as well as minimum energy paths revealing the mechanism and activation energy of transitions can be carried out efficiently. For Fe monolayer on W(110) surface, the model gives magnetic moment as well as exchange coupling between nearest and next-nearest neighbors that are in good agreement with previous density functional theory calculations. When applied to nanoscale Fe islands on this surface, the magnetic moment is predicted to be 10\% larger for atoms at the island rim, explaining in part an experimentally observed trend in the energy barrier for magnetization reversal in small islands. Surprisingly, the magnetic moment of the atoms does not change much along the minimum energy path for the transitions, which for islands containing more than 15 atom rows along either $[001]$ or $[1\bar{1}0]$ directions involves the formation of a thin, temporary domain wall. A noncollinear magnetic state is identified in a $7\times 7$ atomic row Fe island where the magnetic moments are arranged in an antivortex configuration with the central ones pointing out of the $(110)$ plane. This illustrates how the model can describe complicated exchange interactions even though it contains only a few parameters. The minimum energy path between this antivortex state and the collinear ground state is also calculated and the thermal stability of the antivortex state estimated.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Improved initial guess for minimum energy path calculations

A method is presented for generating a good initial guess of a transition path between given initial and final states of a system without evaluation of the energy. An objective function surface is constructed using an interpolation of pairwise distances at each discretization point along the path and the nudged elastic band method then used to find an optimal path on this image dependent pair potential (IDPP) surface. This provides an initial path for the more computationally intensive calculations of the true minimum energy path using some method of choice for evaluating the energy and atomic forces, for example by ab initio or density functional theory. The optimal path on the IDPP surface is significantly closer to the true minimum energy path than a linear interpolation of the Cartesian coordinates and, therefore, reduces the number of iterations needed to reach convergence and averts divergence in the electronic structure calculations when atoms are brought too close to each other in the initial path. The method is illustrated with three examples: (1) rotation of a methyl group in an ethane molecule, (2) an exchange of atoms in an island on a crystal surface, and (3) an exchange of two Si-atoms in amorphous silicon. In all three cases, the computational effort in finding the minimum energy path with DFT was reduced by a factor ranging from 50 % to an order of magnitude by using an IDPP path as the initial path. The time required for parallel computations was reduced even more because of load imbalance when linear interpolation of Cartesian coordinates was used.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Molecular reordering processes on ice (0001) surfaces from long timescale simulations

We report results of long timescale adaptive kinetic Monte Carlo simulations aimed at identifying possible molecular reordering processes on both proton-disordered and ordered (Fletcher) basal plane (0001) surfaces of hexagonal ice. The simulations are based on a force field for flexible molecules and span a time interval of up to 50 μs at a temperature of 100 K, which represents a lower bound to the temperature range of Earth's atmosphere. Additional calculations using both density functional theory and an ab initio based polarizable potential function are performed to test and refine the force field predictions. Several distinct processes are found to occur readily even at this low temperature, including concerted reorientation (flipping) of neighboring surface molecules, which changes the pattern of dangling H-atoms, and the formation of interstitial defects by the downwards motion of upper-bilayer molecules. On the proton-disordered surface, one major surface roughening process is observed that significantly disrupts the crystalline structure. Despite much longer simulation time, such roughening processes are not observed on the highly ordered Fletcher surface which is energetically more stable because of smaller repulsive interaction between neighboring dangling H-atoms. However, a more localized process takes place on the Fletcher surface involving a surface molecule transiently leaving its lattice site. The flipping process provides a facile pathway of increasing proton-order and stabilizing the surface, supporting a predominantly Fletcher-like ordering of low-temperature ice surfaces, but our simulations also show that proton- disordered patches on the surface may induce significant local reconstructions. Further, a subset of the molecules on the Fletcher surface are susceptible to forming interstitial defects.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Effect of Hydrogen Adsorption on the Magnetic Properties of a Surface Nanocluster of Iron

The effect of hydrogen adsorption on the magnetic properties of an Fe$_3$ cluster immersed in a Cu(111) surface has been calculated using densifty functional theory and the results used to parametrize an Alexander-Anderson model which takes into account the interaction of d-electrons with itinerant electrons. A number of adatom configurations containing one to seven H-atoms were analyzed. The sequential addition of hydrogen atoms is found to monotonically reduce the total magnetic moment of the cluster with the effect being strongest when the H-atoms sit at low coordinated sites. Decomposition of the charge density indicates a transfer of 0.4 electrons to each of the H-atoms from both the Fe-atoms and from the copper substrate, irrespective of adsorption site and coverage. The magnetic moment of only the nearest neighbor Fe-atoms is reduced and mainly due to increased population of minority spin d-states. This can be modeled by increased indirect coupling of d-states via the conduction s-band in the Alexander-Anderson model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The effect of the Perdew-Zunger self-interaction correction to density functionals on the energetics of small molecules

Self-consistent calculations using the Perdew-Zunger self-interaction correction (PZ-SIC) to local density and gradient dependent energy functionals are presented for the binding energy and equilibrium geometry of small molecules as well as energy barriers of reactions. The effect of the correction is to reduce binding energy and bond lengths and increase activation energy barriers when bond breaking is involved. The accuracy of the corrected functionals varies strongly, the correction to the binding energy being too weak for the local density approximation but too strong for the gradient dependent functionals considered. For the Perdew, Burke, and Ernzerhof (PBE) functional, a scaling of the PZ-SIC by one half gives improved results on average for both binding energy and bond lengths. The PZ-SIC does not necessarily give more accurate total energy, but it can result in a better cancellation of errors. An essential aspect of these calculations is the use of complex orbitals. A restriction to real orbitals leads to less accurate results as was recently shown for atoms [S. Klüpfel, P. Klüpfel, and H. Jónsson, Phys. Rev. A {84}, 050501 (2011)]. The molecular geometry of radicals can be strongly affected by PZ-SIC. An incorrect, non-linear structure of the \ce{C2H} radical predicted by PBE is corrected by PZ-SIC. The \ce{CH3} radical is correctly predicted to be planar when complex orbitals are used, while it is non-planar when the PZ-SIC calculation is restricted to real orbitals.

0 citations0 reviewsNo code linkedOpen access