Variational quantum Monte Carlo simulations with tensor-network states

We show that the formalism of tensor-network states, such as the matrix product states (MPS), can be used as a basis for variational quantum Monte Carlo simulations. Using a stochastic optimization method, we demonstrate the potential of this approach by explicit MPS calculations for the transverse Ising chain with up to N=256 spins at criticality, using periodic boundary conditions and D*D matrices with D up to 48. The computational cost of our scheme formally scales as ND^3, whereas standard MPS approaches and the related density matrix renromalization group method scale as ND^5 and ND^6, respectively, for periodic systems.Comment: 4+ pages, 2 figures. v2: improved data, comparisons with exact results, to appear in Phys Rev Let

LENS® and SFF: Enabling Technologies for Optimized Structures

Optimized, lightweight, high-strength structures are needed in many applications from aerospace to automotive. In pursuit of such structures, there have been proposed analytical solutions and some specialized FEA solutions for specific structures such as automobile frames. However, generalized 3D optimization methods have been unavailable for use by most designers. Moreover, in the cases where optimized structural solutions are available, they are often hollow, curving, thin wall structures that cannot be fabricated by conventional manufacturing methods. Researchers at Sandia National Laboratories and the University of Rhode Island teamed to solve these problems. The team has been pursuing two methods of optimizing models for generalized loading conditions, and also has been investigating the methods needed to fabricate these structures using Laser Engineered Net Shaping™ (LENS®) and other rapid prototyping methods. These solid freeform fabrication (SFF) methods offer the unique ability to make hollow, high aspect ratio features out of many materials. The manufacturing development required for LENS to make these complex structures has included the addition of rotational axes to Sandia’s LENS machine bringing the total to 5 controlled axes. The additional axes have required new efforts in process planning. Several of the unique structures that are only now possible through the use of SFF technology are shown as part of the discussion of this exciting new application for SFF.Mechanical Engineerin

Work-related psychological health and psychological type among lead elders within the Newfrontiers network of churches in the United Kingdom

Building on a series of recent studies concerned with assessing work-related psychological health and psychological type among various groups of church leaders, this study reports new data provided by 134 Lead Elders within the Newfrontiers network of churches in the United Kingdom who completed the Francis Psychological Type Scales (FPTS) together with the two scales of the Francis Burnout Inventory (FBI) concerned with emotional exhaustion and satisfaction in ministry. Compared with other groups of church leaders, Lead Elders within the Newfrontiers network of churches reported lower levels of emotional exhaustion and higher levels of satisfaction in ministry. Compared with other groups of church leaders, there was a higher proportion of extraverts among Lead Elders within the Newfrontiers network of churches. There was only a weak association between psychological type and burnout

Fcc-bcc transition for Yukawa interactions determined by applied strain deformation

Calculations of the work required to transform between bcc and fcc phases yield a high-precision bcc-fcc transition line for monodisperse point Yukawa (screened-Couloumb) systems. Our results agree qualitatively but not quantitatively with previously published simulations and phenomenological criteria for the bcc-fcc transition. In particular, the bcc-fcc-fluid triple point lies at a higher inverse screening length than previously reported.Comment: RevTex4, 9 pages, 6 figures. Discussion of phase coexistence extended, a few other minor clarifications added, referencing improved. Accepted for publication by Physical Review

Classification of unit-vector fields in convex polyhedra with tangent boundary conditions

A unit-vector field n on a convex three-dimensional polyhedron P is tangent if, on the faces of P, n is tangent to the faces. A homotopy classification of tangent unit-vector fields continuous away from the vertices of P is given. The classification is determined by certain invariants, namely edge orientations (values of n on the edges of P), kink numbers (relative winding numbers of n between edges on the faces of P), and wrapping numbers (relative degrees of n on surfaces separating the vertices of P), which are subject to certain sum rules. Another invariant, the trapped area, is expressed in terms of these. One motivation for this study comes from liquid crystal physics; tangent unit-vector fields describe the orientation of liquid crystals in certain polyhedral cells.Comment: 21 pages, 2 figure

Spin chains and combinatorics: twisted boundary conditions

The finite XXZ Heisenberg spin chain with twisted boundary conditions was considered. For the case of even number of sites $N$, anisotropy parameter -1/2 and twisting angle $2 \pi /3$ the Hamiltonian of the system possesses an eigenvalue $-3N/2$. The explicit form of the corresponding eigenvector was found for $N \le 12$. Conjecturing that this vector is the ground state of the system we made and verified several conjectures related to the norm of the ground state vector, its component with maximal absolute value and some correlation functions, which have combinatorial nature. In particular, the squared norm of the ground state vector is probably coincides with the number of half-turn symmetric alternating sign $N \times N$ matrices.Comment: LaTeX file, 7 page

Dynamic sea surface topography, gravity and improved orbit accuracies from the direct evaluation of SEASAT altimeter data

A method for the simultaneous solution of dynamic ocean topography, gravity and orbits using satellite altimeter data is described. A GEM-T1 based gravitational model called PGS-3337 that incorporates Seasat altimetry, surface gravimetry and satellite tracking data has been determined complete to degree and order 50. The altimeter data is utilized as a dynamic observation of the satellite's height above the sea surface with a degree 10 model of dynamic topography being recovered simultaneously with the orbit parameters, gravity and tidal terms in this model. PGS-3337 has a geoid uncertainty of 60 cm root-mean-square (RMS) globally, with the uncertainty over the altimeter tracked ocean being in the 25 cm range. Doppler determined orbits for Seasat, show large improvements, with the sub-30 cm radial accuracies being achieved. When altimeter data is used in orbit determination, radial orbital accuracies of 20 cm are achieved. The RMS of fit to the altimeter data directly gives 30 cm fits for Seasat when using PGS-3337 and its geoid and dynamic topography model. This performance level is two to three times better than that achieved with earlier Goddard earth models (GEM) using the dynamic topography from long-term oceanographic averages. The recovered dynamic topography reveals the global long wavelength circulation of the oceans with a resolution of 1500 km. The power in the dynamic topography recovery is now found to be closer to that of oceanographic studies than for previous satellite solutions. This is attributed primarily to the improved modeling of the geoid which has occurred. Study of the altimeter residuals reveals regions where tidal models are poor and sea state effects are major limitations

Algorithmic Randomness and Capacity of Closed Sets

We investigate the connection between measure, capacity and algorithmic randomness for the space of closed sets. For any computable measure m, a computable capacity T may be defined by letting T(Q) be the measure of the family of closed sets K which have nonempty intersection with Q. We prove an effective version of Choquet's capacity theorem by showing that every computable capacity may be obtained from a computable measure in this way. We establish conditions on the measure m that characterize when the capacity of an m-random closed set equals zero. This includes new results in classical probability theory as well as results for algorithmic randomness. For certain computable measures, we construct effectively closed sets with positive capacity and with Lebesgue measure zero. We show that for computable measures, a real q is upper semi-computable if and only if there is an effectively closed set with capacity q

Gravitational model improvement at the Goddard Space Flight Center

Major new computations of terrestrial gravitational field models were performed by the Geodynamics Branch of Goddard Space Flight Center (GSFC). This development has incorporated the present state of the art results in satellite geodesy and have relied upon a more consistent set of reference constants than was heretofore utilized in GSFC's GEM models. The solutions are complete in spherical harmonic coefficients out to degree 50 for the gravity field parameters. These models include adjustment for a subset of 66 ocean tidal coefficients for the long wavelength components of 12 major ocean tides. This tidal adjustment was made in the presence of 550 other fixed ocean tidal terms representing 32 major and minor ocean tides and the Wahr frequency dependent solid earth tidal model. In addition 5-day averaged values for Earth rotation and polar motion were derived for the time period of 1980 onward. Two types of models were computed. These are satellite only models relying exclusively on tracking data and combination models which have incorporated satellite altimetry and surface gravity data. The satellite observational data base consists of over 1100 orbital arcs of data on 31 satellites. A large percentage of these observations were provided by third generation laser stations (less than 5 cm). A calibration of the model accuracy of the GEM-T2 satellite only solution indicated that it was a significant improvement over previous models based solely upon tracking data. The rms geoid error for this field is 110 cm to degree and order 36. This is a major advancement over GEM-T1 whose errors were estimated to be 160 cm. An error propagation using the covariances of the GEM-T2 model for the TOPEX radial orbit component indicates that the rms radial errors are expected to be 12 cm. The combination solution, PGS-3337, is a preliminary effort leading to the development of GEM-T3. PGS-3337 has incorporated global sets of surface gravity data and the Seasat altimetry to produce a model complete to (50,50). A solution for the dynamic ocean topography to degree and order 10 was included as part of this adjustment

Evolution of displacements and strains in sheared amorphous solids

The local deformation of two-dimensional Lennard-Jones glasses under imposed shear strain is studied via computer simulations. Both the mean squared displacement and mean squared strain rise linearly with the length of the strain interval $\Delta \gamma$ over which they are measured. However, the increase in displacement does not represent single-particle diffusion. There are long-range spatial correlations in displacement associated with slip lines with an amplitude of order the particle size. Strong dependence on system size is also observed. The probability distributions of displacement and strain are very different. For small $\Delta \gamma$ the distribution of displacement has a plateau followed by an exponential tail. The distribution becomes Gaussian as $\Delta \gamma$ increases to about .03. The strain distributions consist of sharp central peaks associated with elastic regions, and long exponential tails associated with plastic regions. The latter persist to the largest $\Delta \gamma$ studied.Comment: Submitted to J. Phys. Cond. Mat. special volume for PITP Conference on Mechanical Behavior of Glassy Materials. 16 Pages, 8 figure