Bose-Einstein condensation of trapped atoms with dipole interactions

The path integral Monte Carlo method is used to simulate dilute trapped Bose gases and to investigate the equilibrium properties at finite temperatures. The quantum particles have a long-range dipole-dipole interaction and a short-range s-wave interaction. Using an anisotropic pseudopotential for the long-range dipolar interaction and a hard-sphere potential for the short-range s-wave interaction, we calculate the energetics and structural properties as a function of temperature and the number of particles. Also, in order to determine the effects of dipole-dipole forces and the influence of the trapping field on the dipolar condensate, we use two cylindrically symmetric harmonic confinements (a cigar-shaped trap and a disk-shaped trap). We find that the net effect of dipole-dipole interactions is governed by the trapping geometry. For a cigar-shaped trap, the net contribution of dipolar interactions is attractive and the shrinking of the density profiles is observed. For a disk-shaped trap, the net effect of long-range dipolar forces is repulsive and the density profiles expand

Finite-temperature properties of quasi-2D Bose-Einstein condensates

Using the finite-temperature path integral Monte Carlo method, we investigate dilute, trapped Bose gases in a quasi-two dimensional geometry. The quantum particles have short-range, s-wave interactions described by a hard-sphere potential whose core radius equals its corresponding scattering length. The effect of both the temperature and the interparticle interaction on the equilibrium properties such as the total energy, the density profile, and the superfluid fraction is discussed. We compare our accurate results with both the semi-classical approximation and the exact results of an ideal Bose gas. Our results show that for repulsive interactions, (i) the minimum value of the aspect ratio, where the system starts to behave quasi-two dimensionally, increases as the two-body interaction strength increases, (ii) the superfluid fraction for a quasi-2D Bose gas is distinctly different from that for both a quasi-1D Bose gas and a true 3D system, i.e., the superfluid fraction for a quasi-2D Bose gas decreases faster than that for a quasi-1D system and a true 3D system with increasing temperature, and shows a stronger dependence on the interaction strength, (iii) the superfluid fraction for a quasi-2D Bose gas lies well below the values calculated from the semi-classical approximation, and (iv) the Kosterlitz-Thouless transition temperature decreases as the strength of the interaction increases.Comment: 6 pages, 5 figure

Bose-Einstein Condensation Temperature of a Homogeneous Weakly Interacting Bose Gas : PIMC study

Using a finite-temperature Path Integral Monte Carlo simulation (PIMC) method and finite-size scaling, we have investigated the interaction-induced shift of the phase transition temperature for Bose-Einstein condensation of homogeneous weakly interacting Bose gases in three dimensions, which is given by a proposed analytical expression $T_{c} = T_{c}^{0}\{1 + c_{1}an^{1/3}+[c'_{2}\ln(an^{1/3})+c''_{2}]a^{2}n^{2/3} +O(a^{3}n)\}$, where $T_{c}^{0}$ is the critical temperature for an ideal gas, $a$ is the s-wave scattering length, and $n$ is the number density. We have used smaller number densities and more time slices than in the previous PIMC simulations [Gruter {\it et al.}, Phys. Rev. Lett. {\bf 79}, 3549 (1997)] in order to understand the difference in the value of the coefficient $c_{1}$ between their results and the (apparently) other reliable results in the literature. Our results show that $\{(T_{c}-T_{c}^{0})/T_{c}^{0}\}/(an^{1/3})$ depends strongly on the interaction strength $an^{1/3}$ while the previous PIMC results are considerably flatter and smaller than our results. We obtain $c_{1}$ = 1.32 $\pm$ 0.14, in agreement with results from recent Monte Carlo methods of three-dimensional O(2) scalar $\phi^{4}$ field theory and variational perturbation theory

Self-reported Mental Disorders Negatively Influence Surgical Outcomes After Arthroscopic Treatment of Femoroacetabular Impingement.

Background:Femoroacetabular impingement (FAI) is responsible for hip pain and dysfunction, and surgical outcomes depend on multiple factors. The presence of mental disorders negatively influences outcomes of multiple orthopaedic conditions, although the impact on FAI surgery is unclear. Hypothesis:The authors hypothesized that a preoperative self-reported history of mental disorders would negatively influence patient-reported outcome measures after FAI surgery. Study Design:Cohort study; Level of evidence, 3. Methods:A matched-cohort study was performed by reviewing a prospectively collected database of cases of arthroscopic management of FAI with a single surgeon over a 2-year period. Demographics and radiographic parameters were recorded for all patients. Patients completed the Hip Outcome Score-Activity of Daily Living Subscale (HOS-ADL), Hip Outcome Score-Sport-Specific Subscale (HOS-SSS), and modified Harris Hip Score (mHHS) prior to surgery and 2 years after surgery. Unpaired and paired t tests were used to compare results between and within cohorts at baseline and follow-up. Statistical significance was defined as P < .05. Results:The cohort included 301 patients, with 75 and 226 patients reporting and not reporting a history of mental disorders, respectively. Before treatment, all patient-reported outcome measures were significantly lower among patients reporting a history of mental disorders (P < .01 for HOS-ADL, HOS-SSS, and mHHS). Patients in both groups demonstrated significant improvements (P < .0001) in HOS-ADL, HOS-SSS, and mHHS when preoperative outcome measures were compared with follow-up. Patients with reported mental disorders had significantly lower scores after surgery as compared with patients without mental disorders (P < .0001 for HOS-ADL, HOS-SSS, and mHHS). Conclusion:The presence of a reported mental disorder is associated with lower patient-reported outcomes before and after surgical management of FAI. Statistically significant and clinically relevant improvements were observed for patients who reported mental disorders. The magnitude of these improvements was not as large as that for an age- and sex-matched control group without a self-reported mental disorder

Adsorption of para-Hydrogen on Krypton pre-plated graphite

Adsorption of para-Hydrogen on the surface of graphite pre-plated with a single layer of atomic krypton is studied thoretically by means of Path Integral Ground State Monte Carlo simulations. We compute energetics and density profiles of para-hydrogen, and determine the structure of the adsorbed film for various coverages. Results show that there are two thermodynamically stable monolayer phases of para-hydrogen, both solid. One is commensurate with the krypton layer, the other is incommensurate. No evidence is seen of a thermodynamically stable liquid phase, at zero temperature. These results are qualitatively similar to what is seen for for para-hydrogen on bare graphite. Quantum exchanges of hydrogen molecules are suppressed in this system.Comment: 12 pages, 6 figures, to appear in the proceedings of "Advances in Computational Many-Body Physics", Banff, Alberta (Canada), January 13-16 200

Prevalence of Cam Morphology in Females with Femoroacetabular Impingement

Cam and pincer are two common morphologies responsible for femoroacetabular impingement. Previous literature has reported that cam deformity is predominantly a male morphology, while being significantly less common in females. The purpose of this study was to determine the prevalence of cam morphology in female subjects diagnosed with symptomatic FAI. All females presenting to the senior author’s clinic diagnosed with symptomatic FAI between December 2006 and Cam and pincer are two common morphologies responsible for femoroacetabular impingement. Previous literature has reported that cam deformity is predominantly a male morphology, while being significantly less common in females. Cam morphology is commonly assessed with the alpha angle, measured on radiographs. The purpose of this study was to determine the prevalence of cam morphology utilizing the alpha angle in female subjects diagnosed with symptomatic FAI. All females presenting to the senior author’s clinic diagnosed with symptomatic FAI between December 2006 and January 2013 were retrospectively reviewed. Alpha (α) angles were measured on AP (anteroposterior) and lateral (Dunn 90°, cross-table lateral, and/or frog-leg lateral) plain radiographs by two blinded physicians, and the largest measured angle was used. Using Gosvig et al.’s classification, alpha angle was characterized as (pathologic > 57°), borderline (51-56°), subtle (46-50°), very subtle (43-45°), or normal (≤42°). Three hundred and ninety-one patients (438 hips) were analyzed (age 36.2 ± 12.3 years). Among the hips included, 35.6% were normal, 14.6% pathologic, 15.1% borderline, 14.6% subtle, and 20.1% very subtle. There was no correlation between alpha angle and patient age (R = 0.17) or body mass index (BMI) (R = 0.05). The intraclass correlation coefficient (ICC) for α-angle measurements was 0.84. Sixty-four percent of females in this cohort had an alpha angle > 42°. Subtle cam deformity plays a significant role in the pathoanatomy of female patients with symptomatic FAI. As the majority of revision hip arthroscopies are performed due to incomplete cam correction, hip arthroscopists need to be cognizant of and potentially surgically address these subtle lesions

Identification of nonlinear heat transfer laws from boundary observations

We consider the problem of identifying a nonlinear heat transfer law at the boundary, or of the temperature-dependent heat transfer coefficient in a parabolic equation from boundary observations. As a practical example, this model applies to the heat transfer coefficient that describes the intensity of heat exchange between a hot wire and the cooling water in which it is placed. We reformulate the inverse problem as a variational one which aims to minimize a misfit functional and prove that it has a solution. We provide a gradient formula for the misfit functional and then use some iterative methods for solving the variational problem. Thorough investigations are made with respect to several initial guesses and amounts of noise in the input data. Numerical results show that the methods are robust, stable and accurate

Spin-dynamics simulations of the triangular antiferromagnetic XY model

Using Monte Carlo and spin-dynamics methods, we have investigated the dynamic behavior of the classical, antiferromagnetic XY model on a triangular lattice with linear sizes $L \leq 300$. The temporal evolutions of spin configurations were obtained by solving numerically the coupled equations of motion for each spin using fourth-order Suzuki-Trotter decompositions of exponential operators. From space- and time-displaced spin-spin correlation functions and their space-time Fourier transforms we obtained the dynamic structure factor $S({\bf q},w)$ for momentum ${\bf q}$ and frequency $\omega$. Below $T_{KT}$(Kosterlitz-Thouless transition), both the in-plane ($S^{xx}$) and the out-of-plane ($S^{zz}$) components of $S({\bf q},\omega)$ exhibit very strong and sharp spin-wave peaks. Well above $T_{KT}$, $S^{xx}$ and $S^{zz}$ apparently display a central peak, and spin-wave signatures are still seen in $S^{zz}$. In addition, we also observed an almost dispersionless domain-wall peak at high $\omega$ below $T_{c}$(Ising transition), where long-range order appears in the staggered chirality. Above $T_{c}$, the domain-wall peak disappears for all $q$. The lineshape of these peaks is captured reasonably well by a Lorentzian form. Using a dynamic finite-size scaling theory, we determined the dynamic critical exponent $z$ = 1.002(3). We found that our results demonstrate the consistency of the dynamic finite-size scaling theory for the characteristic frequeny $\omega_{m}$ and the dynamic structure factor $S({\bf q},\omega)$ itself.Comment: 8 pages, RevTex, 10 figures, submitted to PR

Comparison of Multi-Sample Variant Calling Methods for Whole Genome Sequencing

Rapid advancement of next-generation sequencing (NGS) technologies has facilitated the search for genetic susceptibility factors that influence disease risk in the field of human genetics. In particular whole genome sequencing (WGS) has been used to obtain the most comprehensive genetic variation of an individual and perform detailed evaluation of all genetic variation. To this end, sophisticated methods to accurately call high-quality variants and genotypes simultaneously on a cohort of individuals from raw sequence data are required. On chromosome 22 of 818 WGS data from the Alzheimer's Disease Neuroimaging Initiative (ADNI), which is the largest WGS related to a single disease, we compared two multi-sample variant calling methods for the detection of single nucleotide variants (SNVs) and short insertions and deletions (indels) in WGS: (1) reduce the analysis-ready reads (BAM) file to a manageable size by keeping only essential information for variant calling ("REDUCE") and (2) call variants individually on each sample and then perform a joint genotyping analysis of the variant files produced for all samples in a cohort ("JOINT"). JOINT identified 515,210 SNVs and 60,042 indels, while REDUCE identified 358,303 SNVs and 52,855 indels. JOINT identified many more SNVs and indels compared to REDUCE. Both methods had concordance rate of 99.60% for SNVs and 99.06% for indels. For SNVs, evaluation with HumanOmni 2.5M genotyping arrays revealed a concordance rate of 99.68% for JOINT and 99.50% for REDUCE. REDUCE needed more computational time and memory compared to JOINT. Our findings indicate that the multi-sample variant calling method using the JOINT process is a promising strategy for the variant detection, which should facilitate our understanding of the underlying pathogenesis of human diseases