Korsholmcarlsen3946

It is found that for the former three cases, our automatic algorithm can reproduce exactly the same MPOs as the optimally hand-crafted ones already known in the literature.Studies using molecular dynamics (MD) have long struggled to simulate the failure modes of materials, predicting unrealistically high ductility and failing to capture brittle fracture. The primary cause of this shortcoming is an inadequate description of bond breaking. While reactive force fields such as ReaxFF show improvements compared to traditional force fields, the charge models used yield unphysical partial charges, especially during dissociation of ionic bonds. This flaw may be remedied by using the atom-condensed Kohn-Sham density functional theory (DFT) approximated to a second order (ACKS2) charge model for determining partial charges. In this work, we present a new ACKS2-enabled Reax force field for fracture simulations of lithium oxide systems, which was obtained by training against an extensive set of DFT, multireference configuration interaction (MRCI), and MRCI+Q reference data using genetic optimization techniques. This new force field significantly improves the bond breaking behavior, but still cannot fully capture the brittle fracture in MD simulations, suggesting more research is needed to improve simulation of brittle fracture.A scheme for the calculation of molecular properties within the framework of unitary coupled-cluster (UCC) theory in both the electronic ground and excited states is presented. The scheme is based on an expectation-value ansatz, similar to the equation-of-motion coupled-cluster method or the intermediate state representation (ISR) approach of the algebraic-diagrammatic construction (ADC) scheme. Due to the UCC ansatz, the resulting equations cannot be given by closed-form expressions but need to be approximated. Explicit expressions for the expectation value of a general one-particle operator correct through second order in perturbation theory have been derived and coded for the electronic ground state as well as for excited states of predominant single-excitation character. The resulting equations are shown to be equivalent to those of the second-order ADC/ISR procedure. As first computational tests, the second-order UCC method (UCC2) and the one employing third-order amplitudes (also eigenvectors) together with the second-order density matrix, denoted as UCC3(2), are applied to the calculation of dipole moments for a series of small closed- and open-shell systems as well as 4-cyanoindole and 2,3-benzofuran and compared to full configuration interaction or experimental results. For the aromatic organic molecules, the UCC2 method is shown to be sufficient for the ground-state dipole moment, whereas the UCC3(2) scheme is superior for excited-state dipole moments.Directional control over surface plasmon polariton (SPP) waves is a prerequisite for the development of miniaturized optical circuitry. Here, the efficacy of single and dual component SPP steering elements is explored through photoemission electron microscopy. Our imaging scheme relies on two-color photoemission and counter-propagating SPP generation, which collectively allow SPPs to be visualized in real space. Wave-vector difference mixing between the two-dimensional nanohole array and photon momenta enables SPP steering with directionality governed by the array lattice constant and input photon direction. In our dual component configuration, separate SPP generation and Bragg diffraction based steering optics are employed. We find that array Bragg planes principally influence the SPP angles through the array band structure, which allows us to visualize both positive and negative refractory waves.A theoretical study of the mechanisms of electroluminescence (EL) generation in photoactive molecules with donor and acceptor centers linked by saturated σ-bonds (molecules of the Aviram-Ratner-type) is presented. The approach is based on the kinetics of single-electron transitions between many-body molecular states. This study shows that the EL polarity arises due to asymmetric coupling of molecular orbitals of the photochromic part of the molecule to the electrodes. The gate voltage controls the power of the EL through the occupancy of the excited singlet state. The shifting of the orbital energies forms a resonant or a non-resonant path for the transmission of electrons through the molecule. The action of the gate voltage is reflected in specific critical voltages. An analytical dependence of the critical voltages on the energies of molecular states involved in the formation of EL, as well as on the gate voltage, was derived for both positive and negative polarities. Conditions under which the gate voltage lowers the absolute value of the bias voltage that is responsible for the activation of the resonance mechanism of EL formation were also established. This is an important factor in control of EL in molecular junctions.We introduce a mean-field theoretical framework for generalizing isotropic pair potentials to anisotropic shapes. This method is suitable for generating pair potentials that can be used in both Monte Carlo and molecular dynamics simulations. We demonstrate the application of this theory by deriving a Lennard-Jones (LJ)-like potential for arbitrary geometries along with a Weeks-Chandler-Anderson-like repulsive variant, showing that the resulting potentials behave very similarly to standard LJ potentials while also providing a nearly conformal mapping of the underlying shape. We then describe an implementation of this potential in the simulation engine HOOMD-blue and discuss the challenges that must be overcome to achieve a sufficiently robust and performant implementation. The resulting potential can be applied to smooth geometries like ellipsoids and to convex polytopes. We contextualize these applications with reference to the existing methods for simulating such particles. The pair potential is validated using standard criteria, and its performance is compared to existing methods for comparable simulations. Finally, we show the results of self-assembly simulations, demonstrating that this method can be used to study the assembly of anisotropic particles into crystal structures.This work proposes to describe open-shell molecules or radicals using the framework of the doubly occupied configuration interaction (DOCI) treatments, so far limited to closed-shell system studies. The proposal is based on considering molecular systems in singlet states generated by adding extra hydrogen atoms located at infinite distance from the target radical system. The energy of this radical is obtained by subtracting the energies of the dissociated hydrogen atoms from that provided by the two-electron reduced density matrix corresponding to the singlet state system in the DOCI space, which is variationally calculated by imposing a set of N-representability conditions. This method is numerically assessed by describing potential energy curves and reduced density matrices in selected ionic and neutral open-shell systems in the doublet spin symmetry ground state.In this paper, we investigate the substrate effect in graphene temperature sensors. Recently, there have been many research studies done on temperature sensors using the nanofabrication technique. However, the sensitivity and response time need to be improved. In this study, we propose a new type of temperature sensor that consists of graphene and Anodic Aluminum Oxide (AAO). In this device, graphene and AAO are used as the sensing material and the substrate, respectively. We characterize the sensitivity and the response time using the experimental results and simulation data. The real-time resistance change of graphene is monitored depending on the temperature, and the response time is also analyzed by COMSOL Multiphysics. To confirm the porous substrate effect, we compare the device performance of the AAO substrate to the performance of the glass substrate. From these results, the suspended graphene on the AAO substrate shows about two times higher sensitivity and a much faster response time than the glass substrate.We carry out molecular dynamics simulations by using an all-atom model to study the nucleation and crystallization of n-alkane droplets under three-dimensional and quasi-two-dimensional conditions. We focus on the development of orientational order of chains from a random state to a neatly ordered one. Two new methods, the map of symmetry breaking and the information entropy of chain orientations, are introduced to characterize the emerge and remelting phenomena of a primary nucleus at the early stage of crystallization. Stepwise nucleation, as well as the surface induced nucleation, of large droplets is observed. We elucidate the kinetic process of the formation of a primary nucleus and the rearrangement of every single molecule involved in a primary nucleus. We found that density fluctuation and orientational preordering are coupled together and occur simultaneously in nucleation. Our results show the pathway of orientational symmetry breaking in the crystallization of n-alkane droplets that are heuristic for the deeper understanding of the crystallization in more complex molecules such as polymers.Internal conversion decay dynamics associated with the potential energy surfaces of three low-lying singlet excited electronic states, S1 (ππ*, A'), S2 (ππ*, A'), and S3 (nπ*, A″), of tropolone are investigated theoretically. Energetic and spatial aspects of conical intersections of these electronic states are explored with the aid of the linear vibronic coupling approach. Symmetry selection rules suggest that non-totally symmetric modes would act as coupling modes between S1 and S3 as well as between S2 and S3. We found that the S1-S2 interstate coupling via totally symmetric modes is very weak. A diabatic vibronic Hamiltonian consisting of 32 vibrational degrees of freedom is constructed to simulate the photoinduced dynamics of S0 → S1 and S0 → S2 transitions. We observe a direct nonadiabatic population transfer from S1 to S3, bypassing S2, during the initial wavepacket propagation on S1. On the other hand, the initial wavepacket evolving on S2 would pass through the S2-S3 and S1-S3 conical intersections before reaching S1. The presence of multiple proton transfer channels on the S1-S2-S3 coupled potential energy surfaces of tropolone is analyzed. Our findings necessitate the treatment of proton tunneling dynamics of tropolone beyond the adiabatic symmetric double well potentials.In this paper, we have analyzed the time series associated with the iterative scheme of a double similarity transformed coupled cluster theory. The coupled iterative scheme to solve the ground state Schrödinger equation is cast as a multivariate time-discrete map, and the solutions show the universal Feigenbaum dynamics. Using recurrence analysis, it is shown that the dynamics of the iterative process is dictated by a small subgroup of cluster operators, mostly those involving chemically active orbitals, whereas all other cluster operators with smaller amplitudes are enslaved. Using synergetics, we will indicate how the master-slave dynamics can suitably be exploited to develop a novel coupled-cluster algorithm in a much reduced dimension.