Theory of Photoinduced Phase Transitions in Molecular Conductors : Interplay Between Correlated Electrons , Lattice Phonons and Molecular Vibrations

Dynamics of photoinduced phase transitions in molecular conductors are reviewed from the perspective of interplay between correlated electrons and phonons. (1) The charge-transfer complex TTF-CA shows a transition from a neutral paraelectric phase to an ionic ferroelectric phase. Lattice phonons promote this photoinduced transition by preparing short-range lattice dimerization as a precursor. Molecular vibrations stabilize the neutral phase so that the ionic phase, when realized, possesses a large ionicity and the Mott character; (2) The organic salts θ-(BEDT-TTF)2RbZn(SCN)4 and α-(BEDT-TTF)2I3 show transitions from a charge-ordered insulator to a metal. Lattice phonons make this photoinduced transition hard for the former salt only. Molecular vibrations interfere with intermolecular transfers of correlated electrons at an early stage; (3) The organic salt κ-(d-BEDT-TTF)2Cu[N(CN)2]Br shows a transition from a Mott insulator to a metal. Lattice phonons modulating intradimer transfer integrals enable photoexcitation-energy-dependent transition pathways through weakening of effective interaction and through introduction of carriers.

Crystals 2012, 2 57 1. Introduction Most of organic conductors are now recognized as strongly correlated electron systems, which show a variety of electronic phases depending on temperature T , pressure P , constituent elements, etc.In order to pursue further possibilities of organic conductors, the feasibility of manipulating the electronic phases on designed spatial and temporal scales is important.Phase transitions are induced in equilibrium by the variation of parameters such as T and P .They can be photoinduced, under nonequilibrium environment, on different time scales ranging from femtoseconds to nanoseconds.
Photoinduced phase transitions are not limited to molecular conductors.However, molecular systems constitute an important class of materials, the dynamic characteristics of which are well investigated both experimentally and theoretically.So, let us summarize the characteristics of their electronic properties.First of all, electric conduction is realized by an overlap between neighboring (highest occupied or lowest unoccupied) molecular orbitals, so that the stacking structure of molecules is important.The overlap is usually confined in particular directions, so that the electric conduction is limited to low-dimensional space.As is well known in condensed matter physics, the instability of the Fermi sea becomes stronger as the spatial dimension becomes lower.Thus, a tendency to form a long-range order, making the system an insulator or a superconductor, is strengthened.Organic conductors can then possess a variety of electronic phases.In some cases, the dimensionality is continuously tuned by changing T , P , and constituent elements [1][2][3].Furthermore, in quasi-two-dimensional molecular conductors, molecular arrangements are continuously changed from squares to triangles, namely, geometric frustration is tuned to modify the phase diagram itself [4,5].
Secondly, molecules, the building blocks of molecular conductors, are large objects and the intermolecular overlaps are generally small, so that the ratios of on-site Coulomb energies to the transfer integrals are large compared with those in transition-metal oxides.This makes electrons correlated strongly with each other.The highest occupied or lowest unoccupied molecular orbitals generally have many nodes, so that the transfer integrals are sensitive to the intermolecular distance and the relative orientation.This makes electrons to couple strongly with lattice phonons.Molecules themselves are soft and the molecular orbitals are sensitive to the bond lengths inside them.This makes electrons to couple strongly with molecular vibrations.Furthermore, these interaction strengths are comparable, which causes the variety of the electronic phases realized in molecular conductors.
As for time evolutions studied theoretically, long-term stochastic evolutions of statistically averaged quantities were treated by master equations [33].Meanwhile, short-term deterministic evolutions of wave functions are treated by the time-dependent Schrödinger equation.Different approximations are introduced to treat wave functions.As the time scale is shortened, electron correlations become important and need to be treated more exactly, as well as the quantum nature of phonons, so that the dimensions of the wave function and the evolution operator are increased in an exponential manner with respect to the system size.Quite recently, with a time resolution of 10 fs, electronic motion and its interference with molecular vibrations have been observed [22].Coherent control of electronic phases may be realized in some molecular material in the near future.
In order to take full advantage of the characteristics of molecular conductors, intramolecular or intradimer degrees of freedom should be exploited.Thanks to the improved time resolution, the evolution of electronic spectra modulated by molecular vibrations is clarified, which gives us information about how the ground state is stabilized by the electron-molecular-vibration couplings [22,27,28].Intramolecular or intradimer degrees of freedom may produce different transition pathways depending on the way by which the system is photoexcited.For instance, the use of molecular degrees of freedom inside a dimer in dimer-Mott insulators enables photoinduced insulator-to-metal transitions through the weakening of effective interaction and the introduction of carriers [15,16].The existence of at least two hierarchies consisting of intra-and inter-molecular degrees of freedom will become a central issue for relaxation processes required for photoinduced phase transitions.
In this context, from the perspective of interplay between intra-and inter-molecular degrees of freedom and interplay among correlated electrons, lattice phonons and molecular vibrations, we review the following photoinduced phase transitions realized in molecular conductors.(1) The photoinduced neutral-to-ionic transition in the mixed-stack charge-transfer complex tetrathiafulvalene-p-chloranil (TTF-CA) is from a paraelectric phase to a ferroelectric phase, and basically from a band insulator to a Mott insulator.It is shown that, in order for the ionic phase to be a typical Mott insulator, electron-molecular-vibration couplings need to be large, stabilizing the neutral phase, which is counterintuitive in the sense that the neutral phase is a band insulator.Lattice phonons are responsible for dimerization and consequently for the ferroelectric ground state.Above the neutral-ionic transition temperature, dimerization is short-range, but it promotes the photoinduced growth of a dimerized ionic domain as a seed; (2) The charge orders in the organic salts θ-(BEDT-TTF) 2 RbZn(SCN) 4 and α-(BEDT-TTF) 2 I 3 [BEDT-TTF=bis(ethylenedithio)tetrathiafulvalene] both show a horizontal-stripe pattern and are quite similar.Their photoinduced dynamics are, however, quite different.It is shown that their slightly different crystal structures make the effects of electron-lattice couplings quite different, because the way by which the charge order is stabilized is different, in spite of the fact that the charge order is basically stabilized by Coulomb interactions.At an early stage, collective motion of electron transfers is observed and it interferes with molecular vibrations; (3) The Mott insulator phase in κ-(d-BEDT-TTF) 2 Cu[N(CN) 2 ]Br can be converted into a metallic phase by photoexcitation.In general, this transition is induced by the weakening of effective interaction relative to the bandwidth or the introduction of carriers away from half filling.Both transition pathways are realized by tuning the energy of photoexcitation, using intradimer and interdimer charge-transfer excitations in the dimer-Mott insulator phase.

Neutral-Ionic Transition in TTF-CA
For TTF-CA, dimerization is essential to the ferroelectric ground state in the ionic phase.To explain the dimerization, different types of electron-lattice couplings have been employed [34][35][36][37].Lattice modulation of transfer integrals triggers the spin-Peierls mechanism for dimerization [34,36,37], while lattice modulation of Coulomb interaction strengths can also cause dimerization [35].Photoinduced neutral-to-ionic transition dynamics has been calculated on the basis of the latter mechanism [38].However, it cannot reproduce the quick growth and oscillation of the lattice dimerization experimentally observed after photoexcitation [27].Therefore, we use a model based on the former mechanism.It is noted that the long-range Coulomb interaction is important for the experimentally observed [39], coherent motion of a macroscopic neutral-ionic domain boundary [26].
The one-dimensional half-filled extended ionic Peierls-Holstein-Hubbard model used here is written as where c † jσ creates an electron with spin σ at site j, n jσ = c † jσ c jσ , and n j = σ n jσ .The parameter t 0 denotes the transfer integral on a regular lattice, ∆ the site energy difference between neighboring orbitals when molecular distortions are absent, U the on-site repulsion strength, and V the nearest-neighbor repulsion strength.The lattice displacement u j at site j modulates the transfer integral between the (j − 1)th and jth orbitals and that between the jth and (j + 1)th orbitals with the coefficient ∓α.The displacement v For the model parameters, we take eV as the unit of energy and use t 0 = 0.17, U = 1.5, and V = 0.6; we vary ∆ around the boundary between the neutral and ionic phases shown in Figure 1.

Neutral paraelectric
We define the strengths of these couplings as βj .The displacements are scaled using α = β (m) j = 1.For simplicity, we set λ β2i−1 = λ β2i = λ β .As for phonons, we take one mode for the donor molecule and two modes for the acceptor molecule in addition to the lattice phonon mode, and use parameters that approximately reproduce the experimentally observed phonon energies [27]: β2i−1 ≡ ω βD = 0.055, and ω (2) β2i ≡ ω βA2 = 0.12.Donor and acceptor molecules are specified by odd and even j's, respectively.For simplicity, we set The time-dependent vector potential A(t) for a pulse of an oscillating electric field is given by where ω pmp is the excitation energy, T pmp is the pulse width, and F is the electric field amplitude.The time-dependent Schrödinger equation for the exact many-electron wave function on the chain of N = 12 sites with periodic boundary condition is numerically solved by expanding the exponential evolution operator with a time slice dt = 0.02 eV −1 to the 15th order and by checking the conservation of the norm [13].The initial state is set in the electronic ground state.The classical equations for the lattice and molecular displacements are solved by the leapfrog method, where the forces are derived from the Hellmann-Feynman theorem.
Crystals 2012, 2 Figure 2 shows the ionicity as a function of the site energy difference ∆ near the phase boundary with different combinations of λ α and λ β .
It is shown that, as λ β increases, the neutral phase is stabilized, and the discontinuity in ionicity is enlarged, by increasing (decreasing) the ionicity in the ionic (neutral) phase on the small-∆ (large-∆) side of the phase boundary.As a consequence, in order for the ionic phase to be a typical Mott insulator with nearly one electron per site, λ β should be so large that the neutral phase is sufficiently stabilized.It is evident that, as λ α increases, the ionic phase is stabilized, and the discontinuity at the transition is suppressed.A finite λ α is necessary for dimerization and the three-dimensional ferroelectric order with a broken inversion symmetry.Both λ α and λ β are large in TTF-CA.
In the case of a large λ α and a large λ β , the neutral phase near the phase boundary is photoexcited with an energy just above the optical gap.The time evolution of the ionicity ρ(t) during and after photoexcitation is plotted in Figure 3.In this particular case, the phase boundary is located between  As ρ(t) increases, the electron density increases for the acceptor molecule and decreases for the donor molecule, so that the displacement increases at the acceptor molecule and decreases at the donor molecule.Thus, −v D (t) basically behaves as (−1) times the cosine function.The ionicity ρ(t) receives a positive feedback from molecular displacements and oscillates in the same phase with them.This −cosine behavior is consistent with the experimental observation [27].
There is a difference between the time profile in Figure 3 and that reported in [27].If we start from the neutral ground state without short-range lattice dimerization as a precursor, the lattice oscillation grows very gradually (not shown) after photoexcitation that breaks inversion symmetry.However, Crystals 2012, 2 63 the experiment shows a quick growth and oscillation of the lattice dimerization.If we start from the neutral ground state with a (small-amplitude) perfect dimerization without thermal fluctuations, an ultrafast charge transfer between neighboring donor and acceptor molecules continues to oscillate without dephasing, which is not realistic.In order to reproduce the time profile reported in [27], we need to include a small amplitude of lattice dimerization to induce a large-amplitude lattice oscillation immediately after photoexcitation.Then, we introduce in the initial state a dimerization (−1) j u j = 0.01, which is much smaller than (−1) j u j = 0.058 of the ground state on the ionic side of the phase boundary.In addition, we need to introduce random numbers according to the Boltzmann distribution at a finite temperature of T = 0.01 eV in u j , v (m) j , uj , and v(m) j of the initial state.Figure 4 shows the transient ionicity ρ(t), the displacement on the donor molecule −v D (t), and the spatially averaged dimerization (−1) j u j (t) in such a case.Now, the phase boundary is located between ∆ = 0.28 (ionic) and ∆ = 0.30 (neutral), and we use ∆ = 0.30.The dimerization (−1) j u j (t) rapidly increases with ρ(t) and oscillates like (−1) times the cosine function.Thus, ρ(t) receives a positive feedback from (−1) j u j already at an early stage as well as from −v D (t).The amplitude of the ρ(t)-oscillation is indeed larger than that without initial dimerization.
In summary, both lattice phonons and molecular vibrations contribute to the thermal and photoinduced neutral-ionic phase transitions.Lattice phonons stabilize the ionic phase through the spin-Peierls mechanism for dimerization and a finite spin gap.Molecular vibrations stabilize the neutral phase through the formation of a bipolaron lattice.The latter produce the large discontinuity in ionicity at the transition.Their contributions are manifested by the photoinduced oscillation of ionicity.The short-range lattice dimerization in the neutral phase promotes the photoinduced growth of a dimerized ionic domain.

Melting of Charge Order in θ-Type and α-Type BEDT-TTF Salts
Charge ordering in molecular conductors has been intensively studied both experimentally [40] and theoretically [41].The Coulomb interaction is a main driving force of charge ordering.It is not due to the Fermi surface nesting, so that it is different from a charge density wave.However, the lattice degrees of freedom must be included in the model to explain the structural deformation at the charge-ordering transition [42,43] and the photoinduced dynamics described below.
We use the two-dimensional three-quarter-filled extended Peierls-Holstein-Hubbard model where c † iσ creates an electron with spin σ at site i, n iσ =c † iσ c iσ , and n i = σ n iσ .The quantity u i,j denotes the intermolecular phonon's displacement, ui,j denotes its time derivative, and b † i creates a quantum Crystals 2012, 2 65 phonon of energy ω b , and g is the electron-molecular-vibration coupling strength.The other notations are standard and are introduced in [21].For instance, t i,j denotes the transfer integral for the bond between the neighboring ith and jth sites.Schematic illustrations of the high-and low-temperature, electronic and lattice structures of the conduction layers in the θand α-type salts are shown in Figure 5.For large systems (N = 12 × 12), we employ the Hartree-Fock approximation for the electronic states and use the parameter values in [20], the results of which are consistent with those for small systems (N = 12) with exact many-electron wave functions (g = 0) [21].When we use exact many-electron-phonon wave functions (g = 0), we use smaller systems (N = 8) and the parameter values in [21,22].Periodic boundary conditions are imposed on all of them.Photoexcitation is introduced in a similar manner to that in the previous section.The time evolution of the wave function and the lattice displacements is obtained by the method described in the previous section.
First, we discuss the different photoinduced melting dynamics of charge orders observed in θ-(BEDT-TTF) 2 RbZn(SCN) 4 and α-(BEDT-TTF) 2 I 3 [19].For this purpose, we ignore the electron-molecularvibration coupling by setting g = 0 for the moment and use large systems (N = 12 × 12).It is already clarified that the mechanisms for stabilizing the charge orders by lattice distortions are different in these two salts [20].In θ-(BEDT-TTF) 2 RbZn(SCN) 4 , the whole charge-rich (charge-poor) stripe is stabilized by strengthening (weakening) the horizontally connected bonds, as schematically shown in Figure 5.In α-(BEDT-TTF) 2 I 3 , the metallic phase without lattice distortion at high temperatures already possesses a charge-rich site B and a charge-poor site C from the kinetic origin.At low temperatures, the charge-rich site A and the charge-poor site A' bridged by the site B are locally stabilized by lattice distortion.Thus, local photoexcitations would easily weaken the charge order in the latter salt, while the charge order in the former salt would be robust.
Then, we investigate the growth of photoinduced domains with weakened charge order after artificially local photoexcitation as indicated by the dashed lines in the lower panels of Figure 5, by concentrating a finite vector potential on the bonds [(i x , i y ) = (7.5, 7), (7.5, 7.5), and (7.5, 8)] that connect four sites within a unit cell.Figure 6 shows, in the parenthesis, the ratio ∆E/∆E c , where ∆E is the increment in the total energy per site after local photoexcitation, and ∆E c is the critical increment above which the charge order is completely melted by global photoexcitation.Crystals 2012, 2

67
In Figure 6, the absolute values of the differences between the hole densities at t = 0 and those at t = 200, 400, and 600 are averaged over the direction parallel to the stripes and denoted by | ∆n i |: where | ψ(t) is the wave function at time t, L x = 12 is the number of sites along the axis parallel to the stripes, and i x (i y ) is the coordinate parallel (perpendicular) to the stripes.It gives a measure of how the photoinduced domain grows in the direction perpendicular to the stripes.In θ-(BEDT-TTF) 2 RbZn(SCN) 4 , the photoinduced domain remains localized near the place of photoexcitation, and hardly grows to the direction perpendicular to the stripes.This property prevents ∆E/∆E c from becoming large.For α-(BEDT-TTF) 2 I 3 , the photoinduced domain expands to the perpendicular direction.This result suggests that a macroscopic domain is much more easily created in the latter salt than in the former salt.
Next, we consider early-stage dynamics.Photoexcitations transfer electrons between neighboring molecules.Transfer integrals that are responsible for it are typically about 0.2 eV, which correspond to 20 fs in the time domain.Meanwhile, periods of lattice oscillations are typically about 1 ps.Thus, we expect that it takes about 1 ps for the lattice effect to appear in the electronic dynamics.Vibrations that have energy scales comparable with intermolecular transfer integrals are those inside a molecule, e.g., C=C stretching modes.Because these molecular vibrations and lattice phonons are coupled, the lattice effect on the electronic dynamics may appear much earlier than 1 ps through their couplings with molecular vibrations.Nonetheless those vibrations which strongly influence the early-stage electronic dynamics are intramolecular ones [22], as schematically shown in Figure 7 and as demonstrated below.

e -
We take an electron-molecular-vibration coupling into account.Because the electronic and vibronic energy scales are comparable, the Born-Oppenheimer approximation breaks down.Thus, we treat molecular vibrations quantum-mechanically (and lattice phonons classically).In numerical calculations, we compared results when molecular vibrations are quantum-mechanically treated and those when they are classically treated and found large differences especially in the electron-molecular-vibration interference pattern.This is in contrast to the case of the neutral-ionic transition in the previous section, where the transition is from an insulator to another insulator.When they are quantum-mechanically treated, the dimension of the wave function becomes huge.Then, we use a weak coupling (g = 0.0625 eV) and a large bare energy for the molecular vibration (ω b = 0.36 eV), which is comparable with the charge-transfer excitation energy in small systems (N = 8).
Figure 8 shows the time evolution of the hole density (red line) and the molecular displacement (green line) at molecule A, where the difference between the equilibrium values at the metal and insulator phases becomes the largest.The four panels sequentially show the evolution, which are separated for the purpose of explanation.Initially (top left panel), the oscillating electric field of pulse width T pmp = 5 fs is applied to the system around t = 0.It directly oscillates the hole density.Then, it forces the molecular displacement to oscillate indirectly through the electron-molecular vibration coupling.The sign of the displacement is so chosen Crystals 2012, 2 69 that it becomes large (small) when the hole density is large (small) in equilibrium.Thus, they are initially in phase.
About 50 fs after the photoexcitation (top right panel), they become out of phase.We calculated their dynamics with different parameters and frequently found their out-of-phase motions on this time scale.The energy supplied by the photoexcitation is much larger than the electron-molecular-vibration coupling.It can only partially absorb the energy and consequently makes the motions of the hole density and the molecular displacement out of phase.Soon after this (bottom left panel), their phases become mismatched, and finally (bottom right panel), the hole density follows the molecular displacement, the frequency of which is renormalized by electronic excitations to be smaller than the bare value.These results are consistent with the experimentally observed behavior [22].In reality, different frequencies of different intramolecular vibrations are coupled, with different strengths, to electrons in the highest occupied molecular orbital, so that the experimentally observed interference pattern is a little more complex.
In summary, both lattice phonons and molecular vibrations contribute to the photoinduced charge-order melting dynamics in the θ-type and α-type BEDT-TTF salts.On the time scale of lattice phonons, the different dynamics in these salts manifest the way by which the charge order is stabilized by lattice phonons is different in these salts.The charge order in the θ-type salt is robust, while that in the α-type salt is fragile.On the time scale of molecular vibrations, i.e., at an early stage, the charge dynamics and the vibrational dynamics interfere with each other.To reproduce the interference pattern theoretically, the quantum nature of molecular vibrations must be properly taken into account.

Mott-Insulator-to-Metal Transition in κ-Type BEDT-TTF Salts
Mott insulators are generally converted into metals either by weakening the effective on-site repulsion or by introducing carriers, as schematically shown in Figure 9.  Photoexcitation may be able to induce a Mott-insulator-to-metal transition via one of these pathways.In most situations, it is realized by the introduction of carriers [9,11,12].Note that carriers introduced by photoexcitations are different from those in a doped Mott insulator in the sense that equal amount of positively charged carriers (empty sites) and negatively charged carriers (doubly occupied sites) are introduced.Their binding through the attractive Coulomb interaction is suppressed by electron correlations, so that the conductivity spectra of photoexcited and doped Mott insulators are similar [10].Recently, it has been suggested that both pathways are realized by tuning the photoexcitation energy [15].It uses the intradimer molecular degrees of freedom in a dimer Mott insulator.
We use the two-dimensional three-quarter-filled extended Peierls-Hubbard model on the anisotropic triangular lattice schematically shown in Figure 10.
where c † iσ creates an electron with spin σ at site i, n iσ = c † iσ c iσ , and n i = σ n iσ .The operator b † ij creates a quantum phonon of energy ω ij , and g ij is the electron-phonon coupling strength.The other notations are standard and are introduced in [16].For instance, the intersite Coulomb repulsion V ij is assumed to be V ij = V 0 / | r i − r j | for the neighboring sites i and j at r i and r j .For simplicity, we consider only one mode for the creation operators b † ij , which modulate the intradimer transfer integrals t b1 .Thus, we have ω ij = ω b1 and g ij = g b1 .We take a high phonon energy ω b1 = 0.05 and a strong electron-phonon coupling g b1 = 0.06 to make the intradimer and interdimer charge-transfer (CT) bands overlap to a large extent.We use exact many-electron-phonon wave functions on small systems (N = 8) with periodic boundary condition and with the number of phonons restricted to a maximum of three at any b1 bond, and the parameter values in [16].Photoexcitation is introduced in a similar manner to that in the previous sections.The time evolution of the wave function is obtained by the method described in the previous sections.
The effective on-site Coulomb energy U dim is evaluated from the energies of the lowest one-and two-hole states for an isolated dimer with a transfer integral t b1 , on-site U and intersite V b1 repulsion strengths on the molecular bases.It is given by which are modulated by the displacement b b1 +b † b1 through the relation t b1 = t (0) b1 −g b1 b b1 +b † b1 .Owing to the molecular degrees of freedom inside a dimer, there are intradimer and interdimer CT excitations, as schematically shown in Figure 10.The force applied to the displacement b b1 + b † b1 depends linearly on the photoinduced difference in the expectation value σ c † iσ c jσ + c † jσ c iσ between sites i and j inside the dimer.It is analytically shown that this difference is insensitive to the photoexcitation energy (i.e., whether charge is transferred inside a dimer or between dimers) [16].As a consequence, any photoexcitation reduces σ c † iσ c jσ + c † jσ c iσ , reduces the magnitude of the intradimer transfer (t b1 < 0), and weakens the effective on-site repulsion U dim .
In order to see the modulation of the effective on-site repulsion U dim in a direct manner, we calculate the expectation value of the displacement b ij + b † ij as a function of time after photoexcitation.Its maximum decrement, −∆ b ij + b † ij , gives the maximum decrement in U dim , −∆U dim through Equation (8).We vary the electric field amplitude F and calculate the increment in the total energy ∆E divided by ω pmp , which corresponds to the number of absorbed photons.We show −∆U dim in Figure 11 as a function of ∆E/ω pmp for different ω pmp .Because the force applied to phonons is similar between the intradimer and interdimer CT processes, the ratio of −∆U dim to ∆E/ω pmp is almost independent of the excitation energy at least in the range of 0.3 < ω pmp < 0.5, which covers the intradimer and interdimer CT excitations.The effective on-site repulsion is therefore confirmed to be weakened to a similar extent irrespective of whether charge is transferred mainly within a dimer or mainly between dimers.
The number of carriers involved in the optical excitations up to ω prb is known to be proportional to the spectral weight obtained by the integration of the conductivity over 0 < ω < ω prb below the charge gap.We calculate the increment in the conductivity ∆σ (ω prb , t) = σ (ω prb , t) − σ (ω prb , −150), where with j being the current operator, a peak-broadening parameter set at 0.005, and It is averaged over 150 < t < 750, ∆σ (ω prb ) = (1/600) 750 150 ∆σ (ω prb , t)dt, and integrated over 0 < ω < ω prb , Figure 12 shows N (ω prb ), as a function of ∆E/ω pmp for different ω pmp .The energy ω prb is set at 0.01, which is well below the charge gap of 0.18 in the ground state.Although the quantity N (ω prb ) increases with ∆E/ω pmp for any ω pmp , its rate depends largely on ω pmp .For any ω prb below the charge gap, N (ω prb ) increases rapidly for ω pmp = 0.3, 0.35, and 0.4 and very slowly for ω pmp = 0.5.For ω pmp = 0.3, 0.35, and 0.4, the rates are close to each other.The number of carriers involved in the low-energy optical excitations is increased efficiently by 0.3 < ω pmp < 0.4, but it is negligibly increased for ω pmp = 0.5.This result shows that carriers introduced by photoexcitations with ω pmp near 0.3 have low excitation energies and are regarded as delocalized.These excitations are characterized as interdimer CT excitations.Although any CT excitation weakens U dim , it requires lattice motion and a long time.Consequently, if a Mott-insulator-to-metal transition is induced, it is mainly through the introduction of carriers.A photoexcitation with ω pmp = 0.5 introduces a negligible number of carriers.As a consequence, if a Mott-insulator-to-metal transition is induced, it is mainly through the weakening of U dim .This excitation is characterized as an intradimer CT excitation.
In summary, phonons modulating the transfer integral inside a dimer contribute to the realization of different transition pathways from the Mott insulator to the metal.In general, this metal-insulator transition can be induced by the weakening of effective interaction relative to the bandwidth or the introduction of carriers away from half filling.The phonons contribute to the weakening of the interaction irrespective of the photoexcitation energy, but this process is slow.The introduction of carriers is fast, but it strongly depends on the photoexcitation energy.

Conclusions
The research field of photoinduced phase transitions is developing rapidly.It is realized by the integration of the development of good target materials, the progress in experimental techniques including the improvement in time resolution of pump-probe spectroscopy, and that in theoretical calculations based on itinerant electron models.A single approach cannot generally cover different time scales of photoinduced dynamics in molecular conductors, where electrons and phonons are strongly correlated and their correlation effects appear differently on respective time scales.Theoretical works need to link closely with experimental studies on different probes with respective time scales, on which the relevant degrees of freedom depend.
Ideally speaking, it is desirable to have a continuous description from early-stage dynamics of individual motions of electrons and phonons to their collective motions, leading to a phase transition and relaxation to the initial phase.It may make it possible for us to construct an efficient transition pathway from the deterministic quantum regime to the stochastic classical regime.The situation is not so simple at present and possibly forever.It is very hard to approach such systems under nonequilibrium environments.In order to expand the possibilities for optical control of electronic phases, however, we need to clarify dynamical characteristics of quantum many-body systems.One of ultimate goals would be coherent control of photoinduced phase transitions.
For instance, electronic motion would generally interfere with various vibrations at an early stage, so that its effects on later dynamics or a possibility for making the transition more efficient should be pursued.In this context, a link between different hierarchies with respective time scales is one of the most important issues.Molecular materials naturally possess hierarchies, i.e., there exist intra-and inter-molecular degrees of freedom.So far, intermolecular degrees of freedom have mainly been focused.Photoinduced phase transitions proceeding on intra-and inter-molecular, electronic and structural stages will become important issues in the near future.From this viewpoint, the photoinduced dynamics in Pd(dmit) 2 metal complexes are currently studied [44].
The difficulty in understanding photoinduced dynamics comes from strong correlations between electrons and phonons.Nonetheless, its possibility for rich phenomena also comes from these strong correlations.Interplay between intra-and inter-molecular degrees of freedom and interplay among correlated electrons, lattice phonons and molecular vibrations will provide key points for developing photoinduced phase transitions in molecular conductors.
the mth mode on the jth molecule modulates the site energy with the coefficient β (m) j .The quantities uj and v(m) j are the time derivatives of u j and v (m) j , respectively.The parameters K α and K (m) βj are their elastic coefficients, and ω α and ω (m) βj are their bare phonon energies, respectively.

Figure 7 .
Figure 7. Schematic electronic and molecular structures of (BEDT-TTF) 2 X that are relevant to early-stage photoexcited dynamics.

Figure 9 .
Figure 9. Schematic diagram of Mott insulator and metal phases spanned by effective on-site repulsion (relative to bandwidth) and band filling.

Figure 11 .
Figure 11.Modulation of effective on-site repulsion −∆U dim , as a function of the number of absorbed photons ∆E/ω pmp , for different ω pmp From[16].Reproduced with permission from JPSJ.