Temperature-extended Jarzynski relation: Application to the numerical calculation of the surface tension

We consider a generalization of the Jarzynski relation to the case where the system interacts with a bath for which the temperature is not kept constant but can vary during the transformation. We suggest to use this relation as a replacement to the thermodynamic perturbation method or the Bennett method for the estimation of the order-order surface tension by Monte Carlo simulations. To demonstrate the feasibility of the method, we present some numerical data for the 3D Ising model

On the center of mass of Ising vectors

We show that the center of mass of Ising vectors that obey some simple constraints, is again an Ising vector.Comment: 8 pages, 3 figures, LaTeX; Claims in connection with disordered systems have been withdrawn; More detailed description of the simulations; Inset added to figure

Geodesics for Efficient Creation and Propagation of Order along Ising Spin Chains

Experiments in coherent nuclear and electron magnetic resonance, and optical spectroscopy correspond to control of quantum mechanical ensembles, guiding them from initial to final target states by unitary transformations. The control inputs (pulse sequences) that accomplish these unitary transformations should take as little time as possible so as to minimize the effects of relaxation and decoherence and to optimize the sensitivity of the experiments. Here we give efficient syntheses of various unitary transformations on Ising spin chains of arbitrary length. The efficient realization of the unitary transformations presented here is obtained by computing geodesics on a sphere under a special metric. We show that contrary to the conventional belief, it is possible to propagate a spin order along an Ising spin chain with coupling strength J (in units of Hz), significantly faster than 1/(2J) per step. The methods presented here are expected to be useful for immediate and future applications involving control of spin dynamics in coherent spectroscopy and quantum information processing

Exact sampling from non-attractive distributions using summary states

Propp and Wilson's method of coupling from the past allows one to efficiently generate exact samples from attractive statistical distributions (e.g., the ferromagnetic Ising model). This method may be generalized to non-attractive distributions by the use of summary states, as first described by Huber. Using this method, we present exact samples from a frustrated antiferromagnetic triangular Ising model and the antiferromagnetic q=3 Potts model. We discuss the advantages and limitations of the method of summary states for practical sampling, paying particular attention to the slowing down of the algorithm at low temperature. In particular, we show that such a slowing down can occur in the absence of a physical phase transition.Comment: 5 pages, 6 EPS figures, REVTeX; additional information at http://wol.ra.phy.cam.ac.uk/mackay/exac

Liquid antiferromagnets in two dimensions

It is shown that, for proper symmetry of the parent lattice, antiferromagnetic order can survive in two-dimensional liquid crystals and even isotropic liquids of point-like particles, in contradiction to what common sense might suggest. We discuss the requirements for antiferromagnetic order in the absence of translational and/or orientational lattice order. One example is the honeycomb lattice, which upon melting can form a liquid crystal with quasi-long-range orientational and antiferromagnetic order but short-range translational order. The critical properties of such systems are discussed. Finally, we draw conjectures for the three-dimensional case.Comment: 4 pages RevTeX, 4 figures include

Time-resolved optical spectroscopy of the pulsating DA white dwarf HS 0507+0434B: New constraints on mode identification and pulsation properties

We present a detailed analysis of time-resolved optical spectra of the ZZ Ceti white dwarf, HS 0507+0434B. Using the wavelength dependence of observed mode amplitudes, we deduce the spherical degree, l, of the modes, most of which have l=1. The presence of a large number of combination frequencies (linear sums or differences of the real modes) enabled us not only to test theoretical predictions but also to indirectly infer spherical and azimuthal degrees of real modes that had no observed splittings. In addition to the above, we measure line-of-sight velocities from our spectra. We find only marginal evidence for periodic modulation associated with the pulsation modes: at the frequency of the strongest mode in the lightcurve, we measure an amplitude of 2.6+/-1.0 km/s, which has a probability of 2% of being due to chance; for the other modes, we find lower values. Our velocity amplitudes and upper limits are smaller by a factor of two compared to the amplitudes found in ZZ Psc. We find that this is consistent with expectations based on the position of HS 0507+0434B in the instability strip. Combining all the available information from data such as ours is a first step towards constraining atmospheric properties in a convectionally unstable environment from an observational perspective.Comment: 16 pages, 12 figs.; accepted for publication in A&

The Yang-Lee zeros of the 1D Blume-Capel model on connected and non-connected rings

We carry out a numerical and analytic analysis of the Yang-Lee zeros of the 1D Blume-Capel model with periodic boundary conditions and its generalization on Feynman diagrams for which we include sums over all connected and non-connected rings for a given number of spins. In both cases, for a specific range of the parameters, the zeros originally on the unit circle are shown to departure from it as we increase the temperature beyond some limit. The curve of zeros can bifurcate and become two disjoint arcs as in the 2D case. We also show that in the thermodynamic limit the zeros of both Blume-Capel models on the static (connected ring) and on the dynamical (Feynman diagrams) lattice tend to overlap. In the special case of the 1D Ising model on Feynman diagrams we can prove for arbitrary number of spins that the Yang-Lee zeros must be on the unit circle. The proof is based on a property of the zeros of Legendre Polynomials.Comment: 19 pages, 5 figure

Deep Markov Random Field for Image Modeling

Markov Random Fields (MRFs), a formulation widely used in generative image modeling, have long been plagued by the lack of expressive power. This issue is primarily due to the fact that conventional MRFs formulations tend to use simplistic factors to capture local patterns. In this paper, we move beyond such limitations, and propose a novel MRF model that uses fully-connected neurons to express the complex interactions among pixels. Through theoretical analysis, we reveal an inherent connection between this model and recurrent neural networks, and thereon derive an approximated feed-forward network that couples multiple RNNs along opposite directions. This formulation combines the expressive power of deep neural networks and the cyclic dependency structure of MRF in a unified model, bringing the modeling capability to a new level. The feed-forward approximation also allows it to be efficiently learned from data. Experimental results on a variety of low-level vision tasks show notable improvement over state-of-the-arts.Comment: Accepted at ECCV 201

Cavopulmonary assist for the failing Fontan circulation: impact of ventricular function on mechanical support strategy

Mechanical circulatory support--either ventricular assist device (VAD, left-sided systemic support) or cavopulmonary assist device (CPAD, right-sided support)--has been suggested as treatment for Fontan failure. The selection of left- versus right-sided support for failing Fontan has not been previously defined. Computer simulation and mock circulation models of pediatric Fontan patients (15-25 kg) with diastolic, systolic, and combined systolic and diastolic dysfunction were developed. The global circulatory response to assisted Fontan flow using VAD (HeartWare HVAD, Miami Lakes, FL) support, CPAD (Viscous Impeller Pump, Indianapolis, IN) support, and combined VAD and CPAD support was evaluated. Cavopulmonary assist improves failing Fontan circulation during diastolic dysfunction but preserved systolic function. In the presence of systolic dysfunction and elevated ventricular end-diastolic pressure (VEDP), VAD support augments cardiac output and diminishes VEDP, while increased preload with cavopulmonary assist may worsen circulatory status. Fontan circulation can be stabilized to biventricular values with modest cavopulmonary assist during diastolic dysfunction. Systemic VAD support may be preferable to maintain systemic output during systolic dysfunction. Both systemic and cavopulmonary support may provide best outcome during combined systolic and diastolic dysfunction. These findings may be useful to guide clinical cavopulmonary assist strategies in failing Fontan circulations

Probability-Changing Cluster Algorithm: Study of Three-Dimensional Ising Model and Percolation Problem

We present a detailed description of the idea and procedure for the newly proposed Monte Carlo algorithm of tuning the critical point automatically, which is called the probability-changing cluster (PCC) algorithm [Y. Tomita and Y. Okabe, Phys. Rev. Lett. {\bf 86} (2001) 572]. Using the PCC algorithm, we investigate the three-dimensional Ising model and the bond percolation problem. We employ a refined finite-size scaling analysis to make estimates of critical point and exponents. With much less efforts, we obtain the results which are consistent with the previous calculations. We argue several directions for the application of the PCC algorithm.Comment: 6 pages including 8 eps figures, to appear in J. Phys. Soc. Jp