1
|
Hasecke L, Mata RA. Optimization of Quantum Nuclei Positions with the Adaptive Nuclear-Electronic Orbital Approach. J Phys Chem A 2024; 128:3205-3211. [PMID: 38619054 PMCID: PMC11056972 DOI: 10.1021/acs.jpca.4c00096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2024] [Revised: 03/08/2024] [Accepted: 03/08/2024] [Indexed: 04/16/2024]
Abstract
The use of multicomponent methods has become increasingly popular over the last years. Under this framework, nuclei (commonly protons) are treated quantum mechanically on the same footing as the electronic structure problem. Under the use of atomic-centered orbitals, this can lead to some complications as the ideal location of the nuclear basis centers must be optimized. In this contribution, we propose a straightforward approach to determine the position of such centers within the self-consistent cycle of a multicomponent calculation, making use of individual proton charge centroids. We test the method on model systems including the water dimer, a protonated water tetramer, and a porphine system. Comparing to numerical gradient calculations, the adaptive nuclear-electronic orbital (NEO) procedure is able to converge the basis centers to within a few cents of an Ångström and with less than 0.1 kcal/mol differences in absolute energies. This is achieved in one single calculation and with a small added computational effort of up to 80% compared to a regular NEO- self-consistent field run. An example application for the human transketolase proton wire is also provided.
Collapse
Affiliation(s)
- Lukas Hasecke
- Institute of Physical Chemistry, University
of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Ricardo A. Mata
- Institute of Physical Chemistry, University
of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| |
Collapse
|
2
|
Li TE, Paenurk E, Hammes-Schiffer S. Squeezed Protons and Infrared Plasmonic Resonance Energy Transfer. J Phys Chem Lett 2024; 15:751-757. [PMID: 38226772 DOI: 10.1021/acs.jpclett.3c03112] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2024]
Abstract
Unusual nuclear quantum effects may emerge near noble metal nanostructures such as squeezed vibrational states in molecular junctions and plasmonic resonance energy transfer in the infrared domain. Herein, nuclear quantum effects near heavy metals are studied by nuclear-electronic orbital density functional theory (NEO-DFT) with an effective core potential. For a quantum proton sandwiched between a pair of gold tips modeled by two Au6 clusters, NEO-DFT calculations suggest that the quantum proton density can be squeezed as the tip distance decreases. For an HF molecule placed near a one-dimensional Au nanowire composed of up to 34 Au atoms, real-time NEO time-dependent density functional theory (RT-NEO-TDDFT) shows that the infrared plasmonic motion within the Au nanowire may resonantly transfer electronic energy to the HF proton vibrational stretch mode. Overall, these calculations illustrate the advantages of the NEO approach for probing nuclear quantum effects, such as squeezed proton vibrational states and infrared plasmonic resonance energy transfer.
Collapse
Affiliation(s)
- Tao E Li
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Eno Paenurk
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | |
Collapse
|
3
|
Zhang Y, Wang Y, Xu X, Chen Z, Yang Y. Vibrational Spectra of Highly Anharmonic Water Clusters: Molecular Dynamics and Harmonic Analysis Revisited with Constrained Nuclear-Electronic Orbital Methods. J Chem Theory Comput 2023; 19:9358-9368. [PMID: 38096546 DOI: 10.1021/acs.jctc.3c01037] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2023]
Abstract
Vibrational spectroscopy is widely used to gain insights into structural and dynamic properties of chemical, biological, and materials systems. Thus, an efficient and accurate method to simulate vibrational spectra is desired. In this paper, we justify and employ a microcanonical molecular simulation scheme to calculate the vibrational spectra of three challenging water clusters: the neutral water dimer (H4O2), the protonated water trimer (H7O3+), and the protonated water tetramer (H9O4+). We find that with the accurate description of quantum nuclear delocalization effects through the constrained nuclear-electronic orbital framework, including vibrational mode coupling effects through molecular dynamics simulations can additionally improve the vibrational spectrum calculations. In contrast, without the quantum nuclear delocalization picture, conventional ab initio molecular dynamics may even lead to less accurate results than harmonic analysis.
Collapse
Affiliation(s)
- Yuzhe Zhang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Yiwen Wang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Xi Xu
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
| | - Zehua Chen
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
4
|
Hasecke L, Mata RA. Nuclear Quantum Effects Made Accessible: Local Density Fitting in Multicomponent Methods. J Chem Theory Comput 2023; 19:8223-8233. [PMID: 37920900 PMCID: PMC10687858 DOI: 10.1021/acs.jctc.3c01055] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2023] [Revised: 10/18/2023] [Accepted: 10/19/2023] [Indexed: 11/04/2023]
Abstract
The simulation of nuclear quantum effects (NQEs) is crucial for an accurate description of systems and processes involving light nuclei, such as hydrogen atoms. Within the last years, the importance of those effects has been highlighted for a vast range of systems with tremendous implications in chemistry, biology, physics, and materials sciences. However, while electronic structure theory methods have become routine tools for quantum chemical investigations, there is still a lack of approaches to address NQEs that are computationally accessible and straightforward to use. To address this, we present the first combination of the nuclear-electronic orbital Hartree-Fock approach with both local and density fitting approximations (LDF-NEO-HF). This results in a low-order scaling approach that enables the inclusion of NQEs for large systems within a fraction of a day and for small to medium size systems in minutes. Moreover, we demonstrate the qualitative accuracy and robustness of our approach to retrieve NQEs for three real-use cases motivated by chemical, biological, and materials science applications.
Collapse
Affiliation(s)
- Lukas Hasecke
- Institute of Physical Chemistry, University of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Ricardo A. Mata
- Institute of Physical Chemistry, University of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| |
Collapse
|
5
|
Chow M, Li TE, Hammes-Schiffer S. Nuclear-Electronic Orbital Quantum Mechanical/Molecular Mechanical Real-Time Dynamics. J Phys Chem Lett 2023; 14:9556-9562. [PMID: 37857272 PMCID: PMC11401051 DOI: 10.1021/acs.jpclett.3c02275] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2023]
Abstract
Simulating the nuclear-electronic quantum dynamics of large-scale molecular systems in the condensed phase is key for studying biologically and chemically important processes such as proton transfer and proton-coupled electron transfer reactions. Herein, the real-time nuclear-electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) approach is combined with a hybrid quantum mechanical/molecular mechanical (QM/MM) strategy to enable the accurate description of coupled nuclear-electronic quantum dynamics in the presence of heterogeneous environments such as solvent or proteins. The densities of the electrons and quantum protons are propagated in real time, while the other nuclei are propagated classically on the instantaneous electron-proton vibronic surface. This approach is applied to phenol bound to lysozyme, intramolecular proton transfer in malonaldehyde, and nonequilibrium excited-state intramolecular proton transfer in o-hydroxybenzaldehyde. These examples illustrate that the RT-NEO-TDDFT framework, coupled with an atomistic representation of the environment, allows the simulation of condensed-phase systems that exhibit significant nuclear quantum effects.
Collapse
Affiliation(s)
- Mathew Chow
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Tao E Li
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | |
Collapse
|
6
|
Xu X. Constrained Nuclear-Electronic Orbital Density Functional Theory with a Dielectric Continuum Solvent Model. J Phys Chem A 2023. [PMID: 37470267 DOI: 10.1021/acs.jpca.3c02507] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/21/2023]
Abstract
Solvent effects are crucial for simulating chemical and biological processes in solutions. The continuum solvation model is widely used for incorporating solvent effects with different levels of theoretical descriptions of solutes. For solutes and solutions containing hydrogen atoms, nuclear quantum effects can also be nonnegligible for reliable simulations. In this work, we couple our recently developed constrained nuclear-electronic orbital density functional theory with a dielectric continuum solvation model to cover nuclear quantum effects and solvent effects simultaneously. This approach is applied to the formate ion, where an anomalous solvatochromic shift in C-H stretch frequency was reported in experiments. By using this new approach to account for nuclear quantum effects and solvent effects, we show that the vibrational frequency of the C-H stretch and the solvatochromic shift are accurately described.
Collapse
Affiliation(s)
- Xi Xu
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
| |
Collapse
|
7
|
Wang Y, Chen Z, Yang Y. Calculating Vibrational Excited State Absorptions with Excited State Constrained Minimized Energy Surfaces. J Phys Chem A 2023. [PMID: 37335973 DOI: 10.1021/acs.jpca.3c01420] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2023]
Abstract
The modeling and interpretation of vibrational spectra are crucial for studying reaction dynamics using vibrational spectroscopy. Most prior theoretical developments focused on describing fundamental vibrational transitions while fewer developments focused on vibrational excited state absorptions. In this study, we present a new method that uses excited state constrained minimized energy surfaces (CMESs) to describe vibrational excited state absorptions. The excited state CMESs are obtained similarly to the previous ground state CMES development in our group but with additional wave function orthogonality constraints. Using a series of model systems, including the harmonic oscillator, Morse potential, double-well potential, quartic potential, and two-dimensional anharmonic potential, we demonstrate that this new procedure provides good estimations of the transition frequencies for vibrational excited state absorptions. These results are significantly better than those obtained from harmonic approximations using conventional potential energy surfaces, demonstrating the promise of excited state CMES-based methods for calculating vibrational excited state absorptions in real systems.
Collapse
Affiliation(s)
- Yiwen Wang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Zehua Chen
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
8
|
Solaris J, Krueger TD, Chen C, Fang C. Photogrammetry of Ultrafast Excited-State Intramolecular Proton Transfer Pathways in the Fungal Pigment Draconin Red. Molecules 2023; 28:3506. [PMID: 37110741 PMCID: PMC10144053 DOI: 10.3390/molecules28083506] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2023] [Revised: 04/11/2023] [Accepted: 04/13/2023] [Indexed: 04/29/2023] Open
Abstract
Proton transfer processes of organic molecules are key to charge transport and photoprotection in biological systems. Among them, excited-state intramolecular proton transfer (ESIPT) reactions are characterized by quick and efficient charge transfer within a molecule, resulting in ultrafast proton motions. The ESIPT-facilitated interconversion between two tautomers (PS and PA) comprising the tree fungal pigment Draconin Red in solution was investigated using a combination of targeted femtosecond transient absorption (fs-TA) and excited-state femtosecond stimulated Raman spectroscopy (ES-FSRS) measurements. Transient intensity (population and polarizability) and frequency (structural and cooling) dynamics of -COH rocking and -C=C, -C=O stretching modes following directed stimulation of each tautomer elucidate the excitation-dependent relaxation pathways, particularly the bidirectional ESIPT progression out of the Franck-Condon region to the lower-lying excited state, of the intrinsically heterogeneous chromophore in dichloromethane solvent. A characteristic overall excited-state PS-to-PA transition on the picosecond timescale leads to a unique "W"-shaped excited-state Raman intensity pattern due to dynamic resonance enhancement with the Raman pump-probe pulse pair. The ability to utilize quantum mechanics calculations in conjunction with steady-state electronic absorption and emission spectra to induce disparate excited-state populations in an inhomogeneous mixture of similar tautomers has broad implications for the modeling of potential energy surfaces and delineation of reaction mechanisms in naturally occurring chromophores. Such fundamental insights afforded by in-depth analysis of ultrafast spectroscopic datasets are also beneficial for future development of sustainable materials and optoelectronics.
Collapse
|
9
|
Li TE, Hammes-Schiffer S. Electronic Born-Oppenheimer approximation in nuclear-electronic orbital dynamics. J Chem Phys 2023; 158:114118. [PMID: 36948810 DOI: 10.1063/5.0142007] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/24/2023] Open
Abstract
Within the nuclear-electronic orbital (NEO) framework, the real-time NEO time-dependent density functional theory (RT-NEO-TDDFT) approach enables the simulation of coupled electronic-nuclear dynamics. In this approach, the electrons and quantum nuclei are propagated in time on the same footing. A relatively small time step is required to propagate the much faster electronic dynamics, thereby prohibiting the simulation of long-time nuclear quantum dynamics. Herein, the electronic Born-Oppenheimer (BO) approximation within the NEO framework is presented. In this approach, the electronic density is quenched to the ground state at each time step, and the real-time nuclear quantum dynamics is propagated on an instantaneous electronic ground state defined by both the classical nuclear geometry and the nonequilibrium quantum nuclear density. Because the electronic dynamics is no longer propagated, this approximation enables the use of an order-of-magnitude larger time step, thus greatly reducing the computational cost. Moreover, invoking the electronic BO approximation also fixes the unphysical asymmetric Rabi splitting observed in previous semiclassical RT-NEO-TDDFT simulations of vibrational polaritons even for small Rabi splitting, instead yielding a stable, symmetric Rabi splitting. For the intramolecular proton transfer in malonaldehyde, both RT-NEO-Ehrenfest dynamics and its BO counterpart can describe proton delocalization during the real-time nuclear quantum dynamics. Thus, the BO RT-NEO approach provides the foundation for a wide range of chemical and biological applications.
Collapse
Affiliation(s)
- Tao E Li
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| | | |
Collapse
|
10
|
Culpitt T, Peters LDM, Tellgren EI, Helgaker T. Time-dependent nuclear-electronic orbital Hartree-Fock theory in a strong uniform magnetic field. J Chem Phys 2023; 158:114115. [PMID: 36948801 DOI: 10.1063/5.0139675] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2023] Open
Abstract
In an ultrastrong magnetic field, with field strength B ≈ B0 = 2.35 × 105 T, molecular structure and dynamics differ strongly from that observed on the Earth. Within the Born-Oppenheimer (BO) approximation, for example, frequent (near) crossings of electronic energy surfaces are induced by the field, suggesting that nonadiabatic phenomena and processes may play a more important role in this mixed-field regime than in the weak-field regime on Earth. To understand the chemistry in the mixed regime, it therefore becomes important to explore non-BO methods. In this work, the nuclear-electronic orbital (NEO) method is employed to study protonic vibrational excitation energies in the presence of a strong magnetic field. The NEO generalized Hartree-Fock theory and time-dependent Hartree-Fock (TDHF) theory are derived and implemented, accounting for all terms that result as a consequence of the nonperturbative treatment of molecular systems in a magnetic field. The NEO results for HCN and FHF- with clamped heavy nuclei are compared against the quadratic eigenvalue problem. Each molecule has three semi-classical modes owing to the hydrogen-two precession modes that are degenerate in the absence of a field and one stretching mode. The NEO-TDHF model is found to perform well; in particular, it automatically captures the screening effects of the electrons on the nuclei, which are quantified through the difference in energy of the precession modes.
Collapse
Affiliation(s)
- Tanner Culpitt
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Laurens D M Peters
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Erik I Tellgren
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Trygve Helgaker
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| |
Collapse
|
11
|
Samsonova I, Tucker GB, Alaal N, Brorsen KR. Hydrogen-Atom Electronic Basis Sets for Multicomponent Quantum Chemistry. ACS OMEGA 2023; 8:5033-5041. [PMID: 36777583 PMCID: PMC9910068 DOI: 10.1021/acsomega.2c07782] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/06/2022] [Accepted: 01/13/2023] [Indexed: 06/18/2023]
Abstract
Multicomponent methods are a conceptually simple way to include nuclear quantum effects into quantum chemistry calculations. In multicomponent methods, the electronic molecular orbitals are described using the linear combination of atomic orbitals approximation. This requires the selection of a one-particle electronic basis set which, in practice, is commonly a correlation-consistent basis set. In multicomponent method studies, it has been demonstrated that large electronic basis sets are required for quantum hydrogen nuclei to accurately describe electron-nuclear correlation. However, as we show in this study, much of the need for large electronic basis sets is due to the correlation-consistent electronic basis sets not being optimized to describe nuclear properties and electron-nuclear correlation. Herein, we introduce a series of correlation-consistent electronic basis sets for hydrogen atoms called cc-pVnZ-mc with additional basis functions optimized to reproduce multicomponent density functional theory protonic densities. These new electronic basis sets are shown to yield better protonic densities with fewer electronic basis functions than the standard correlation-consistent basis sets and reproduce other protonic properties such as proton affinities and protonic excitation energies, even though they were not optimized for these purposes. The cc-pVnZ-mc basis sets should enable multicomponent many-body calculations on larger systems due to the improved computational efficiency they provide for a given level of accuracy.
Collapse
|
12
|
Chen Z, Yang Y. Incorporating Nuclear Quantum Effects in Molecular Dynamics with a Constrained Minimized Energy Surface. J Phys Chem Lett 2023; 14:279-286. [PMID: 36595586 DOI: 10.1021/acs.jpclett.2c02905] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The accurate incorporation of nuclear quantum effects in large-scale molecular dynamics (MD) simulations remains a significant challenge. Recently, we combined constrained nuclear-electronic orbital (CNEO) theory with classical MD and obtained a new approach (CNEO-MD) that can accurately and efficiently incorporate nuclear quantum effects into classical simulations. In this Letter, we provide the theoretical foundation for CNEO-MD by developing an alternative formulation of the equations of motion for MD. In this new formulation, the expectation values of quantum nuclear positions evolve classically on an effective energy surface that is obtained from a constrained energy minimization procedure when solving for the quantum nuclear wave function, thus enabling the incorporation of nuclear quantum effects in classical MD simulations. For comparison with other existing approaches, we examined a series of model systems and found that this new MD approach is significantly more accurate than the conventional way of performing classical MD and generally outperforms centroid MD and ring-polymer MD in describing vibrations in these model systems.
Collapse
Affiliation(s)
- Zehua Chen
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin53706, United States
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin53706, United States
| |
Collapse
|
13
|
Yu Q, Roy S, Hammes-Schiffer S. Nonadiabatic Dynamics of Hydrogen Tunneling with Nuclear-Electronic Orbital Multistate Density Functional Theory. J Chem Theory Comput 2022; 18:7132-7141. [PMID: 36378867 DOI: 10.1021/acs.jctc.2c00938] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Proton transfer reactions play a critical role in many chemical and biological processes. The development of computationally efficient approaches to describe the quantum dynamics of proton transfer, which often involves hydrogen tunneling, is challenging. Herein, the nuclear-electronic orbital multistate density functional theory (NEO-MSDFT) method is combined with both Ehrenfest and surface hopping nonadiabatic dynamics methods to describe hydrogen tunneling. The NEO-MSDFT method treats the transferring hydrogen nucleus quantum mechanically on the same level as the electrons and incorporates both static and dynamical correlation by mixing localized NEO-DFT solutions with a nonorthogonal configuration interaction scheme. The other nuclei are propagated on the NEO-MSDFT vibronic surfaces during the Ehrenfest or surface hopping dynamics. These methods are applied to proton transfer in malonaldehyde as a prototypical hydrogen tunneling system. The inclusion of vibronically nonadiabatic effects is found to significantly impact the proton transfer time and tunneling dynamics. This approach is applicable to a wide range of other proton transfer reactions.
Collapse
Affiliation(s)
- Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Saswata Roy
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | |
Collapse
|