Biomedical subjects

B Lasorne

Publications and source records attributed to B Lasorne.

3 recordsLinked to original sources

PubMed2005

Wave packets in a bifurcating region of an energy landscape: Diels-Alder dimerization of cyclopentadiene.

Quantum dynamics in a valley ridge inflection (VRI) point region is analyzed in the case of the Diels-Alder endo-dimerization of cyclopentadiene pointed out recently by [Caramella et al., J. Am. Chem. Soc. 124, 1130 (2002)]. The VRI point is located along the reaction path connecting the bispericyclic symmetrical transition structure put in evidence by Caramella et al. and the transition state of the Cope rearrangement. Dynamics is carried out by using constrained Hamiltonian methodology. The active coordinates are the first formed C-C bond length and the difference between the two other C-C bond lengths which achieve the dimerization as 4+2 or 2+4 adducts. A two-dimensional (2D) minimum-energy surface have been computed at the Becke 3 Lee-Yong-Parr6-31G* level. The energy landscape can be classified as an uphill ridge-pitchfork VRI bifurcation according to a recent classification of bifurcation events [W. Quapp, J. Mol. Struct. 695-696, 95 (2004)]. Dynamics does not describe the thermal reaction but concerns wave packets which could be prepared by pulse reagents, i.e., by coherent control. We analyze how the shape and initial location on the ground potential-energy surface are linked to the synchronous or asynchronous mechanism of the final step after the first transition state. We use a one-dimensional model of optimum control theory to check the feasibility of such a coherent preparation. The wave-packet evolution in the VRI domain is well explained by semiclassical predictions even with the negative curvature of the unstable ridge. Finally, a crude model of dissipation has been introduced to test the stability of the 2D predictions.

Journal Article↗

PubMed2004

Cumulative isomerization probability studied by various transition state wave packet methods including the MCTDH algorithm. Benchmark: HCN-->CNH isomerization.

The 3D cumulative isomerization probability N(E) for the transfer of a light particle between two atoms is computed by one time-independent and two time-dependent versions of the transition state wave packet (TSWP) method. The time-independent method is based on the direct expansion of the microcanonical projection operator on Chebyshev polynomials. In the time-dependent TSWP methods, the propagations are carried out by the split operator scheme and the multiconfiguration time-dependent Hartree (MCTDH) algorithm. This is the very first implementation of the TSWP method in the Heidelberg MCTDH package [G. W. Worth, M. H. Beck, A. Jackle, and H.-D. Meyer, The MCDTH package, Version 8.2 (2000); H.-D Meyer, Version 8.3 (2002). See http://www.pci.uni-heidelberg.de/tc/usr/mctdh/]. The benchmark is the HCN-->CNH isomerization for zero total angular momentum. Particular insights are given into the tunneling region. In larger systems, the time-dependent version of TSWP making use of the MCTDH algorithm will permit to treat more and more modes quantum mechanically, for very accurate results. Therefore, it was important to calibrate the implementation. Besides, we also assess the efficiency of a reduced dimensionality approach by comparing the new exact 3D calculations of N(E) for the HCN-->CNH isomerization with results obtained via 1D or 2D active subspaces. This suggests that, it should be possible to take directly benefit of the present 3D approaches, adapted for triatomic Jacobi coordinates to compute N(E) for H-transfer in larger systems, via three active coordinates. The prerequisite is then the simplification of the reduced 3D kinetic energy operator with rigid constraint to take the form corresponding to a pseudo triatomic system in Jacobi coordinates with effective masses. This last step is checked in the methoxy radical and malonaldehyde. Finally, different ways to obtain reliable eigenvectors of the flux operator associated with a dividing surface are revisited.

Journal Article↗

PubMed2004

Nonadiabatic interactions in wave packet dynamics of the bromoacetyl chloride photodissociation.

The competitive photodissociation of bromoacetyl chloride BrCH2COCl in the first 1A" state (S1) by 248 nm photons is investigated by nonadiabatic wave packet simulations. We show that the preferential breaking of the stronger C-Cl bond (alpha to the excited carbonyl) over the weaker C-Br bond (beta) could be explained by a diabatic trapping or nonadiabatic recrossing as previously proposed. Our energy resolved flux analysis agrees fairly well with the experimental branching ratio (C-Cl:C-Br=1.0:0.4). Even if this does not prove the mechanism, this at least prevents to discard it. A reduced dimensionality approach based on constrained Hamiltonian is used. The nonadiabatic dissociation is studied in the two C-O/C-X (X=Br, Cl) subspaces to emphasize the role of the C-O vibration upon [nO-->piCO*] excitation. The internal torsion and wagging dihedral angles are frozen at their Franck-Condon value, according to preliminary dynamical tests. The other inactive coordinates are optimized at the trans and Cs constrained geometry in the first excited state. Corresponding 2D cuts in the potential energy surfaces have been computed at the CASSCF level. The nonadiabatic kinetic couplings are highly peaked along an avoided crossing seam in both cases. A two-state diabatic model with a constant potential coupling is proposed in the two C-O/C-X subspaces. The inclusion of the C-O stretching in the active coordinates improves the value of the branching ratio over our previous 1D computation.

Journal Article↗