Uncertainty in cardiac myofiber orientation and stiffnesses dominate the variability of left ventricle deformation response.

Computational cardiac modelling is a mature area of biomedical computing and is currently evolving from a pure research tool to aiding in clinical decision making. Assessing the reliability of computational model predictions is a key factor for clinical use, and uncertainty quantification (UQ) and sensitivity analysis are important parts of such an assessment. In this study, we apply UQ in computational heart mechanics to study uncertainty both in material parameters characterizing global myocardial stiffness and in the local muscle fiber orientation that governs tissue anisotropy. The uncertainty analysis is performed using the polynomial chaos expansion (PCE) method, which is a nonintrusive meta-modeling technique that surrogates the original computational model with a series of orthonormal polynomials over the random input parameter space. In addition, in order to study variability in the muscle fiber architecture, we model the uncertainty in orientation of the fiber field as an approximated random field using a truncated Karhunen-Loéve expansion. The results from the UQ and sensitivity analysis identify clear differences in the impact of various material parameters on global output quantities. Furthermore, our analysis of random field variations in the fiber architecture demonstrate a substantial impact of fiber angle variations on the selected outputs, highlighting the need for accurate assignment of fiber orientation in computational heart mechanics models.

Comparative investigation of pure and mixed rare gas atoms on coronene molecules.

Clusters formed by the combination of rare gas (RG) atoms of He, Ne, Ar, and Kr on coronene have been investigated by means of a basin-hopping algorithm and path integral Monte Carlo calculations at T = 2 K. Energies and geometries have been obtained and the role played by the specific RG-RG and RG-coronene interactions on the final results is analysed in detail. Signatures of diffuse behavior of the He atoms on the surface of the coronene are in contrast with the localization of the heavier species, Ar and Kr. The observed coexistence of various geometries for Ne suggests the motion of the RG atoms on the multi-well potential energy surface landscape offered by the coronene. Therefore, the investigation of different clusters enables a comparative analysis of localized versus non-localized features. Mixed Ar-He-coronene clusters have also been considered and the competition of the RG atoms to occupy the docking sites on the molecule is discussed. All the obtained information is crucial to assess the behavior of coronene, a prototypical polycyclic aromatic hydrocarbon clustering with RG atoms at a temperature close to that of interstellar medium, which arises from the critical balance of the interactions involved.

Examination of the Feynman-Hibbs Approach in the Study of NeN-Coronene Clusters at Low Temperatures.

Feynman-Hibbs (FH) effective potentials constitute an appealing approach for investigations of many-body systems at thermal equilibrium since they allow us to easily include quantum corrections within standard classical simulations. In this work we apply the FH formulation to the study of NeN-coronene clusters (N = 1-4, 14) in the 2-14 K temperature range. Quadratic (FH2) and quartic (FH4) contributions to the effective potentials are built upon Ne-Ne and Ne-coronene analytical potentials. In particular, a new corrected expression for the FH4 effective potential is reported. FH2 and FH4 cluster energies and structures-obtained from energy optimization through a basin-hopping algorithm as well as classical Monte Carlo simulations-are reported and compared with reference path integral Monte Carlo calculations. For temperatures T > 4 K, both FH2 and FH4 potentials are able to correct the purely classical calculations in a consistent way. However, the FH approach fails at lower temperatures, especially the quartic correction. It is thus crucial to assess the range of applicability of this formulation and, in particular, to apply the FH4 potentials with great caution. A simple model of N isotropic harmonic oscillators allows us to propose a means of estimating the cutoff temperature for the validity of the method, which is found to increase with the number of atoms adsorbed on the coronene molecule.

Coronene molecules in helium clusters: Quantum and classical studies of energies and configurations.

Coronene-doped helium clusters have been studied by means of classical and quantum mechanical (QM) methods using a recently developed He-C24H12 global potential based on the use of optimized atom-bond improved Lennard-Jones functions. Equilibrium energies and geometries at global and local minima for systems with up to 69 He atoms were calculated by means of an evolutive algorithm and a basin-hopping approach and compared with results from path integral Monte Carlo (PIMC) calculations at 2 K. A detailed analysis performed for the smallest sizes shows that the precise localization of the He atoms forming the first solvation layer over the molecular substrate is affected by differences between relative potential minima. The comparison of the PIMC results with the predictions from the classical approaches and with diffusion Monte Carlo results allows to examine the importance of both the QM and thermal effects.

Quantum Features of Anionic Species He*⁻ and He2*⁻ in Small He(N) Clusters.

We present variational calculations on systems containing a few boson helium atoms attached to electronically excited atomic and molecular helium anions He*⁻ and He2*⁻ and characterize their structures and energetics. Previously reported high-level ab initio results [Huber, S. E.; Mauracher, A. Mol. Phys. 2014, 112, 794] to describe the interactions between excited (metastable) anions and a neutral He atom have been employed. For the case of the atomic species He*⁻, the corresponding interaction with He suggests large anharmonicity effects due to the presence of a deep well of ∼17,500 cm⁻¹ at short distances, together with a more external shallow secondary well of ∼4 cm⁻¹, both supporting bound levels. Moreover, when a sum of pairwise interactions is assumed to describe the full PES corresponding to the presence of several neutral He atoms, geometrical constraints already predict the complete solvation of the anionic impurity by six helium atoms, giving rise to a bipyramidal structure. In turn, for the anisotropic weak interaction He-He2*⁻, where the anionic dimer is considered as a rigid rotor, the obtained structures show the tendency of the helium atoms to pack themselves together and largely far away from the dopant, thereby confirming the heliophobic character of He2*⁻.

Reactive scattering calculations for (87)Rb+(87)RbHe→Rb2((3)Σ(u)(+),v)+He from ultralow to intermediate energies.

We investigate atom-diatom reactive collisions, as a preliminary step,in order to assess the possibility of forming Rb(2) molecules in their lowest triplet electronic state by cold collisions of rubidium atoms on the surface of helium nanodroplets [corrected]. A simple model related to the well-known Rosen treatment of linear triatomic molecules [N. Rosen, J. Chem. Phys. 1, 319 (1933)] in relative coordinates is used, allowing to estimate reactive probabilities for different values of the total angular momentum. The best available full dimensional potential energy surface [Guillon et al., J. Chem. Phys. 136, 174307 (2012)] is employed through the calculations. Noticeable values of the probabilities in the ultracold regime, which numerically fulfill the Wigner threshold law, support the feasibility of the process. The rubidium dimer is mainly produced at high vibrational states, and the reactivity is more efficient for a bosonic helium partner than when the fermion species is considered.

A configurational study of helium clusters doped with He(∗-) and He2(∗-).

Helium clusters doped with electronically excited atomic and molecular helium anions He(∗-) and He2(∗-) at T = 0.4 K are studied by means of path integral Monte Carlo calculations. Geometry and energetics of the systems with up to 32 solvating He atoms are characterised. The interactions between the anions and the neutral He atoms have been described by fitting previously reported ab initio points to analytical expressions. The HeN-He(∗-) clusters with N > 6 display a structure defined by a bipyramid which completely solvates the atomic anion, whereas the rest of surrounding He atoms form a dimple around that initial cage. On the contrary, the structures observed for the HeN-He2(∗-) clusters clearly show the dopant located outside the helium droplet, thereby confirming the heliophobic character of He2(∗-).