Metropolis Monte Carlo on the Lefschetz thimble: application to a one-plaquette model

We propose a new algorithm based on the Metropolis sampling method to perform Monte Carlo integration for path integrals in the recently proposed formulation of quantum field theories on the Lefschetz thimble. The algorithm is based on a mapping between the curved manifold defined by the Lefschetz thimble of the full action and the flat manifold associated with the corresponding quadratic action. We discuss an explicit method to calculate the residual phase due to the curvature of the Lefschetz thimble. Finally, we apply this new algorithm to a simple one-plaquette model where our results are in perfect agreement with the analytic integration. We also show that for this system the residual phase does not represent a sign problem

Waste Heat Recovery in Food and Drinks Industry (Abstract only)

Most baking processes in the food manufacturing sector involve use of gas-fired ovens. Only about one-third of the total energy used in these ovens adds value to the final product. The remaining two-thirds is discharged with the exhaust gases at 150-250o C and thus represents an opportunity for heat recovery. However, the low temperature range, fouling and presence of corrosive materials in the exhaust streams make heat recovery technically challenging and uneconomical. The existing low grade heat recovery technolgies mostly use gas to liquid heat transfer to produce hot water for use in other areas of the manufacturing plant. The performance of these systems is governed by hot water demand in the factory and is therefore not recommended if there are frequent fluctuations in demand or if a more efficient technology, such as combined heat and power, is already in place. This study involves design, manufacturing and testing of a novel low-temperature gas to gas heat recovery system using an array of heat pipe heat exchangers, for industrial-scale baking ovens at a large confectionary manufacturing plant. Unlike gas to liquid heat transfer, a gas to gas heat transfer system provides direct savings in oven fuel consumption, independent of the hot water and other energy demands elsewhere in the plant. The heat recovery potential of the system is estimated using a thermodynamic model developed based on energy and mass balance for the ovens. The design enables recovery of up to 50% of the energy available through the exhaust stack, increasing the energy efficiency of the overall process to 60% and reducing food manufacturing costs by one third

Quantum Monte Carlo with Coupled-Cluster wave functions

We introduce a novel many body method which combines two powerful many body techniques, viz., quantum Monte Carlo and coupled cluster theory. Coupled cluster wave functions are introduced as importance functions in a Monte Carlo method designed for the configuration interaction framework to provide rigorous upper bounds to the ground state energy. We benchmark our method on the homogeneous electron gas in momentum space. The importance function used is the coupled cluster doubles wave function. We show that the computational resources required in our method scale polynomially with system size. Our energy upper bounds are in very good agreement with previous calculations of similar accuracy, and they can be systematically improved by including higher order excitations in the coupled cluster wave function.Comment: Submitted to Physical Review Letter

Odd-particle systems in the shell model Monte Carlo: circumventing a sign problem

We introduce a novel method within the shell model Monte Carlo approach to calculate the ground-state energy of a finite-size system with an odd number of particles by using the asymptotic behavior of the imaginary-time single-particle Green's functions. The method circumvents the sign problem that originates from the projection on an odd number of particles and has hampered direct application of the shell model Monte Carlo method to odd-particle systems. We apply this method to calculate pairing gaps of nuclei in the iron region. Our results are in good agreement with experimental pairing gaps

Importance of vegetation in tsunami mitigation: evidence from large eddy simulations with fluid-structure interactions

Communities worldwide are increasingly interested in nature-based solutions like coastal forests for the mitigation of coastal risks. Still, it remains unclear how much protective benefit vegetation provides, particularly in the limit of highly energetic flows after tsunami impact. The present thesis, using a three-dimensional incompressible computational fluid dynamics model with a fluid-structure interaction approach, aims to quantify how energy reflection and dissipation vary with different degrees of rigidity and vegetation density of a coastal forest. In this study, tree trunks are represented as cylinders, and the elastic modulus of hardwood trees such as pine or oak is used to characterize the rigidity of these cylinders. To capture tsunami bore propagation in onshore, dam break flow is used over the wet surface in the numerical studies. After validating numerical code against experimental studies, multi-cylinder configurations are incorporated and Froude Number is used to scale the flow parameters and vegetation flow parameter (VFP) to scale the tree parameters such as elastic modulus, the diameter of the trunk, etc. Numerical tests are conducted for different cylinder diameters, densities, and elastic moduli. The numerical results show that energy reflection increases with rigidity only for a single cylinder. In the presence of multiple cylinders, the difference in energy reflection created by varying rigidity diminishes as the number of cylinders increases. Instead of rigidity, the blockage area created by the presence of multiple tree trunks is found to dominate energy reflection. As tree trunks are deformed by the hydrodynamic forces, they alter the flow field around them, causing turbulent kinetic energy generation in the wake region. As a consequence, trees dissipate flow energy, highlighting the importance of coastal forests in reducing the onshore energy flux of tsunamis by means of both reflection and dissipation

Number-conserving theory of nuclear pairing gaps: a global assessment

We study odd-even mass staggering of nuclei, also called pairing gaps, using a Skyrme self-consistent mean-field theory and a numerically exact treatment of the pairing Hamiltonian. We find that the configuration-space Monte Carlo method proposed by Cerf and Martin offers a practical computational procedure to carry out the numerical solutions in large-dimensional model spaces. Refitting the global strength of the pairing interaction for 443 neutron pairing gaps in our number-conserving treatment, we find the correction to the pairing correlation energies and pairing gaps to have rms values of 0.6 MeV and 0.12 MeV, respectively. The exact treatment provides a significant improvement in the fit to experimental gaps, although it is partially masked by a larger rms error due to deficiencies in other aspects of the theory such as the mean-field energy functional.Comment: 11 pages, 9 figure