Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics.

The exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron-ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron-nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.

Direct Dynamics with Nuclear-Electronic Orbital Density Functional Theory.

Direct dynamics simulations of chemical reactions typically require the selection of a method for generating the potential energy surfaces and a method for the dynamical propagation of the nuclei on these surfaces. The nuclear-electronic orbital (NEO) framework avoids this Born-Oppenheimer separation by treating specified nuclei on the same level as the electrons with wave function methods or density functional theory (DFT). The NEO approach is particularly applicable to proton, hydride, and proton-coupled electron transfer reactions, where the transferring proton(s) and all electrons are treated quantum mechanically. In this manner, the zero-point energy, density delocalization, and anharmonicity of the transferring protons are inherently and efficiently included in the energies, optimized geometries, and dynamics.This Account describes how various NEO methods can be used for direct dynamics simulations on electron-proton vibronic surfaces. The strengths and limitations of these approaches are discussed, and illustrative examples are presented. The NEO-DFT method can be used to simulate chemical reactions on the ground state vibronic surface, as illustrated by the application to hydride transfer in C4H9+. The NEO multistate DFT (NEO-MSDFT) method is useful for simulating ground state reactions in which the proton density becomes bilobal during the dynamics, a characteristic of hydrogen tunneling, as illustrated by proton transfer in malonaldehyde. The NEO time-dependent DFT (NEO-TDDFT) method produces excited electronic, vibrational, and vibronic surfaces. The application of linear-response NEO-TDDFT to H2 and H3+, as well as the partially and fully deuterated counterparts, shows that this approach produces accurate fundamental vibrational excitation energies when all nuclei and all electrons are treated quantum mechanically. Moreover, when only specified nuclei are treated quantum mechanically, this approach can be used to optimize geometries on excited state vibronic surfaces, as illustrated by photoinduced single and double proton transfer systems, and to conduct adiabatic dynamics on these surfaces. The real-time NEO-TDDFT method provides an alternative approach for simulating nonequilibrium nuclear-electronic dynamics of such systems. These various NEO methods can be combined with nonadiabatic dynamics methods such as Ehrenfest and surface hopping dynamics to include the nonadiabatic effects between the quantum and classical subsystems. The real-time NEO-TDDFT Ehrenfest dynamics simulation of excited state intramolecular proton transfer in o-hydroxybenzaldehyde illustrates the power of this type of combined approach. The field of multicomponent quantum chemistry is in the early stages, and the methods discussed herein provide the foundation for a wide range of promising future directions to be explored. An appealing future direction is the expansion of the real-time NEO-TDDFT method to describe the dynamics of all nuclei and electrons on the same level. Direct dynamics simulations using NEO wave function methods such as equation-of-motion coupled cluster or multiconfigurational approaches are also attractive but computationally expensive options. The further development of NEO direct dynamics methods will enable the simulation of the nuclear-electronic dynamics for a vast array of chemical and biological processes that extend beyond the Born-Oppenheimer approximation.

Psychological, physical and social factors influence decision to return to sport after revision ACL reconstruction with BPTB graft.

PURPOSE: There is limited evidence in literature regarding the patient-reported factors that influence their return to sport (RTS) in revision anterior cruciate ligament reconstruction (ACLR). The medium-term results of a prospective consecutive cohort of patients undergoing single- and two-stage revision ACLR with bone patellar tendon bone graft (BPTB) and patient-reported factors that influence their decision to return to sport are presented in this study. METHODS: Seventy-two patients were included in this prospective study. Single- or two-stage revision with BPTB graft was performed based on pre-operative planning. Iliac crest bone graft was used. Pre-operative and follow-up Lysholm and Tegner activity scores and RTS, level of sport and patient-reported factors affecting RTS were recorded. The mean follow-up was 9 years (SD 2.7 years). RESULTS: Single-stage revision ACLR was performed in 61 patients. In 11 patients (15%), revision ACLR was performed in two stages. There was a significant improvement in Lysholm score from mean 51.1 to 86.7 (p < 0.001). The incidence of re-rupture in this cohort was 0%. The median Tegner score was 6 (range 2-9). Twenty-five patients (34.7%) did not return to any sport at final follow-up. Twenty-nine (40.2%) patients returned to their pre-injury level of sport. Fear of reinjury (79%, p < 0.001) and persistent knee symptoms (35.8%, p = 0.03) were the most common factors limiting RTS in non-returners. CONCLUSION: Psychological and social factors may have an influence on RTS in addition to physical factors. LEVEL OF EVIDENCE: Level III.

Subcutaneous implantable cardioverter defibrillator: Can it overtake its transvenous counterpart.

Over the past decade, the emergence of the subcutaneous implantable cardioverter defibrillator (S-ICD) has provided cardiologists with an option to provide both primary or secondary prevention treatment of sudden cardiac death (SCD) without the associated risks that come with the use of intracardiac leads. S-ICD may prove to be a useful option in those who are young, have thromboembolic risk, immunodeficiency states, unfavorable anatomy due to adult congenital heart disease (ACHD). This article reviews the existing literature to determine whether S-ICD can prove to be a safe alternative in comparison to Transvenous implantable cardioverter defibrillator (TV-ICD) and in which patient population should S-ICD be considered over TV-ICD.

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations.

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

Synthesis of LnII -in-Cryptand Complexes by Chemical Reduction of LnIII -in-Cryptand Precursors: Isolation of a NdII -in-Cryptand Complex.

Lanthanide triflates have been used to incorporate NdIII and SmIII ions into the 2.2.2-cryptand ligand (crypt) to explore their reductive chemistry. The Ln(OTf)3 complexes (Ln=Nd, Sm; OTf=SO3 CF3 ) react with crypt in THF to form the THF-soluble complexes [LnIII (crypt)(OTf)2 ][OTf] with two triflates bound to the metal encapsulated in the crypt. Reduction of these LnIII -in-crypt complexes using KC8 in THF forms the neutral LnII -in-crypt triflate complexes [LnII (crypt)(OTf)2 ]. DFT calculations on [NdII (crypt)]2+ ], the first NdII cryptand complex, assign a 4f4 electron configuration to this ion.

Isolation of a Square-Planar Th(III) Complex: Synthesis and Structure of [Th(OC6H2tBu2-2,6-Me-4)4]1.

Reduction of Th(OC6H2tBu2-2,6-Me-4)4 using either KC8 or Li in THF forms a new example of a crystallographically characterizable Th(III) complex in the salts [K(THF)5(Et2O)][Th(OC6H2tBu2-2,6-Me-4)4] and [Li(THF)4][Th(OC6H2tBu2-2,6-Me-4)4]. Surprisingly, in each structure the four aryloxide ligands are arranged in a square-planar geometry, the first example of this coordination mode for an f element complex. The Th(III) ion and four oxygen donor atoms are coplanar to within 0.05 Å with O-Th-O angles of 89.27(8) to 92.02(8)° between cis ligands. The ligands have Th-O-C(ipso) angles of 173.9(2) to 178.6(4)°, and the aryl rings make angles of 58.5 to 65.1° with the ThO4 plane. The effect of the eight tert-butyl substituents in generating the unusual structure through packing and/or dispersion forces is discussed. EPR spectroscopy reveals an axial signal consistent with a metal-based radical in a planar complex. DFT calculations yield a C4-symmetric structure that accommodates a low-lying SOMO of 6dz2 character with 7s Rydberg admixture.

Predicting Reaction Products and Automating Reactive Trajectory Characterization in Molecular Simulations with Support Vector Machines.

A machine learning-based methodology for the prediction of chemical reaction products, along with automated elucidation of mechanistic details via phase space analysis of reactive trajectories, is introduced using low-dimensional heuristic models and then applied to ab initio computer simulations of the photodissociation of acetaldehyde, an important chemical system in atmospheric chemistry. Our method is centered around training Support Vector Machines (SVMs) to identify optimal separatrices that delineate the regions of phase space that lead to different chemical reaction products. In contrast to the more common "black box" type machine learning methodologies for analyzing chemical simulation data, this SVM-based methodology allows for mechanistic insight to be gleaned from further analysis of the positioning of the phase space points used to train the SVM with respect to the separatrices. For example, a pair of phase space points that are in close proximity to each other but on opposite sides of a separatrix may be situated on opposite sides of a transition state, while phase space points occurring early in a simulation that are distant from a separatrix are likely to belong to trajectories that are highly biased toward the product state associated with the basin of attraction to which they belong. In addition to inferring mechanistic details about multiple-pathway chemical reactions, our method can also be used to increase reactive trajectory sampling efficiency in molecular simulations via rejection sampling, with trajectories leading to undesired product states being identified and terminated early in the simulation rather than being carried to completion.

Multistate hybrid time-dependent density functional theory with surface hopping accurately captures ultrafast thymine photodeactivation.

We report an efficient analytical implementation of first-order nonadiabatic derivative couplings between arbitrary Born-Oppenheimer states in the hybrid time-dependent density functional theory (TDDFT) framework using atom-centered basis functions. Our scheme is based on quadratic response theory and includes orbital relaxation terms neglected in previous approaches. Simultaneous computation of multiple derivative couplings and energy gradients enables efficient multistate nonadiabatic molecular dynamics simulations in conjunction with Tully's fewest switches surface hopping (SH) method. We benchmark the thus obtained multistate TDDFT-SH scheme by simulating ultrafast decay of UV-photoexcited thymine, for which accurate gas-phase data from ultrafast spectroscopy experiments are available. The calculations predict a fast 153 fs decay from the bright S2 to the dark S1 excited state, followed by a much slower 14 ps S1 deactivation to the ground state; statistical uncertainties were estimated using bootstrap sampling. These results agree well with the experimentally observed time constants of 100-200 fs and 5-7 ps, respectively, unlike previous multiconfigurational self-consistent field and second-order algebraic diagrammatic construction calculations. Furthermore, our results support the S1-trapping hypothesis [J. J. Szymczak et al., J. Phys. Chem. A, 2009, 113, 12686-12693]. For thymine, the computational cost of a single TDDFT-SH time-step including the lowest 3 states, all couplings and gradients, is ∼5 times larger than the cost of a single Born-Oppenheimer dynamics time step for the ground state in our implementation. Thus, ps nonadiabatic dynamics simulations using multistate hybrid TDDFT-SH for systems with up to ∼100 atoms are possible without drastic approximations on single workstation nodes. Our implementation will be made available through Turbomole.

5-Methoxyquinoline Photobasicity Is Mediated by Water Oxidation.

A mechanism explaining the photobasicity of 5-methoxyquinoline (5-MeOQ) is proposed on the basis of nonadiabatic molecular dynamics simulations using time-dependent hybrid density functional theory (TDDFT) and fewest switches surface hopping (SH) and analysis of existing ultrafast spectroscopy experiments. According to the TDDFT-SH simulations, the rate-determining step is hole transfer from photoexcited 5-MeOQ to adjacent water molecules within ∼5 ps followed by rapid electron-coupled proton transfer and deactivation to the ground state. This fast redox-catalyzed proton transfer mechanism is consistent with simple thermodynamic arguments, correlated wavefunction calculations, and recent isotope substitution time-resolved fluorescence experiments. Although energetically feasible, direct protonation of 5-MeOQ in the S1 state appears to be too slow to contribute significantly to 5-MeOQ photobasicity in aqueous solution.

Unphysical divergences in response theory.

Transition densities between excited states are key for nonlinear theoretical spectroscopy and multi-state non-adiabatic molecular dynamics (NAMD) simulations. In the framework of response theory, these transition densities are accessible from poles of the quadratic response function. It was shown recently that the thus obtained transition densities within time-dependent Hartree-Fock (TDHF) and adiabatic time-dependent density functional theory (TDDFT) exhibit unphysical divergences when the difference in excitation energy of the two states of interest matches another excitation energy. This unphysical behavior is a consequence of spurious poles in the quadratic response function. We show that the incorrect pole structure of the quadratic response is not limited to TDHF and adiabatic TDDFT, but is also present in many other approximate many-electron response functions, including those from coupled cluster and multiconfigurational self-consistent field response theory. The divergences appear in regions of the potential energy surface where the ground state is perfectly well behaved, and they are frequently encountered in NAMD simulations of photochemical reactions. The origin of the divergences is traced to an incorrect instantaneous time-dependence of the effective Hamiltonian. The implications for computations of frequency-dependent response properties are considerable and call into question the validity of conventional approximate many-electron response theories beyond linear response.

Curing the Divergence in Time-Dependent Density Functional Quadratic Response Theory.

The adiabatic approximation in time-dependent density functional theory is known to give an incorrect pole structure in the quadratic response function, leading to unphysical divergences in excited state-to-state transition probabilities and hyperpolarizabilties. We find the form of the exact quadratic response kernel and derive a practical and accurate approximation that cures the divergence. We demonstrate our results on excited state-to-state transition probabilities of a model system and of the LiH molecule.

Non-linear Hall effect in multi-Weyl semimetals.

In the presence of time reversal symmetry, a non-linear Hall effect can occur in systems without an inversion symmetry. One of the prominent candidates for detection of such Hall signals are Weyl semimetals. In this article, we investigate the Berry curvature induced second and third order Hall effect in multi-Weyl semimetals with topological chargesn=1,2,3. We use low energy effective models to obtain general analytical expressions and discover the presence of a large Berry curvature dipole (BCD) in multi-Weyl semimetals, compared to usual (n = 1) Weyl semimetals. We also study the BCD in a realistic tight-binding lattice model and observe two different kinds of variation with increasing topological charge-these can be attributed to different underlying Berry curvature components. We provide estimates of the signatures of second harmonic of Hall signal in multi-Weyl semimetals, which can be detected experimentally. Furthermore, we predict the existence of a third order Hall signal in multi-Weyl semimetals. We derive the analytical expressions of Berry connection polarizability tensor, which is responsible for third order effects, using a low energy model and estimate the measurable conductivity. Our work can help guide experimental discovery of Berry curvature multipole physics in multi-Weyl semimetals.

Nonadiabatic Dynamics of Hydrogen Tunneling with Nuclear-Electronic Orbital Multistate Density Functional Theory.

Proton transfer reactions play a critical role in many chemical and biological processes. The development of computationally efficient approaches to describe the quantum dynamics of proton transfer, which often involves hydrogen tunneling, is challenging. Herein, the nuclear-electronic orbital multistate density functional theory (NEO-MSDFT) method is combined with both Ehrenfest and surface hopping nonadiabatic dynamics methods to describe hydrogen tunneling. The NEO-MSDFT method treats the transferring hydrogen nucleus quantum mechanically on the same level as the electrons and incorporates both static and dynamical correlation by mixing localized NEO-DFT solutions with a nonorthogonal configuration interaction scheme. The other nuclei are propagated on the NEO-MSDFT vibronic surfaces during the Ehrenfest or surface hopping dynamics. These methods are applied to proton transfer in malonaldehyde as a prototypical hydrogen tunneling system. The inclusion of vibronically nonadiabatic effects is found to significantly impact the proton transfer time and tunneling dynamics. This approach is applicable to a wide range of other proton transfer reactions.

The Role of 18F-Flortaucipir (AV-1451) in the Diagnosis of Neurodegenerative Disorders.

Tau protein plays a vital role in maintaining the structural and functional integrity of the nervous system; however, hyperphosphorylation or abnormal phosphorylation of tau protein plays an essential role in the pathogenesis of several neurodegenerative disorders. The development of radioligand such as the 18F-flortaucipir (AV-1451) has provided us with the opportunity to assess the underlying tau pathology in various etiologies of dementia. For the purpose of this article, we aimed to evaluate the utility of 18F-AV-1451 in the differential diagnosis of various neurodegenerative disorders. We used PubMed to look for the latest, peer-reviewed, and informative articles. The scope of discussion included the role of 18F-AV-1451 positron emission tomography (PET) to aid in the diagnosis of Alzheimer's disease (AD), frontotemporal dementia (FTD), dementia with Lewy bodies (DLB), and Parkinson's disease with dementia (PDD). We also discussed if the tau burden identified by neuroimaging correlated well with the clinical severity and identified the various challenges of 18F-AV-1451 PET. We concluded that although the role of 18F-AV-1451 seems promising in the neuroimaging of AD, the benefit appears uncertain when it comes to the non-Alzheimer's tauopathies. More research is required to identify the off-target binding sites of 18F-AV-1451 to determine its clinical utility in the future.

Analytical Gradients for Nuclear-Electronic Orbital Time-Dependent Density Functional Theory: Excited-State Geometry Optimizations and Adiabatic Excitation Energies.

The computational investigation of photochemical processes often entails the calculation of excited-state geometries, energies, and energy gradients. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons, thereby including the associated nuclear quantum effects and non-Born-Oppenheimer behavior into quantum chemistry calculations. The multicomponent density functional theory (NEO-DFT) and time-dependent DFT (NEO-TDDFT) methods allow efficient calculations of ground and excited states, respectively. Herein, the analytical gradients are derived and implemented for the NEO-TDDFT method and the associated Tamm-Dancoff approximation (NEO-TDA). The programmable equations for these analytical gradients as well as the NEO-DFT analytical Hessian are provided. The NEO approach includes the anharmonic zero-point energy (ZPE) and density delocalization associated with the quantum protons as well as vibronic mixing in geometry optimizations and energy calculations of ground and excited states. The harmonic ZPE associated with the other nuclei can be computed via the NEO Hessian. This approach is used to compute the 0-0 adiabatic excitation energies for a set of nine small molecules with all protons quantized, exhibiting slight improvement over the conventional electronic approach. Geometry optimizations of two excited-state intramolecular proton-transfer systems, [2,2'-bipyridyl]-3-ol and [2,2'-bipyridyl]-3,3'-diol, are performed with one and two quantized protons, respectively. The NEO calculations for these systems produce electronically excited-state geometries with stronger intramolecular hydrogen bonds and similar relative stabilities compared to conventional electronic methods. This work provides the foundation for nonadiabatic dynamics simulations of fundamental processes such as photoinduced proton transfer and proton-coupled electron transfer.

Distraction osteogenesis in the pediatric population.

OBJECTIVES: Distraction osteogenesis has been described routinely in the mandible for the advancement of bony segments. Complications, though infrequent, may include postoperative infection, implant extrusion, nonunion of the bony segments, device malfunction, cranial nerve paresis, and premature consolidation. METHODS: Seventy-eight distractions of the mandible were performed over 10 years. The technique in placement of these internal microdistraction devices involves making intraoral and extraoral incisions and beginning distraction after a latency period of 3 days. Following this latency period, distraction occurs at 1 mm/day. RESULTS: In this series of patients, distraction was accomplished successfully. There was a 2.6% rate of wound infection in this series; 2.6% also had bony nonunion; 3.8% had premature bony consolidation; and 3.8% had facial nerve complications. In the three patients with facial nerve paresis, we followed the facial nerve clinically and each paresis resolved within 6 months. CONCLUSIONS: Distraction of the mandible may be accomplished in the pediatric population. As with any intervention, inherent perioperative complications may still arise.

Transcanal cochlear implantation under monitored anesthesia care.

OBJECTIVE: Given the associated risk of general anesthesia in elderly patients with cardiovascular disease, the authors set out to determine the feasibility of transcanal cochlear implantation under local anesthesia with monitored anesthesia care. METHODS: A 70-year-old man with a history of coronary artery bypass grafting, diabetes mellitus, and an American Society of Anesthesiologists Class III cardiac status underwent cochlear implantation under local with monitored anesthesia care. RESULT: With the described technique and regimen of intravenous remifentanil and dexmedetomidine, the patient tolerated the 60-minute procedure without tachycardia, hyper- or hypotension, or cardiac ischemia. CONCLUSION: Cochlear implantation using the pericanal electrode technique performed under local anesthesia with monitored anesthesia care is possible in patients at risk for undergoing general anesthesia for cochlear implantation.

Ultrasound-aided fixation of a biodegradable cranial fixation system: uses in pediatric neurosurgery.

OBJECT: Bioresorbable implant systems have been used in neurosurgery for the rigid fixation of cranial and facial bones. A relatively recent advancement has been the fixation of these implants using an ultrasonic device. The experience with such a device in neurosurgical practice has been limited. The authors report on their experience with ultrasound-aided fixation of bioresorbable implants in pediatric neurosurgical practice. METHODS: The study consisted of 2 parts. The retrospective portion consisted of a chart review of pertinent clinical information, complications, and outcomes after the use of a commercially available ultrasound-aided bioresorbable implant system (SonicWeld Rx, KLS Martin L.P.). Follow-up was obtained in all patients via clinical examination or telephone interview. The prospective portion of the study consisted of video analysis of the implantation technique in a routine craniotomy. Implantation times were measured, and delays during treatment were noted. RESULTS: Over a period of 2 years, 28 consecutive patients underwent placement of these implants for bone fixation during craniotomies or craniofacial reconstructions. The only complication was seen in a child with Crouzon syndrome, who had a wound infection caused by Serratia sepsis from a central venous line infection. There were no repeated operations for implant-related swelling, and no cases of premature plate resorption, bone instability, or settling. In vivo, the average time required to implant a resorbable pin with this system was 22 seconds. CONCLUSIONS: The use of a bioresorbable implant system with ultrasound-aided pin fixation in pediatric neurosurgery cases achieved adequate stability with few complications. This system was easy to use and provided rapid fixation of implants.