  Disordered molecular systems such as amorphous catalysts, organic thin films,
electrolyte solutions, and water are at the cutting edge of computational
exploration today. Traditional simulations of such systems at length-scales
relevant to experiments in practice require a compromise between model accuracy
and quality of sampling. To remedy the situation, we have developed an approach
based on generative machine learning called the Morphological Autoregressive
Protocol (MAP) which provides computational access to mesoscale disordered
molecular configurations at linear cost at generation for materials in which
structural correlations decay sufficiently rapidly. The algorithm is
implemented using an augmented PixelCNN deep learning architecture that we
previously demonstrated produces excellent results in 2 dimensions (2D) for
mono-elemental molecular systems. Here, we extend our implementation to
multielemental 3D and demonstrate performance using water as our test system in
two scenarios: 1. liquid water, and 2. a sample conditioned on the presence of
a rare motif. We trained the model on small-scale samples of liquid water
produced using path-integral molecular dynamics simulation including nuclear
quantum effects under ambient conditions. MAP-generated water configurations
are shown to accurately reproduce the properties of the training set and to
produce stable trajectories when used as initial conditions in classical and
quantum dynamical simulations. We expect our approach to perform equally well
on other disordered molecular systems while offering unique advantages in
situations when the disorder is quenched rather than equilibrated.
