No Arabic abstract
We present a new method for calculating the product yield of a radical pair recombination reaction in the presence of a weak time-dependent magnetic field. This method successfully circumvents the computational difficulties presented by a direct solution of the Liouville-von Neumann equation for a long-lived radical pair containing many hyperfine-coupled nuclear spins. Using a modified formulation of Floquet theory, treating the time-dependent magnetic field as a perturbation, and exploiting the slow radical pair recombination, we show that one can obtain a good approximation to the product yield by considering only nearly-degenerate sub-spaces of the Floquet space. Within a significant parameter range, the resulting method is found to give product yields in good agreement with exact quantum mechanical results for a variety of simple model radical pairs. Moreover it is considerably more efficient than the exact calculation, and it can be applied to radical pairs containing significantly more nuclear spins. This promises to open the door to realistic theoretical investigations of the effect of radiofrequency electromagnetic radiation on the photochemically induced radical pair recombination reactions in the avian retina which are believed to be responsible for the magnetic compass sense of migratory birds.
We construct a density-functional formalism adapted to uniform external magnetic fields that is intermediate between conventional Density Functional Theory and Current-Density Functional Theory (CDFT). In the intermediate theory, which we term LDFT, the basic variables are the density, the canonical momentum, and the paramagnetic contribution to the magnetic moment. Both a constrained-search formulation and a convex formulation in terms of Legendre--Fenchel transformations are constructed. Many theoretical issues in CDFT find simplified analogues in LDFT. We prove results concerning $N$-representability, Hohenberg--Kohn-like mappings, existence of minimizers in the constrained-search expression, and a restricted analogue to gauge invariance. The issue of additivity of the energy over non-interacting subsystems, which is qualitatively different in LDFT and CDFT, is also discussed.
Magnetic graphene nanoribbons (GNRs) have become promising candidates for future applications, including quantum technologies. Here, we characterize magnetic states hosted by chiral graphene nanoribbons (chGNRs). The substitution of a hydrogen atom at the chGNR edge by a ketone group effectively adds one p_z electron to the {pi}-electron network, thus producing an unpaired {pi} radical. A closely related scenario occurs for regular ketone-functionalized chGNRs in which one oxygen atom is missing. Two such radical states can interact via exchange coupling and we study those interactions as a function of their relative position, which includes a remarkable dependence on the chirality, as well as on the nature of the surrounding GNR, i.e., with or without ketone functionalization. In addition, we determine the parameters whereby this type of systems with oxygen heteroatoms can be adequately described within the widely used mean-field Hubbard model. Altogether, we provide new insights to both theoretically model and devise GNR-based nanostructures with tunable magnetic properties.
We present an efficient implementation of the second- and third-order single-reference algebraic diagrammatic construction theory for electron attachment (EA) and ionization (IP) energies and spectra (EA/IP-ADC(n), n = 2, 3). Our new EA/IP-ADC program features spin adaptation for closed-shell systems, density fitting for efficient handling of the two-electron integral tensors, as well as vectorized and parallel implementation of tensor contractions. We demonstrate capabilities of our efficient implementation by applying the EA/IP-ADC(n) (n = 2, 3) methods to compute the photoelectron spectrum of the TEMPO radical, as well as the vertical and adiabatic electron affinities of TEMPO and two DNA base pairs (guanine-cytosine and adenine-thymine). The spectra and electron affinities computed using large diffuse basis sets with up to 1028 molecular orbitals are found to be in a good agreement with the best available results from the experiment and theoretical simulations.
Radical pair recombination reactions are known to be sensitive to the application of both low and high magnetic fields. The application of a weak magnetic field reduces the singlet yield of a singlet-born radical pair, whereas the application of a strong magnetic field increases the singlet yield. The high field effect arises from energy conservation: when the magnetic field is stronger than the sum of the hyperfine fields in the two radicals, ${rm S}to {rm T}_{pm}$ transitions become energetically forbidden, thereby reducing the number of pathways for singlet to triplet interconversion. The low field effect arises from symmetry breaking: the application of a weak magnetic field lifts degeneracies among the zero field eigenstates and increases the number of pathways for singlet to triplet interconversion. However, the details of this effect are more subtle, and have not previously been properly explained. Here we present a complete analysis of the low field effect in a radical pair containing a single proton, and in a radical pair in which one of the radicals contains a large number of hyperfine-coupled nuclear spins. We find that the new transitions that occur when the field is switched on are between ${rm S}$ and ${rm T}_0$ in both cases, and not between ${rm S}$ and ${rm T}_{pm}$ as has previously been claimed. We then illustrate this result by using it in conjunction with semiclassical spin dynamics simulations to account for the observation of a biphasic--triphasic--biphasic transition with increasing magnetic field strength in the magnetic field effect on the time-dependent survival probability of a photoexcited carotenoid-porphyrin-fullerene radical pair.
An implementation of coupled-cluster (CC) theory to treat atoms and molecules in finite magnetic fields is presented. The main challenges stem from the magnetic-field dependence in the Hamiltonian, or, more precisely, the appearance of the angular momentum operator, due to which the wave function becomes complex and which introduces a gauge-origin dependence. For this reason, an implementation of a complex CC code is required together with the use of gauge-including atomic orbitals to ensure gauge-origin independence. Results of coupled-cluster singles--doubles--perturbative-triples (CCSD(T)) calculations are presented for atoms and molecules with a focus on the dependence of correlation and binding energies on the magnetic field.