While a detailed study of DNA solvation, including energetic and dynamic effects, is outside the scope of the paper, we demonstrate that the machine learning model can capture the interaction of a representative static snapshot for a relatively large (18 base pairs) solvated DNA system. To see this effect with our model, we use a similar strategy as in the density difference plots above to visualize the effect of solvent by taking the difference in the electrostatic potentials for the solvent-polarized (Fig. ) and gas-phase DNA structures (Fig. ). The solvent-polarized system contains over 4000 atoms, which is well-outside the range of a conventional quantum calculation, and the gas-phase system contains over 1200 atoms. For the gas-phase system, the electrostatic potential is calculated at a density isovalue of 0.001 a.u, corresponding to an average distance on the potential surface of about 1.8 Å from the nearest DNA atom. After taking the difference in the electrostatic potentials of the solvent-polarized and gas-phases DNA systems (V DN A,polarized -V DN A,gas ), the average of the resulting electrostatic potential is shifted to zero. The resulting difference is plotted in Fig. .