Coupled-Channels Density Matrix (CCDM)

• Coupled-Channels Density Matrix (CCDM) is an open quantum system approach that models low-energy nuclear reactions by partitioning essential intrinsic states from a complex many-body environment.
• It employs a Lindblad-type master equation to incorporate both dissipation and quantum decoherence, providing a more complete description than traditional absorptive potentials.
• Applications to reactions such as 12C+12C and 16O+154Sm illustrate how CCDM affects tunnelling, fusion, and angular distributions by modifying interference and flux loss.

Searching arXiv for the specified CCDM papers and closely related work. arxiv_search(query="Coupled-channels density-matrix low-energy nuclear reaction dynamics Diaz-Torres", max_results=10, sort_by="relevance") Coupled-Channels Density Matrix (CCDM) denotes an open-quantum-system formulation of low-energy nuclear reaction dynamics in which the relative motion of projectile and target, together with a truncated set of intrinsic collective states, is treated as a reduced system, while excluded many-body configurations are treated as an environment. In this framework the density operator evolves according to a Liouville–von Neumann equation supplemented by Lindblad dissipative terms, so that dissipation and quantum decoherence are incorporated simultaneously rather than through an absorptive potential alone. The approach was developed by Diaz-Torres to quantify how decoherence modifies scattering, tunnelling, fusion, and inelastic observables in systems such as $^{12}\mathrm{C}+{}^{12}\mathrm{C}$ and $^{16}\mathrm{O}+{}^{154}\mathrm{Sm}$ (Diaz-Torres, 2010, Diaz-Torres, 2010).

1. Open-system partitioning in low-energy nuclear reactions

At low energy, near the ground state, collective modes such as vibration and rotation dominate the intrinsic structure of atomic nuclei. Conventional coupled-channels calculations therefore employ a truncated basis of low-lying collective states $\{|\phi_i\rangle\}$ as the relevant intrinsic degrees of freedom. CCDM retains this coupled-channels logic but reformulates it in explicitly open-system terms: the reduced system $S$ consists of the relative projectile–target motion plus a small set of “important” intrinsic states, whereas the environment $E$ consists of all remaining many-body states, including high-lying single-particle excitations, doorway states, transfer channels, continuum breakup states, and compound-nucleus molecular resonances (Diaz-Torres, 2010).

The total Hilbert space is written as

$H_{\mathrm{total}} = H_S \otimes H_E,$

with total Hamiltonian

$$H = H_S \otimes 1_E + 1_S \otimes H_E + V_{\mathrm{int}},$$

where $V_{\mathrm{int}}$ couples the reduced system to the environment. In this partitioning, the nuclear collision is not treated as an isolated coupled-channels problem but as an open quantum system whose reduced dynamics can become mixed through environmental coupling.

This construction directly addresses a limitation of the standard truncated coupled-channels basis. Excluded states are often essential, but in traditional treatments their influence is represented through complex potentials. CCDM was formulated around the question of whether this is a complete description of open-system dynamics, especially when the loss of phase coherence among channels may itself be observable (Diaz-Torres, 2010).

2. Master-equation structure and coupled-channels representation

The formal starting point is the Liouville–von Neumann equation for the full density operator $\rho_{\mathrm{tot}}(t)$,

$$i\hbar \frac{d\rho_{\mathrm{tot}}}{dt} = [H,\rho_{\mathrm{tot}}],$$

followed by a partial trace over the environment,

$^{16}\mathrm{O}+{}^{154}\mathrm{Sm}$0

Under the standard Born–Markov and secular approximations, the reduced density operator obeys a Lindblad-type master equation,

$^{16}\mathrm{O}+{}^{154}\mathrm{Sm}$1

where $^{16}\mathrm{O}+{}^{154}\mathrm{Sm}$2 includes Lamb-shift corrections and the dissipator has the canonical Lindblad form

$^{16}\mathrm{O}+{}^{154}\mathrm{Sm}$3

The operators $^{16}\mathrm{O}+{}^{154}\mathrm{Sm}$4 mediate decoherence and dissipation channels; examples given in the formulation include projectors onto collective-state subspaces and radial-position-dependent operators that absorb flux into the environment (Diaz-Torres, 2010).

In a product basis $^{16}\mathrm{O}+{}^{154}\mathrm{Sm}$5, the reduced density matrix is expanded as

$^{16}\mathrm{O}+{}^{154}\mathrm{Sm}$6

which yields coupled evolution equations of the form

$^{16}\mathrm{O}+{}^{154}\mathrm{Sm}$7

Here $^{16}\mathrm{O}+{}^{154}\mathrm{Sm}$8 is the radial kinetic operator and $^{16}\mathrm{O}+{}^{154}\mathrm{Sm}$9 are the usual coupled-channels couplings (Diaz-Torres, 2010).

The 2010 extension introduced full angular-momentum couplings by expanding the density operator in the basis

$\{|\phi_i\rangle\}$0

tensored with a discretized radial grid $\{|\phi_i\rangle\}$1. The coherent Hamiltonian is written as $\{|\phi_i\rangle\}$2, with matrix elements

$\{|\phi_i\rangle\}$3

In this representation, $\{|\phi_i\rangle\}$4 is the target spin, $\{|\phi_i\rangle\}$5 the relative orbital angular momentum, $\{|\phi_i\rangle\}$6 the total angular momentum, $\{|\phi_i\rangle\}$7 the bare channel potential, $\{|\phi_i\rangle\}$8 the coupling potential, and $\{|\phi_i\rangle\}$9 the intrinsic excitation energy. The dissipative Liouvillian is then added through physically motivated jump operators associated with local decay rates $S$0, including couplings to auxiliary environmental states $S$1 (Diaz-Torres, 2010).

3. Relation to conventional coupled-channels and optical-model descriptions

CCDM reduces to the conventional coupled-channels formalism when all Lindblad operators are set to zero, equivalently when $S$2. In that limit the reduced density operator undergoes unitary coupled-channels evolution. For a pure-state ansatz $S$3, one recovers the usual coupled-channels Schrödinger equation for amplitudes $S$4,

$S$5

CCDM is therefore not a replacement for coupled channels in the narrow sense; it is an extension that preserves the coherent limit exactly (Diaz-Torres, 2010).

The principal distinction from the traditional optical-model treatment lies in the status of excluded channels. In standard coupled-channels or optical-model calculations one often adds an imaginary potential $S$6 to simulate loss of flux into unobserved channels. This accounts for dissipation, but it leaves the reduced description in a coherent pure-state form. In the Lindblad formulation, the corresponding losses arise from the jump operators and lead to a genuine mixed state with

$S$7

Thus, absorption and decoherence are not identified with each other: dissipation removes probability from the reduced system, while decoherence suppresses off-diagonal density-matrix elements $S$8 with $S$9 and thereby destroys channel interference (Diaz-Torres, 2010).

A common misconception is that an imaginary potential provides a complete effective description of open-system dynamics. CCDM was formulated precisely to test that assumption. The explicit Lindblad treatment shows that flux loss can coexist with preserved purity in the purely absorptive case, whereas environmental couplings can additionally induce true mixing of the reduced state (Diaz-Torres, 2010, Diaz-Torres, 2010).

4. Decoherence diagnostics and reaction observables

A central diagnostic in CCDM is the coherence ratio

$E$0

which equals $E$1 for a pure state and decreases below $E$2 as decoherence develops. This scalar measure was used in the time-dependent calculations to discriminate situations in which dissipation is present but coherence is largely preserved from those in which environmental couplings rapidly destroy coherence (Diaz-Torres, 2010).

In the density-matrix formulation, asymptotic observables are extracted after propagation to times $E$3 or, in the numerical implementation, to a final time $E$4 at which outgoing wave packets have separated. The overview formulation gives representative expressions such as the elastic probability,

$E$5

and the fusion or tunnelling probability,

$E$6

Inelastic cross sections are similarly expressed through projectors onto the relevant channels. When decoherence suppresses $E$7 for $E$8, interference among alternative pathways is reduced, and the energy dependence of $E$9 is modified (Diaz-Torres, 2010).

The full angular-momentum implementation also provides asymptotic scattering observables in terms of the projected density matrix. After projection onto a target state $H_{\mathrm{total}} = H_S \otimes H_E,$0 and a scattering direction $H_{\mathrm{total}} = H_S \otimes H_E,$1, the differential probability $H_{\mathrm{total}} = H_S \otimes H_E,$2 is expressed through spherical harmonics, Clebsch–Gordan coefficients, and the asymptotic matrix

$H_{\mathrm{total}} = H_S \otimes H_E,$3

while total state populations are obtained by angular integration and summation over $H_{\mathrm{total}} = H_S \otimes H_E,$4 (Diaz-Torres, 2010).

Because CCDM is implemented with incident wave packets, energy-resolved information is not taken as primitive but reconstructed. The 2010 formulation introduced a window-operator method $H_{\mathrm{total}} = H_S \otimes H_E,$5 and defined energy-binned reflection matrices

$H_{\mathrm{total}} = H_S \otimes H_E,$6

thus permitting direct comparison with stationary coupled-channels results while retaining decoherence effects (Diaz-Torres, 2010).

5. Model calculations and representative applications

A one-dimensional illustration summarized in the 2010 overview contrasts three situations: coherent evolution, optical absorption through an imaginary pocket, and Lindblad dissipation. In that example, wave-packet scattering off a barrier with an imaginary pocket preserves purity, with

$H_{\mathrm{total}} = H_S \otimes H_E,$7

at all times, whereas Lindblad dissipation yields $H_{\mathrm{total}} = H_S \otimes H_E,$8 together with a marked reduction and reshaping of the tunnelling probability $H_{\mathrm{total}} = H_S \otimes H_E,$9. The result isolates decoherence from mere absorption and shows that the two need not have the same phenomenological consequences (Diaz-Torres, 2010).

For $$H = H_S \otimes 1_E + 1_S \otimes H_E + V_{\mathrm{int}},$$0, a realistic two-center shell-model calculation summarized in the same work indicates that when two oblate $$H = H_S \otimes 1_E + 1_S \otimes H_E + V_{\mathrm{int}},$$1 nuclei overlap, their single-particle molecular levels depend sensitively on relative alignment. Different alignments open distinct “giant-resonance-like” environments, including butterfly and twisting modes, which couple to low-lying rotational states and induce decoherence. In this setting, separate Lindblad operators $$H = H_S \otimes 1_E + 1_S \otimes H_E + V_{\mathrm{int}},$$2 may be introduced for different alignment-dependent environments. The resulting loss of coupled-channels coherence among long-range Coulomb-excited rotational channels is described as smoothing out narrow molecular-resonance peaks in the fusion excitation function and offering a mechanism for the resonant structures in the astrophysical $$H = H_S \otimes 1_E + 1_S \otimes H_E + V_{\mathrm{int}},$$3-factor for $$H = H_S \otimes 1_E + 1_S \otimes H_E + V_{\mathrm{int}},$$4 (Diaz-Torres, 2010).

The most explicit application in the technical development is the low-energy collision $$H = H_S \otimes 1_E + 1_S \otimes H_E + V_{\mathrm{int}},$$5. In that calculation the projectile $$H = H_S \otimes 1_E + 1_S \otimes H_E + V_{\mathrm{int}},$$6 is treated as inert, while the deformed target $$H = H_S \otimes 1_E + 1_S \otimes H_E + V_{\mathrm{int}},$$7 is represented by the $$H = H_S \otimes 1_E + 1_S \otimes H_E + V_{\mathrm{int}},$$8, $$H = H_S \otimes 1_E + 1_S \otimes H_E + V_{\mathrm{int}},$$9, and $V_{\mathrm{int}}$0 states of its ground-state band, with energies $V_{\mathrm{int}}$1, $V_{\mathrm{int}}$2, and $V_{\mathrm{int}}$3. These are coupled by a macroscopic deformed Woods–Saxon form with all-order nuclear couplings and Coulomb couplings to second order in $V_{\mathrm{int}}$4 and first order in $V_{\mathrm{int}}$5. The coherent couplings $V_{\mathrm{int}}$6 are generated by the code CCFULL (Diaz-Torres, 2010).

Two dissipative environments are then added. The first is a “fusion” environment coupling all channels to an internal pocket via a Fermi-shaped decay function $V_{\mathrm{int}}$7 of depth $V_{\mathrm{int}}$8 at pocket radius $V_{\mathrm{int}}$9. The second is a “surface” environment coupling only the target ground state to auxiliary states through a Gaussian decay of width $\rho_{\mathrm{tot}}(t)$0 at a nuclear-contact radius of approximately $\rho_{\mathrm{tot}}(t)$1, modelling complex transfer processes. The calculated coherence ratio $\rho_{\mathrm{tot}}(t)$2 shows that the fusion environment alone leaves $\rho_{\mathrm{tot}}(t)$3, whereas adding the surface environment causes a rapid drop in $\rho_{\mathrm{tot}}(t)$4, indicating strong decoherence. At incident energy $\rho_{\mathrm{tot}}(t)$5, inclusion of the surface environment shifts minima in the $\rho_{\mathrm{tot}}(t)$6, $\rho_{\mathrm{tot}}(t)$7, and $\rho_{\mathrm{tot}}(t)$8 angular distributions by a few degrees and reduces the amplitude of the interference pattern. The reported total populations further show that decoherence from the surface environment reduces $\rho_{\mathrm{tot}}(t)$9 and $$i\hbar \frac{d\rho_{\mathrm{tot}}}{dt} = [H,\rho_{\mathrm{tot}}],$$0 excitation probabilities and the fusion yield while increasing the elastic fraction (Diaz-Torres, 2010).

6. Physical implications, scope, and interpretation

In the CCDM framework, quantum decoherence is a dynamical process by which a pure coupled-channels superposition evolves into a mixed state, reducing interference among fusion, fission, and inelastic pathways. The formalism was presented as the first fully quantum-mechanical framework that simultaneously couples relative motion and intrinsic states, treats excluded many-body modes as an environment, and describes both dissipation and decoherence on the same footing (Diaz-Torres, 2010).

The immediate physical consequence is that reaction observables may depend not only on channel couplings and absorptive strength but also on the extent to which environmental interactions suppress coherence between amplitudes. This is most explicit for sub-barrier tunnelling and fusion, where the coherent coupled-channels mechanism typically enhances penetration. When decoherence is active, that enhancement is diminished, and the energy dependence of tunnelling or fusion probabilities is correspondingly altered (Diaz-Torres, 2010).

The $$i\hbar \frac{d\rho_{\mathrm{tot}}}{dt} = [H,\rho_{\mathrm{tot}}],$$1 calculations further show that decoherence can modify angular distributions through the suppression of partial-wave interference and can redistribute asymptotic yields between elastic, inelastic, and fusion sectors. This suggests that transfer-like surface environments may be important whenever coherent coupled-channels calculations reproduce mean trends but miss detailed interference structure or overestimate inelastic and fusion probabilities (Diaz-Torres, 2010).

Within the 2010 program, the approach was highlighted as relevant both to the astrophysically important $$i\hbar \frac{d\rho_{\mathrm{tot}}}{dt} = [H,\rho_{\mathrm{tot}}],$$2 system and to the broader objective of a unified quantum dynamical description of reaction processes involving weakly-bound nuclei (Diaz-Torres, 2010). A plausible implication is that CCDM is best understood not merely as a numerical variant of coupled channels, but as a conceptual extension that separates three effects often conflated in phenomenological practice: coherent channel coupling, dissipative probability loss, and genuine decoherence.