Codes to reproduce data for:

by Muhammad Rizwanur Rahman, Li Shen, James P. Ewen, David M. Heyes, Daniele Dini, and, Edward R. Smith (under review). Initial preprint: arXiv preprint arXiv:2311.00419

by Muhammad Rizwanur Rahman, Li Shen, James P. Ewen, Benjamin Collard, David M. Heyes, Daniele Dini, and, Edward R. Smith. Journal of Chemical Physics (Communication), April 2023. https://aip.scitation.org/doi/10.1063/5.0149974

** This repository will be further developed.

Flowmol MD solver was used for the simulations. This is an open source platform developed by one of the authors (E. R. Smith) and can be found here: https://github.com/edwardsmith999/flowmol#flowmol Once Flowmol is installed:

1. Run MD_01_equilibration.in with NVT ensemble.

2. Depending on the poking technique:

   (i) Poked-hole: run MD_02_poke_a_hole.in for few time steps (i.e, 1000) with the final_state generated from step 1.

   (ii) Cut-hole: use MD_02_cut_a_hole.py to cut a hole of target radius and target density. This python code takes the final_state file from the equilibration state and removes molecules according to the user preference:

   Remove molecules (currently this is based on a similar setup to the external force) where inputs are position [x,y,z] of the hole center, radius and then direction (x=0, y=1, z=2), ends[bottom,top] of the film, and target density i.e.: fs.remove_molecules([0.,0.,0.],8,0,[-6.5,6.5], 0.01). The python code writes a new initial_state file to be used in the production phase.

   NB. MD_02_cut_a_hole.py requires the final_state (from the equilibration phase) file in the same directory.

3. Use the final_state from step 2 as a restart file to run MD_03_production.in with NVE ensemble for sufficiently long time.

4. For spontaneous hole formation, skip step 2, and start production phase without externally poking or cutting a hole in the film. To attain rupture within computationally feasible time, start with thinner films.

5. To save finalstates frequently for the investigation of teh memory window, set the RESCUE_SNAPSHOT_FREQ to a lower value.

1. Func_hole_growth.py and Func_surface_functions.py contain the required subroutines called from other .py files.
2. hole-growth.py tracks the hole radius over time.
3. surface-tension.py measures the local surface tension of the film.
4. surface-tension-mapped-pressure-tensor.py measures local surface tension with pressuree-tensor rotated along the rim profile. Also tracks the evolution of the film profile with time.
5. radial-averaged-density.py creates images of the radial averaged density field setting the origin at the hole center (0,0,0).
6. Temperature-Velocity-radial-averaged.py plots the radial averaged velocity and temperature field.