# Spin-Boson Hamiltonian and Optical Absorption of Molecular Dimers

## I Introduction

The spectral properties of the spin-boson Hamiltonian have frequently been adressed in the past and in the last years investigations on the spectral features of classical nonintegrability in the quantum case were performed (see e.g.  and references therein). The interest in the spectral and dynamical properties of the spin-boson Hamiltonian is due to its numerous applications which include a large variety of phenomena in molecular and solid state physics with one realization being the optical and transfer properties of excitonic-vibronic coupled dimers. Thus the spectral and related optical properties of the corresponding dimer Hamiltonian have intensively been investigated over the years, as examples we point to the work exposed in . Most of these investigations were based on a Hamiltonian corresponding to a symmetric dimer configuration. However, from the point of view of a general asymmetric situation the symmetric case is singular, e.g. the quantum states possess parity properties in the symmetric case, which are not present in general. It should be noted that dimers often constitute subunits of more complicated molecular aggregates with an asymmetric structure, which breaks the symmetry of the dimer configuration (see e.g. ). Therefore the aim of this paper is to go beyond the symmetric case and to connect the properties of the eigenstates of a generalized and symmetry-broken spin-boson Hamiltonian with the line spectrum of absorption bands of asymmetric excitonic-vibronic coupled molecular dimers.

We combine a phase space analysis of the eigenstates, which is based on the method of Husimi projections , with the calculation of the matrix elements of optical transitions. Husimi projections have repeatedly been used in order to establish the features of nonintegrability of classical systems in the corresponding quantum case . Analyzing the Husimi projections of the eigenstates of the spin-boson Hamiltonian in relation to optical transitions we indicate a bridge between the eigenstate analysis of nonintegrable systems and observables in an optical experiment. In particular, we show how the intensity variations in the Husimi projections are related to the spectral randomness in the fine structure of absorption bands. Such spectral randomness in the eigenstates of the spin-boson Hamiltonian is also known as incipience of quantum chaos . In this case random features of the spectrum are just appearing while the system is still far from displaying universal spectral fluctuations described by random matrix theory and known for excited states of spectra of polyatomic molecules .

A central point in our eigenstate analysis will be to find out to what extent the adiabatic reference systems of the spin-boson Hamiltonian are present in its exact eigenstates and how the appearance of spectral randomness can be interpreted as a mixing of such reference systems. Following  we show that Bloch projections are a useful quantitative characteristic to describe this mixing. Computing the Bloch projections of the eigenstates, it is possible to distinguish the spectral region where the adiabatic branches of the spectrum are still intact from the region with a substantial mixing of adiabatic reference systems. Furthermore, by performing projections of the numerically obtained eigenstates onto the ground state we obtain the details of the fine structure of the absorption bands in the spectral regions as indicated by the Franck-Condon principle. We note that the asymmetric case has some special features, which are not present in the symmetric situation: In the asymmetric case the final states of an optical transition can be located in the higher energetic regions of the corresponding adiabatic potentials. Then, as will be shown below for two representative cases, the fine structure of the spectrum and of the absorption bands can greatly vary from regular to irregular arrangements of absorption lines.

The paper is organized as follows: In Sec.II the model and the basic equations are presented. In Sec.III the eigenstate analysis for two representative cases of the asymmetric spin-boson Hamiltonian is performed and in Sec.IV this analysis is connected with the properties of the absorption bands in the corresponding spectral regions. Finally, in Sec.V our conclusions are summarized.

## II Model

We consider a dimer system described by the Hamiltonian

$$\widehat{H}=\underset{n}{}(ϵ_n+\gamma _nq_n)|nn|+\underset{\stackrel{n,m}{nm}}{}V_{nm}|nm|+\underset{n}{}\frac{1}{2}(p_n^2+\omega _n^2q_n^2),$$
(1)   
where $`|n`$ ($`|m`$) are the excited states of the molecular monomers constituting the dimer ($`n,m=1,2`$) with the excitonic energies $`ϵ_n`$, the coupling constants to intramolecular vibrations $`\gamma _n`$ and the transfer matrix element $`V_{nm}`$. The variables $`q_n`$, $`p_n`$ and $`\omega _n`$ represent coordinate, momentum and frequency of the vibrations, respectively. The coupling of the excited dimer states to an incident light wave with electrical field vector $`\stackrel{}{E}_n(t)`$ is given by

$$\widehat{H}_{int}=\underset{n}{}\stackrel{}{\mu }_n\stackrel{}{E}_n(t)|n0|,$$
(2)   
where $`\stackrel{}{\mu }_n`$ and $`\stackrel{}{E}_n(t)`$ are the optical transition matrix element and electric field at the nth monomer, respectively, and $`|0`$ is the excitonic vacuum. We consider an asymmetric configuration of the dimer system by assuming different monomer energies $`ϵ_n`$ and optical transition matrix elements $`\stackrel{}{\mu }_n`$. In order to keep the number of asymmetry parameters minimal we assume symmetric transfer rates, $`V_{nm}=V_{mn}=V`$, $`V>0`$, and equal coupling constants and frequencies, i.e. $`\gamma _1=\gamma _2=\gamma `$ and $`\omega _1=\omega _2=\omega `$. Then the coupling can be reduced to a single vibrational mode by introducing the relative displacement of the vibrations.

The Hamiltonian (1) can be represented as an operator in the space of two-dimensional vectors $`C=\left(\genfrac{}{}{0pt}{}{c_1}{c_2}\right)`$ using the standard Pauli spin matrices $`\widehat{\sigma }_i`$ ($`i=x,y,z`$). Passing to dimensionsless variables by measuring the energy in units of $`2V`$, $`H=H/2V`$, one obtains from (1) the spin-boson Hamiltonian

$$\widehat{H}=ϵ_+𝟙\frac{\mathrm{𝟙}}{\mathrm{𝟚}}\widehat{\sigma }_𝕩+\frac{\mathrm{𝟙}}{\mathrm{𝟚}}(\widehat{}^\mathrm{𝟚}+𝕣^\mathrm{𝟚}\widehat{}^\mathrm{𝟚})𝟙+(\sqrt{\frac{𝕡}{\mathrm{𝟚}}}𝕣\widehat{}+ϵ_{})\widehat{\sigma }_𝕫.$$
(3)   
In (3), the dimensionless relative displacement and corresponding momentum, to which the exciton is coupled, are given by $`\widehat{Q}=\sqrt{2V}(\widehat{q}_1\widehat{q}_2)`$ and $`\widehat{P}=(\widehat{p}_1\widehat{p}_2)/\sqrt{2V}`$. The dimensionless parameters of the spin-boson Hamiltonian (3) are related to the dimer system (1) by

$`p`$ $`=`$ $`{\displaystyle \frac{\gamma ^2}{2V\omega ^2}},`$ (4)   
$`r`$ $`=`$ $`{\displaystyle \frac{\omega }{2V}},`$ (5)   
$`ϵ_\pm `$ $`=`$ $`{\displaystyle \frac{ϵ_1\pm ϵ_2}{4V}}.`$ (6)   
The parameter $`p`$ expresses the coupling strength, and $`r`$ is the adiabatic parameter measuring the relative strength of quantum effects of the subsystems. The sum $`ϵ_+`$ represents the center of the excitation energy, whereas the difference $`ϵ_{}`$ is the asymmetry parameter. For $`ϵ_{}=0`$, one obtains the symmetric spin-boson Hamiltonian with conservation of total parity (given by the operation $`\widehat{Q}\widehat{Q}`$, $`\widehat{\sigma }_z\widehat{\sigma }_z`$). For $`ϵ_{}0`$, this symmetry is broken.

The adiabatic reference systems associated with (3) are introduced by considering the eigenvalue problem of the adiabatic part $`\widehat{h}_{\mathrm{ad}}`$,

$$\widehat{h}_{\mathrm{ad}}=\frac{1}{2}\widehat{\sigma }_x+(\sqrt{\frac{p}{2}}rQ+ϵ_{})\widehat{\sigma }_z+\frac{r^2}{2}Q^2𝟙.$$
(7)   
The $`Q`$-dependent eigenvalues and eigenstates of $`\widehat{h}_{\mathrm{ad}}`$ are easily obtained. In particular, one finds the two adiabatic potentials $`U_{\mathrm{ad}}^\pm (Q)`$ for the upper (+) and lower (-) states

$$U_{\mathrm{ad}}^\pm (Q)=\frac{r^2}{2}Q^2\pm \sqrt{\frac{1}{4}+\left(ϵ_{}+\sqrt{\frac{p}{2}}rQ\right)^2},$$
(8)   
from which the Hamiltonians of the adiabatic reference systems are obtained

$$H^\pm (Q)=\frac{1}{2}P^2+U_{\mathrm{ad}}^\pm (Q).$$
(9)   
A key point in our analysis of the eigenstates and absorption bands will be the extent to which the exact eigenstates of the Hamiltonian (3) can be viewed as a mixing of adiabatic reference systems connected with the adiabatic Hamiltonians (9). As is shown below, the extent of this mixing can be controlled by analyzing Husimi projections and by computing Bloch projections from the numerically obtained eigenstates. Husimi projections constitute phase space representations of the eigenstates which can be compared with the phase space orbits of the adiabatic reference systems (9). Bloch projections constitute another independent indicator of the mixing of the adiabatic branches in the spectrum. For the adiabatic branches these projections are calculated from the eigenstates of the Hamiltonian (7)

$`|\phi _{\mathrm{ad}}^+(Q)`$ $`=`$ $`{\displaystyle \frac{1}{\sqrt{2}}}\left(\sqrt{1+A(Q)}|\sqrt{1A(Q)}|\right)`$ (11)   
$`|\phi _{\mathrm{ad}}^{}(Q)`$ $`=`$ $`{\displaystyle \frac{1}{\sqrt{2}}}\left(\sqrt{1A(Q)}|+\sqrt{1+A(Q)}|\right)`$ (12)   
where

$$A(Q)=\frac{2ϵ_{}+\sqrt{2p}rQ}{\sqrt{1+(2ϵ_{}+\sqrt{2p}rQ)^2}}.$$
(13)   
Computing the expectation values $`x=\widehat{\sigma }_x`$ of the Pauli matrix $`\widehat{\sigma }_x`$ for the adiabatic eigenstates (II), one obtains the Bloch projections of the adiabatic branches

$$\phi _{\mathrm{ad}}^\pm (Q)|\widehat{\sigma }_x|\phi _{\mathrm{ad}}^\pm (Q)=\frac{1}{\sqrt{1+(2ϵ_{}+\sqrt{2p}rQ)^2}},$$
(14)   
According to eq.(14) there is a distinction in the sign of the Bloch projections for the adiabatic eigenstates. Computing the Bloch projections from the numerically obtained eigenstates this distinction in the sign will be used as an indicator for the case in which the spectrum resolves into adiabatic branches.

## III Eigenstate analysis

Our analysis is based on a diagonalization of the spin-boson Hamiltonian using a basis of product states of spin states with harmonic oscillator eigenstates. In this basis the eigenstates are obtained in the form

$$|\lambda =\underset{m}{}\left(c_m^{(\lambda )}|+c_m^{(\lambda )}|\right)|m=\underset{m}{}|s_m^{(\lambda )}|m,$$
(15)   
where $`|=\left(\genfrac{}{}{0pt}{}{1}{0}\right)`$, $`\left(|=\left(\genfrac{}{}{0pt}{}{0}{1}\right)\right)`$ denote the spin up (down) states, $`|m`$ the harmonic oscillator eigenstates and $`c_{zm}^{(\lambda )}`$ ($`z=,`$) the expansion coefficients of the $`\lambda `$ eigenstate. For a fixed $`m`$ we will also use the spin representation $`|s_m^{(\lambda )}`$ for the up and down coefficients $`c_{zm}^{(\lambda )}`$ ($`z=,`$) for a given oscillator state as indicated in the second part of (15). The matrix dimension in the numerical diagonalization was $`N=4000`$, for the eigenstate analysis the first $`1100`$ eigenstates were used.

For the phase space representation of the eigenstates Husimi projections are used. Husimi projections $`h_\lambda (\alpha (Q,P);s)`$ are defined by projecting an eigenstate $`|\lambda `$ onto a set of coherent oscillator states $`|\alpha (Q,P)`$, which scan the oscillator phase plane of the vibrational subsystem, while the spin projection $`|s=c_{}|+c_{}|`$ is kept fixed,

$$h_\lambda (\alpha (Q,P);s)=|\lambda |\alpha (Q,P),s|^2,$$
(16)   
and $`|\alpha (Q,P),s=|\alpha (Q,P)|s`$. Using the explicit form of $`|\lambda `$ in (15) one finds

$$h_\lambda (\alpha (Q,P);s)=\left|\underset{m}{}(c_m^{(\lambda )}c_{}+c_m^{(\lambda )}c_{})m|\alpha (Q,P)\right|^2,$$
(17)   
where $`m|\alpha (Q,P)=\frac{\alpha ^m}{\sqrt{m!}}\mathrm{e}^{\frac{|\alpha |^2}{2}}`$ and $`\alpha (Q,P)=\sqrt{\frac{r}{2}}\widehat{Q}+\frac{i}{\sqrt{2r}}\widehat{P}`$.

For the Bloch projection of an eigenstate,

$$x_\lambda =\lambda \left|\widehat{\sigma }_x\right|\lambda ,$$
(18)   
one finds from the expansion (15)

$$x_\lambda =\underset{m}{}\left(c_m^{(\lambda )}c_m^{(\lambda )}+c_m^{(\lambda )}c_m^{(\lambda )}\right).$$
(19)   
Before turning to the eigenstate analysis for definite parameter sets we note that the appearance of spectral irregularities can be expected in the high energetic region of the overlap of both adiabatic potentials only, where a sufficient number of states of both potentials mix. It is obvious that such a region can be probed by the final states of an optical transition in the asymmetric case only: In the symmetric case, as is evident from the adiabatic potentials (8), the optical transitions would necessarily terminate in the ground state region of the upper potential, where only a few states of the upper potential can mix with the lower potential. Therefore when investigating the structure of the spectrum and eigenstates of the spin-boson Hamiltonian we paid special attention to a systematic change of the asymmetry parameter $`ϵ_{}`$ combined with the coupling parameter $`p`$ in order to produce appropriate configurations of the adiabatic potentials. Such configurations were found for relatively high coupling $`p`$ and asymmetry values $`ϵ_{}`$. Below we present the eigenstate analysis for two typical asymmetric configurations in the adiabatic parameter region ($`r<1`$) with the parameter sets

A: $`p=4`$, $`r=0.1`$ and $`ϵ_{}=5`$ and

B: $`p=20`$, $`r=0.1`$ and $`ϵ_{}=10`$.

The parameter set B corresponds to a stronger coupling value $`p`$ and a larger asymmetry parameter $`ϵ_{}`$, as compared to set A. To indicate the location of the absorption bands, which will be analyzed in section IV, the adiabatic potentials of set B and the ground state of the system are shown in Fig. 1.

In Fig. 2 we start with the eigenstate analysis for the parameter set A by displaying the Bloch projections $`x(E_\lambda )`$ computed from the eigenstates according to eq.(19). In the low energy region one finds one adiabatic branch in the $`x(E_\lambda )`$ dependence associated with the lower adiabatic potential (for which $`x_\lambda >0`$ in accordance with the sign in (14)). The overlap of the two potentials is marked by the appearance of a second adiabatic branch in the $`x(E_\lambda )`$ dependence. The mixing of the eigenstates of the two adiabatic potentials is seen by the disappearance of both adiabatic branches and indicated by a broad band of Bloch projections in the high energy region. The presence of three characteristic regions in the Bloch projections with one adiabatic branch, two adiabatic branches and the mixing region is reflected in the Husimi projections of the eigenstates. In the regions where the adiabatic branches are intact one obtains regular Husimi projections concentrated around the phase space orbits of the integrable adiabatic reference Hamiltonians (9). Typical examples of such regular projections corresponding to equal spin projections $`c_{}=c_{}`$ for the energetic overlap region with two adiabatic branches present are shown in Fig. 3(a) and Fig. 3(b). In Fig. 3(a) the Husimi projection corresponding to an eigenstate in the upper adiabatic potential is shown, whereas in Fig. 3(b) a projection associated with an eigenstates of the lower adiabatic potential is displayed. Sequences of such projections clearly attributable to the coexistence of two independent adiabatic branches were observed until the mixing region is reached. Then Husimi projections incorporating the phase space features of both adiabatic potentials are observed, as shown in Fig. 3(c).

A closer inspection of these projections in the mixing region shows that although the distribution is concentrated on the regular phase space orbits of the adiabatic potentials, the intensity between these parts varies in an irregular and random way, when Husimi projections for a sequence of eigenstates are considered. The structure of the eigenstates which shows up in the random shift of intensity in the Husimi projections is a source of spectral randomness in the absorption bands. In the case of the parameter set A, however, this region of an irregular and random shift of the phase space distribution is not reached by the absorption processes. In this case the absorption band is located in the region, where the two adiabatic branches are still intact and no spectral mixing between their eigenstates occurs.

In the case of the parameter set B one finds a situation, where the spectral window, which is cut by the transition matrix elements of the absorption process, reaches the region of an irregular mixing of adiabatic states. The Husimi projections for selected eigenstates and equal spin projections which correspond to such a region and which are relevant as final states of an optical transition are shown in the set of Fig. 4(a)-(c). One observes a random variation of the intensity in the Husimi distributions of the eigenstates between the orbits of the two adiabatic potentials. As we show in the next section, as a consequence of this random variation intensities of lines for optical transitions terminating in these states vary randomly.

## IV Absorption

In the calculation of the optical absorption the dimer system is assumed to be small compared to the optical wavelength. Then the electric field vector is approximately equal at the monomer sites of the dimer configuration, i.e. $`\stackrel{}{E}_n(t)=\stackrel{}{E}(t)`$. Introducing the projections $`\mu _n`$ of the monomer transition matrix elements onto the direction of the field polarization, i.e. $`\stackrel{}{E}(t)\stackrel{}{\mu }_n=E(t)\mu _n`$, ($`n=1,2`$), the interaction of the dimer with the incident light wave is represented in the form

$$\widehat{H}_{int}=E(t)\widehat{M}_{int},$$
(20)   
where $`\widehat{M}_{int}`$ is the interaction operator

$$\widehat{M}_{int}=\underset{n}{}|n\mu _n0|.$$
(21)   
In the ground state there is no excitonic-vibronic interaction, i.e. the ground state of the system $`|g`$ is given by the direct product of the excitonic vacuum $`|0=|exc0`$ and the vibrational ground state. For the ground state of the vibrational subsystem we assume the zero temperature case, where it is in its lowest state $`|m=0`$, with $`|m=0`$ the $`m=0`$ Hermitian polynomial of the undisplaced $`Q`$ oscillator. Then $`|g=|exc0|m=0`$. For the optical transition matrix element $`M_{\lambda g}=\lambda |\widehat{M}_{int}|g`$ of the interaction operator (21), which corresponds to an absorption process from the ground state $`g`$ into an excited state $`\lambda `$, one then obtains using eq.(15) the expression

$$M_{\lambda g}=\mu _1c_0^{(\lambda )}+\mu _2c_0^{(\lambda )}.$$
(22)   
The matrix element (22) is representable in the form of a particular projection of the spin vector associated with the $`m=0`$ vibrational state in (15), $`|s_0^{(\lambda )}`$. Introducing the angles

$$\mathrm{cos}\alpha =\frac{\mu _1}{\sqrt{\mu _1^2+\mu _2^2}},\mathrm{sin}\alpha =\frac{\mu _2}{\sqrt{\mu _1^2+\mu _2^2}}$$
(23)   
and a spin vector $`|s_\mu `$ defined by

$$|s_\mu =\left(\genfrac{}{}{0pt}{}{\mathrm{cos}\alpha }{\mathrm{sin}\alpha }\right),$$
(24)   
one represents (22) in the form

$$M_{\lambda g}=\sqrt{\mu _1^2+\mu _2^2}s_\mu |s_0^{(\lambda )}.$$
(25)   
Measuring the square of (25), i.e. the absorption strength of an optical transition from $`|g`$ to $`|\lambda `$, in units of $`(\mu _1^2+\mu _2^2)`$ one finds for the dimensionsless absorption strength $`Q_{\lambda g}`$

$$Q_{\lambda g}=\frac{M_{\lambda g}^2}{\mu _1^2+\mu _2^2}=(s_\mu |s_0^{(\lambda )})^2.$$
(26)   
Comparing eq.(26) with the expression for the Husimi projection (17) it follows that $`Q_{\lambda g}`$ is equal to a Husimi projection of the $`\lambda `$-eigenstate of (3) taken at the phase space point $`Q=P=0`$ with the spin projection being fixed at $`s_\mu `$ : For $`Q=P=0`$ one obtains $`\alpha (Q,P)=0`$, which selects the $`m=0`$ term in the sum of (17), inserting the value $`s_\mu `$ for the spin vector $`s`$, one finds that $`h_\lambda (0;s)`$ reduces to the r.h.s. of eq.(26), i.e.

$$Q_{\lambda g}=h_\lambda (0;s_\mu ).$$
(27)   
Eq. (27) is in line with the Franck-Condon principle for optical transitions: The transition is vertical from the ground state region at $`Q=P=0`$ in the phase space representation and probes the intensity of the final state by its Husimi distribution in the same region. The optical transition matrix elements of the monomers define the spin projection $`s_\mu `$. For the special case of a symmetric dimer with $`\mu _1=\mu _2`$ the transition occurs into the symmetric combination of monomer states. Differences in the optical transition matrix elements, i.e. $`\mu _1\mu _2`$, introduce a second asymmetry parameter besides that of the site energy asymmetry $`ϵ_{}`$. These differences in the optical transition matrix elements express the optical asymmetry of the dimer and enter the optical dimer matrix element (22) through the spin projections (23).

In Fig. 5(a) the absorption bands $`Q_{\lambda g}`$ calculated from the numerically obtained eigenstates for the parameter set A are represented as stick spectra. Absorption bands are located in accordance with the Franck-Condon principle. This is seen by comparing their positions with the final states following from the adiabatic potentials in (8) by setting $`Q=0`$. One finds an energetically lower band for optical transitions terminating in the lower adiabatic potential and a higher band for transitions occuring into the region of the overlap of the two potentials. The shape of the lower band on the left side of Fig. 5(a) is completely regular. The lines of the higher band on the right side of Fig. 5(a) are also regularly arranged. A closer inspection of the higher band, however, reveals that between the lines shown in Fig. 5(a) lines with a much weaker intensity are embedded, which form another regular band. This is evident from the inset, Fig. 5(b), where parts of the lines of the weak and strong bands are shown together. In order to make the weak band visible in the inset a much smaller intensity scale is used. The two bands, into which the absorption spectrum of the high energy part resolves, are easily identified as to belong to the states of the two coexisting adiabatic branches analyzed in the preceding section for the parameter set A: The weak band with relatively small intensities is due to optical transitions terminating in the high energy states of the lower potential (the sign of the Bloch variables of the final states is $`x_\lambda >0`$ in accordance with the sign of the lower adiabatic branch (14)), whereas the strong band with much greater intensities is due to optical transitions into the low energy states of the upper potential (for these states $`x_\lambda <0`$, again in accordance with the sign of the adiabatic branch (14)).

In Fig. 6 the absorption bands are shown for the parameter set B. The bands are located in the spectral regions of the Franck-Condon energies indicated by the arrows in Fig. 1. The lower absorption band, shown on the left side of the Fig. 6, is regular like the lower band in the Fig. 5(a). The higher band, displayed on the right side of Fig. 6, however, is completely irregular and cannot be resolved into independent subbands as in the case of the parameter set A. This is a consequence of the mixing of the adiabatic reference systems and the irregular structure of the eigenstates analyzed by the Husimi projections in the preceding section. In particular, these projections (see Fig. 4) show a random variation of the intensity in the $`Q=P=0`$ region, which according to eq.(27) is relevant for the final states of the optical transitions. This random variation is probed by optical transition matrix elements and results in an irregular pattern of lines in the fine structure of the high energy absorption band.

The absorption bands in Fig. 5 and Fig. 6 were calculated for equal optical transition matrix elements at the monomers $`\mu _1=\mu _2`$, which correspond to spin projections of the eigenstates with equal components. The changes occuring in the lower and upper absorption bands in the limiting cases of an optical asymmetry with $`\mu _10,\mu _2=0`$ and $`\mu _1=0,\mu _20`$ for the case of the parameter set B are compared in the Figs. 7 (a) and (b), respectively. These limiting cases correspond to an optical asymmetry, when one of the monomers constituting the dimer is optically active only. Then according to eq.(22) the optical transition matrix element is determined by either the spin up $`c_0^{(\lambda )}`$ or the spin down $`c_0^{(\lambda )}`$ coefficients of the $`\lambda `$ eigenstate. As is seen from Fig. 7, which is representative for the case of a positive sign of the excitation energy asymmetry $`ϵ_{}`$ ($`ϵ_{}=10>0`$ in the case shown in Fig. 7), the lower and upper bands show a redistribution of their intensities with the optical asymmetry: For the given sign of the excitation energy asymmetry $`ϵ_{}>0`$ ($`ϵ_{}<0`$) the intensity of the upper band is increased (decreased) compared to the symmetric case, $`\mu _1=\mu _2`$, displayed in Fig. 6. The lower band behaves in the opposite way, the peak intensities of both bands differing by a factor of about $`10^3`$. In particular, maximum absorption is reached for the high energy band, when the monomer with the higher excitation energy is optically active only, i.e. $`ϵ_{}>0`$ with $`\mu _10`$ and $`\mu _2=0`$. Reversing the optical activity of the monomers but still considering the same asymmetry in the excitation energies, i.e. $`ϵ_{}>0`$ with $`\mu _1=0`$ and $`\mu _20`$, one finds maximum absorption for the low energy band (if $`ϵ_{}<0`$, the quantities $`\mu _1`$ and $`\mu _2`$ have to be interchanged in the above consideration). As in the case of equal transition matrix elements the lower bands shown on the left sides of Fig. 7(a) and Fig. 7(b) display a regular fine structure, whereas the fine structure of the upper bands, displayed on the right sides of Fig. 7(a) and Fig. 7(b) is irregular. We note that the appearance of the weak bands is due to the coupling of the excited states of both monomers: The monomer, which is optically not active, is still present in the absorption spectrum due to the coupling of the excited states in the dimer system. In the case displayed on the right side of Fig. 7(b), e.g. an optically nonactive monomer with a high excitation energy attached to an optically active monomer with a lower excitation energy produces a weak but irregular band in the dimer system. A closer inspection of the fine structure of the irregular bands for the limiting cases of an optical asymmetry, in Fig. 7(a) and Fig. 7(b), together with the irregular band for the symmetric case, shown in Fig. 6, reveals that the sequence of lines with a strong and small absorption strength is almost identical in all bands, there is only a small overall change in the line intensities in all three representations. The reason behind this behavior becomes evident after an inspection of the ratio of the spin projections

$$r^{(\lambda )}=c_0^{(\lambda )}/c_0^{(\lambda )},$$
(28)   
which is presented in Fig. 8 as a function of the eigenstate energy $`E_\lambda `$ for the spectral regions corresponding to the lower and upper absorption bands, respectively. As is evident from Fig. 8 $`r^{(\lambda )}`$ is a smooth function of the energy $`E_\lambda `$ in the spectral regions of both bands. In particular, for the case of the irregular upper absorption bands, for which in a sequence of eigenstates the spin projection coefficients $`c_0^{(\lambda )}`$ and $`c_0^{(\lambda )}`$ vary in an irregular way, their ratio $`r^{(\lambda )}`$ considered between neighbouring eigenstates is practically identical. This correlation between the spin up and spin down ($`c_0^{(\lambda )}`$ and $`c_0^{(\lambda )}`$) coefficients results in an almost smooth interpolation of the fine structure of the absorption bands between the limiting cases displayed in the Fig. 7. The smoothness of the function $`r^{(\lambda )}`$ can be used to express the absorption strength for arbitrary optical asymmetry by one of the limiting cases. For example, representing the down coefficients $`c_0^{(\lambda )}`$ by the up coefficients $`c_0^{(\lambda )}`$ in (22) using eq. (28), one obtains

$$Q_{\lambda g}=\frac{(\mu _1+\mu _2r^{(\lambda )})^2}{\mu _1^2+\mu _2^2}Q_{\lambda g}^+,$$
(29)   
where $`Q_{\lambda g}^+=[c_0^{(\lambda )}]^2`$ is the absorption strength for the case $`\mu _10`$ and $`\mu _2=0`$. The smoothness of the ratio of the different spin projections $`r^{(\lambda )}`$ together with eq. (29) allows an interpolation between absorption bands for different optical asymmetries. In particular, using this ratio in the prefactor in front of $`Q_{\lambda g}^+`$ the relative intensities of absorption bands corresponding to different optical asymmetries can be determined.

## V Conclusions

The transition from the coexistence to a mixing of the adiabatic branches in the eigenstates of a symmetry-broken spin-boson Hamiltonian can be controlled by Bloch and Husimi projections and it shows up in the structure of absorption bands in excitonic-vibronic coupled dimer systems with asymmetric adiabatic potential configurations. In the mixing region of the adiabatic branches phase space distributions of the eigenstates are concentrated on the phase space orbits of both adiabatic potentials and their intensity varies in a random and irregular way. This random variation in the intensitiy distribution of the eigenstates is probed by the optical transition matrix elements and it is a source of spectral randomness in the fine structure of the absorption bands. The asymmetry of the optical transition matrix elements of the dimer configuration is representable by suitably chosen spin projections of the spin part of the eigenstates. In the limiting cases of the optical asymmetry, in which one of the monomers constituting the dimer is optically active only, the intensities of the absorption lines are determined by either the spin up or spin down coefficients. For the states relevant for the optical absorption correlated spin down and spin up coefficients are obtained. As a result the shape of the fine structure of the irregular absorption bands remains almost the same for different optical asymmetries. Finally we point out that in the case of a symmetry-broken spin-boson Hamiltonian optical absorption processes can probe the high energy regions of the overlap of adiabatic potentials where it is possible to observe a transition from regular to irregular absorption spectra for asymmetric molecular dimer configurations.