Single-photon absorption and emission from a natural photosynthetic complex.

Photosynthesis is generally assumed to be initiated by a single photon1-3 from the Sun, which, as a weak light source, delivers at most a few tens of photons per nanometre squared per second within a chlorophyll absorption band1. Yet much experimental and theoretical work over the past 40 years has explored the events during photosynthesis subsequent to absorption of light from intense, ultrashort laser pulses2-15. Here, we use single photons to excite under ambient conditions the light-harvesting 2 (LH2) complex of the purple bacterium Rhodobacter sphaeroides, comprising B800 and B850 rings that contain 9 and 18 bacteriochlorophyll molecules, respectively. Excitation of the B800 ring leads to electronic energy transfer to the B850 ring in approximately 0.7 ps, followed by rapid B850-to-B850 energy transfer on an approximately 100-fs timescale and light emission at 850-875 nm (refs. 16-19). Using a heralded single-photon source20,21 along with coincidence counting, we establish time correlation functions for B800 excitation and B850 fluorescence emission and demonstrate that both events involve single photons. We also find that the probability distribution of the number of heralds per detected fluorescence photon supports the view that a single photon can upon absorption drive the subsequent energy transfer and fluorescence emission and hence, by extension, the primary charge separation of photosynthesis. An analytical stochastic model and a Monte Carlo numerical model capture the data, further confirming that absorption of single photons is correlated with emission of single photons in a natural light-harvesting complex.

Molecular Mechanics Simulations and Improved Tight-Binding Hamiltonians for Artificial Light Harvesting Systems: Predicting Geometric Distributions, Disorder, and Spectroscopy of Chromophores in a Protein Environment.

We present molecular mechanics and spectroscopic calculations on prototype artificial light harvesting systems consisting of chromophores attached to a tobacco mosaic virus (TMV) protein scaffold. These systems have been synthesized and characterized spectroscopically, but information about the microscopic configurations and geometry of these TMV-templated chromophore assemblies is largely unknown. We use a Monte Carlo conformational search algorithm to determine the preferred positions and orientations of two chromophores, Coumarin 343 together with its linker and Oregon Green 488, when these are attached at two different sites (104 and 123) on the TMV protein. The resulting geometric information shows that the extent of disorder and aggregation properties and therefore the optical properties of the TMV-templated chromophore assembly are highly dependent on both the choice of chromophores and the protein site to which they are bound. We use the results of the conformational search as geometric parameters together with an improved tight-binding Hamiltonian to simulate the linear absorption spectra and compare with experimental spectral measurements. The ideal dipole approximation to the Hamiltonian is not valid because the distance between chromophores can be very small. We found that using the geometries from the conformational search is necessary to reproduce the features of the experimental spectral peaks.

A deterministic alternative to the full configuration interaction quantum Monte Carlo method.

Development of exponentially scaling methods has seen great progress in tackling larger systems than previously thought possible. One such technique, full configuration interaction quantum Monte Carlo, is a useful algorithm that allows exact diagonalization through stochastically sampling determinants. The method derives its utility from the information in the matrix elements of the Hamiltonian, along with a stochastic projected wave function, to find the important parts of Hilbert space. However, the stochastic representation of the wave function is not required to search Hilbert space efficiently, and here we describe a highly efficient deterministic method that can achieve chemical accuracy for a wide range of systems, including the difficult Cr2 molecule. We demonstrate for systems like Cr2 that such calculations can be performed in just a few cpu hours which makes it one of the most efficient and accurate methods that can attain chemical accuracy for strongly correlated systems. In addition our method also allows efficient calculation of excited state energies, which we illustrate with benchmark results for the excited states of C2.

Rovibrational spectra for the HCCCN*HCN and HCN*HCCCN binary complexes in 4He droplets.

Rovibrational spectra are measured for the HCCCN*HCN and HCN*HCCCN binary complexes in helium droplets at low temperature. Though no Q-branch is observed in the infrared spectrum of the linear HCN*HCCCN dimer, which is consistent with previous experimental results obtained for other linear molecules, a prominent Q-branch is found in the corresponding infrared spectrum of the HCCCN*HCN complex. This Q-branch, which is reminiscent of the spectrum of a parallel band of a prolate symmetric top, implies that some component of the total angular momentum is parallel to the molecular axis. The appearance of this particular spectroscopic feature is analyzed here in terms of a nonsuperfluid helium density induced by the molecular interactions. Finite temperature path integral Monte Carlo simulations are performed using potential energy surfaces calculated with second-order Möller-Plesset perturbation theory, to investigate the structural and superfluid properties of both HCCCN*HCN(4He)N and HCN*HCCCN(4He)N clusters with N < or = 200. Explicit calculation of local and global nonsuperfluid densities demonstrates that this difference in the rovibrational spectra of the HCCCN*HCN and HCN*HCCCN binary complexes in helium can be accounted for by local differences in the superfluid response to rotations about the molecular axis, i.e., different parallel nonsuperfluid densities. The parallel and perpendicular nonsuperfluid densities are found to be correlated with the locations and strengths of extrema in the dimer interaction potentials with helium, differences between which derive from the variable extent of polarization of the CN bond in cyanoacetylene and the hydrogen-bonded CH unit in the two isomers. Calculation of the corresponding helium moments of inertia and effective rotational constants of the binary complexes yields overall good agreement with the experimental values.