Beside a Garden Wall

With Ukulele arrangement. Contains advertisements and/or short musical examples of pieces being sold by publisher.https://digitalcommons.library.umaine.edu/mmb-vp/7172/thumbnail.jp

Butterfly Hysteresis and Slow Relaxation of the Magnetization in (Et4N)3Fe2F9: Manifestations of a Single-Molecule Magnet

(Et4N)3Fe2F9 exhibits a butterfly--shaped hysteresis below 5 K when the magnetic field is parallel to the threefold axis, in accordance with a very slow magnetization relaxation in the timescale of minutes. This is attributed to an energy barrier Delta=2.40 K resulting from the S=5 dimer ground state of [Fe2F9]^{3-} and a negative axial anisotropy. The relaxation partly occurs via thermally assisted quantum tunneling. These features of a single-molecule magnet are observable at temperatures comparable to the barrier height, due to an extremely inefficient energy exchange between the spin system and the phonons. The butterfly shape of the hysteresis arises from a phonon avalanche effect.Comment: 18 pages, 5 eps figures, latex (elsart

Sequential Quasi-Monte Carlo

We derive and study SQMC (Sequential Quasi-Monte Carlo), a class of algorithms obtained by introducing QMC point sets in particle filtering. SQMC is related to, and may be seen as an extension of, the array-RQMC algorithm of L'Ecuyer et al. (2006). The complexity of SQMC is $O(N \log N)$, where $N$ is the number of simulations at each iteration, and its error rate is smaller than the Monte Carlo rate $O_P(N^{-1/2})$. The only requirement to implement SQMC is the ability to write the simulation of particle $x_t^n$ given $x_{t-1}^n$ as a deterministic function of $x_{t-1}^n$ and a fixed number of uniform variates. We show that SQMC is amenable to the same extensions as standard SMC, such as forward smoothing, backward smoothing, unbiased likelihood evaluation, and so on. In particular, SQMC may replace SMC within a PMCMC (particle Markov chain Monte Carlo) algorithm. We establish several convergence results. We provide numerical evidence that SQMC may significantly outperform SMC in practical scenarios.Comment: 55 pages, 10 figures (final version

Pulmonary vein flow split effects in patient-specific simulations of left atrial flow

Disruptions to left atrial (LA) blood flow, such as those caused by atrial fibrillation (AF), can lead to thrombosis in the left atrial appendage (LAA) and an increased risk of systemic embolism. LA hemodynamics are influenced by various factors, including LA anatomy and function, and pulmonary vein (PV) inflow conditions. In particular, the PV flow split can vary significantly among and within patients depending on multiple factors. In this study, we investigated how changes in PV flow split affect LA flow transport, focusing for the first time on blood stasis in the LAA, using a high-fidelity patient-specific computational fluid dynamics (CFD) model. We use an Immersed Boundary Method, simulating the flow in a fixed, uniform Cartesian mesh and imposing the movement of the LA walls with a moving Lagrangian mesh generated from 4D Computerized Tomography images. We analyzed LA anatomies from eight patients with varying atrial function, including three with AF and either a LAA thrombus or a history of Transient Ischemic Attacks (TIAs). Using four different flow splits (60/40% and 55/45% through right and left PVs, even flow rate, and same velocity through each PV), we found that flow patterns are sensitive to PV flow split variations, particularly in planes parallel to the mitral valve. Changes in PV flow split also had a significant impact on blood stasis and could contribute to increased risk for thrombosis inside the LAA, particularly in patients with AF and previous LAA thrombus or a history of TIAs. Our study highlights the importance of considering patient-specific PV flow split variations when assessing LA hemodynamics and identifying patients at increased risk for thrombosis and stroke. This knowledge is relevant to planning clinical procedures such as AF ablation or the implementation of LAA occluders.This work was partially supported by Comunidad de Madrid (Synergy Grant Y2018/BIO-4858 PREFI-CM), Spanish Research Agency (AEI, grant number PID2019-107279RB-I00), Instituto de Salud Carlos III (grant numbers PI15/02211-ISBITAMI and DTS/1900063-ISBIFLOW), and by the EU-European Regional Development Fund. Funding for open access charge: Universidad de Málaga /CBU

The Hubble flow around the Local Group

We use updated data on distances and velocities of galaxies in the proximity of the Local Group (LG) in order to establish properties of the local Hubble flow. For 30 neighbouring galaxies with distances 0.7 < D_LG < 3.0 Mpc, the Local flow is characterized by the Hubble parameter H_loc = (78+/-2) km/(s*Mpc), the mean-square peculiar velocity sigma_v = 25 km/s, corrected for errors of radial velocity measurements (~4 km/s) and distance measurements (~10 km/s), as well as the radius of the zero-velocity surface R_0 = (0.96+/-0.03) Mpc. The minimum value for sigma_v is achieved when the barycenter of the LG is located at the distance D_c = (0.55+/-0.05) D_M31 towards M31 corresponding to the Milky Way-to-M31 mass ratio M_MW / M_M31 ~ 4/5. In the reference frame of the 30 galaxies at 0.7 - 3.0 Mpc, the LG barycenter has a small peculiar velocity ~(24+/-4) km/s towards the Sculptor constellation. The derived value of R_0 corresponds to the total mass M_T(LG) = (1.9+/-0.2) 10^12 M_sun with Omega_m = 0.24 and a topologically flat universe, a value in good agreement with the sum of virial mass estimates for the Milky Way and M31.Comment: 14 pages, 6 figures, 1 table. Accepted for publication in MNRA

Efficient multi-fidelity computation of blood coagulation under flow

Clot formation is a crucial process that prevents bleeding, but can lead to severe disorders when imbalanced. This process is regulated by the coagulation cascade, a biochemical network that controls the enzyme thrombin, which converts soluble fibrinogen into the fibrin fibers that constitute clots. Coagulation cascade models are typically complex and involve dozens of partial differential equations (PDEs) representing various chemical species’ transport, reaction kinetics, and diffusion. Solving these PDE systems computationally is challenging, due to their large size and multi-scale nature. We propose a multi-fidelity strategy to increase the efficiency of coagulation cascade simulations. Leveraging the slower dynamics of molecular diffusion, we transform the governing PDEs into ordinary differential equations (ODEs) representing the evolution of species concentrations versus blood residence time. We then Taylor-expand the ODE solution around the zero-diffusivity limit to obtain spatiotemporal maps of species concentrations in terms of the statistical moments of residence time, , and provide the governing PDEs for . This strategy replaces a high-fidelity system of N PDEs representing the coagulation cascade of N chemical species by N ODEs and p PDEs governing the residence time statistical moments. The multi-fidelity order (p) allows balancing accuracy and computational cost providing a speedup of over N/p compared to high-fidelity models. Moreover, this cost becomes independent of the number of chemical species in the large computational meshes typical of the arterial and cardiac chamber simulations. Using a coagulation network with N = 9 and an idealized aneurysm geometry with a pulsatile flow as a benchmark, we demonstrate favorable accuracy for low-order models of p = 1 and p = 2. The thrombin concentration in these models departs from the high-fidelity solution by under 20% (p = 1) and 2% (p = 2) after 20 cardiac cycles. These multi-fidelity models could enable new coagulation analyses in complex flow scenarios and extensive reaction networks. Furthermore, it could be generalized to advance our understanding of other reacting systems affected by flow.MGH, MGV and OF have been partially supported by the Spanish Research Agency and the European Regional Development Fund, under grant number PID2019-107279RB-I00. MGH, MGV, PML, JB and OF have been partially supported by the Comunidad de Madrid and the European Regional Development Fund, under grant number Y2018/BIO-4858 PREFI-CM, and by the Instituto de Salud Carlos III and the European Regional Development Fund, under grant numbers PI15/02211-ISBITAMI and DTS/1900063-ISBIFLOW. AG, EMcV, AK and JCdA have been partially supported by the US National Institutes of Health, under grant 1R01HL160024. JCdA has been partially supported by the US National Insitutes of Health, under grant number 1R01HL158667

Pulmonary vein flow split effects in patient-specific simulations of left atrial flow

Disruptions to left atrial (LA) blood flow, such as those caused by atrial fibrillation (AF), can lead to thrombosis in the left atrial appendage (LAA) and an increased risk of systemic embolism. LA hemodynamics are influenced by various factors, including LA anatomy and function, and pulmonary vein (PV) inflow conditions. In particular, the PV flow split can vary significantly among and within patients depending on multiple factors. In this study, we investigated how changes in PV flow split affect LA flow transport, focusing for the first time on blood stasis in the LAA, using a high-fidelity patient-specific computational fluid dynamics (CFD) model. We use an Immersed Boundary Method, simulating the flow in a fixed, uniform Cartesian mesh and imposing the movement of the LA walls with a moving Lagrangian mesh generated from 4D Computerized Tomography images. We analyzed LA anatomies from eight patients with varying atrial function, including three with AF and either a LAA thrombus or a history of Transient Ischemic Attacks (TIAs). Using four different flow splits (60/40% and 55/45% through right and left PVs, even flow rate, and same velocity through each PV), we found that flow patterns are sensitive to PV flow split variations, particularly in planes parallel to the mitral valve. Changes in PV flow split also had a significant impact on blood stasis and could contribute to increased risk for thrombosis inside the LAA, particularly in patients with AF and previous LAA thrombus or a history of TIAs. Our study highlights the importance of considering patient-specific PV flow split variations when assessing LA hemodynamics and identifying patients at increased risk for thrombosis and stroke. This knowledge is relevant to planning clinical procedures such as AF ablation or the implementation of LAA occluders.This work was partially supported by Comunidad de Madrid (Synergy Grant Y2018/BIO-4858 PREFI-CM), Spanish Research Agency (AEI, grant number PID2019-107279RB-I00), Instituto de Salud Carlos III (grant numbers PI15/02211-ISBITAMI and DTS/1900063-ISBIFLOW), and by the EU-European Regional Development Fund . Funding for open access charge: Universidad de Málaga / CBUA

Magnetar Driven Bubbles and the Origin of Collimated Outflows in Gamma-ray Bursts

We model the interaction between the wind from a newly formed rapidly rotating magnetar and the surrounding supernova shock and host star. The dynamics is modeled using the two-dimensional, axisymmetric thin-shell equations. In the first ~10-100 seconds after core collapse the magnetar inflates a bubble of plasma and magnetic fields behind the supernova shock. The bubble expands asymmetrically because of the pinching effect of the toroidal magnetic field, just as in the analogous problem of the evolution of pulsar wind nebulae. The degree of asymmetry depends on E_mag/E_tot. The correct value of E_mag/E_tot is uncertain because of uncertainties in the conversion of magnetic energy into kinetic energy at large radii in relativistic winds; we argue, however, that bubbles inflated by newly formed magnetars are likely to be significantly more magnetized than their pulsar counterparts. We show that for a ratio of magnetic to total power supplied by the central magnetar L_mag/L_tot ~ 0.1 the bubble expands relatively spherically. For L_mag/L_tot ~ 0.3, however, most of the pressure in the bubble is exerted close to the rotation axis, driving a collimated outflow out through the host star. This can account for the collimation inferred from observations of long-duration gamma-ray bursts (GRBs). Outflows from magnetars become increasingly magnetically dominated at late times, due to the decrease in neutrino-driven mass loss as the young neutron star cools. We thus suggest that the magnetar-driven bubble initially expands relatively spherically, enhancing the energy of the associated supernova, while at late times it becomes progressively more collimated, producing the GRB.Comment: 14 pages, 8 figures, accepted for publication in MNRA