Source author record

Abhishek Mukherjee

Abhishek Mukherjee appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Forest density is more effective than tree rigidity at reducing the onshore energy flux of tsunamis: Evidence from Large Eddy Simulations with Fluid-Structure Interactions

Communities around the world are increasingly interested in nature-based solutions to mitigation of coastal risks like coastal forests, but it remains unclear how much protective benefits vegetation provides, particularly in the limit of highly energetic flows after tsunami impact. The current study, 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. We represent tree trunks as cylinders and use the elastic modulus of hardwood trees such as pine or oak to characterize the rigidity of these cylinders. 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, we find that the blockage area created by the presence of multiple tree trunks dominates 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 coastal forests reducing the onshore energy flux of tsunamis by means of both reflection and dissipation.

preprint2015arXiv

Constraining the nuclear energy density functional with quantum Monte Carlo calculations

We study the problem of an impurity in fully polarized (spin-up) low density neutron matter with the help of an accurate quantum Monte Carlo method in conjunction with a realistic nucleon-nucleon interaction derived from chiral effective field theory at next-to-next-to-leading-order. Our calculations show that the behavior of the proton spin-down impurity is very similar to that of a polaron in a fully polarized unitary Fermi gas. We show that our results can be used to put tight constraints on the time-odd parts of the energy density functional, independent of the time-even parts, in the density regime relevant to neutron-rich nuclei and compact astrophysical objects such as neutron stars and supernovae.

preprint2014arXiv

Lefschetz thimble Monte Carlo for many body theories: application to the repulsive Hubbard model away from half filling

Recently, a new method, based on stochastic integration on the surfaces of steepest descent of the action, was introduced to tackle the sign problem in quantum field theories. We show how this method can be used in many body theories to perform fully non-perturbative calculations of quantum corrections about mean field solutions. We discuss an explicit algorithm for implementing our method, and present results for the repulsive Hubbard model away from half-filling at intermediate temperatures. Our results are consistent with those from the recent state of the art cluster dynamical mean field theory calculations.

preprint2014arXiv

Quantum Monte Carlo calculations of neutron matter with non-local chiral interactions

We present fully non-perturbative quantum Monte Carlo calculations with non-local chiral effective field theory (EFT) interactions for the ground state properties of neutron matter. The equation of state, the nucleon chemical potentials and the momentum distribution in pure neutron matter up to one and a half times the nuclear saturation density are computed with a newly optimized chiral EFT interaction at next-to-next-to-leading order. This work opens the way to systematic order by order benchmarking of chiral EFT interactions, and \emph{ab initio} prediction of nuclear properties while respecting the symmetries of quantum chromodynamics.

preprint2013arXiv

Configuration-interaction Monte Carlo method and its application to the trapped unitary Fermi gas

We develop a quantum Monte Carlo method to estimate the ground-state energy of a fermionic many-particle system in the configuration-interaction shell model approach. The fermionic sign problem is circumvented by using a guiding wave function in Fock space. The method provides an upper bound on the ground-state energy whose tightness depends on the choice of the guiding wave function. We argue that the antisymmetric geminal product class of wave functions is a good choice for guiding wave functions. We demonstrate our method for the trapped two-species fermionic cold atom system in the unitary regime of infinite scattering length using the particle-number projected Hartree-Fock-Bogoliubov wave function as the guiding wave function. We estimate the ground-state energy and energy-staggering pairing gap as a function of the number of particles. Our results compare favorably with exact numerical diagonalization results and with previous coordinate-space Monte Carlo calculations.

preprint2013arXiv

Level densities of nickel isotopes: microscopic theory versus experiment

We apply a spin-projection method to calculate microscopically the level densities of a family of nickel isotopes $^{59-64}$Ni using the shell model Monte Carlo approach in the complete $pfg_{9/2}$ shell. Accurate ground-state energies of the odd-mass nickel isotopes, required for the determination of excitation energies, are determined using the Green's function method recently introduced to circumvent the odd particle-number sign problem. Our results are in excellent agreement with recent measurements based on proton evaporation spectra and with level counting data at low excitation energies. We also compare our results with neutron resonance data, assuming equilibration of parity and a spin-cutoff model for the spin distribution at the neutron binding energy, and find good agreement with the exception of $^{63}$Ni.

preprint2013arXiv

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.

preprint2013arXiv

Monte Carlo simulations on the Lefschetz thimble: taming the sign problem

We present the first practical Monte Carlo calculations of the recently proposed Lefschetz thimble formulation of quantum field theories. Our results provide strong evidence that the numerical sign problem that afflicts Monte Carlo calculations of models with complex actions can be softened significantly by changing the domain of integration to the Lefschetz thimble or approximations thereof. We study the interacting complex scalar field theory (relativistic Bose gas) in lattices of size up to 8^4 using a computationally inexpensive approximation of the Lefschetz thimble. Our results are in excellent agreement with known results. We show that - at least in the case of the relativistic Bose gas - the thimble can be systematically approached and the remaining residual phase leads to a much more tractable sign problem (if at all) than the original formulation. This is especially encouraging in view of the wide applicability - in principle - of our method to quantum field theories with a sign problem. We believe that this opens up new possibilities for accurate Monte Carlo calculations in strongly interacting systems of sizes much larger that previously possible.

preprint2013arXiv

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.

preprint2012arXiv

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.

preprint2010arXiv

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.