In this paper we present a theoretical and computational study of the temporal dynamics and energetics of Coulomb explosion of (CD4)(n) and (CH4)(n) (n=55-4213) molecular heteroclusters in ultraintense (I=10(16)-10(19) W cm(-2)) laser fields, addressing the manifestation of electron dynamics, together with nuclear energetic and kinematic effects on the heterocluster Coulomb instability. The manifestations of the coupling between electron and nuclear dynamics were explored by molecular dynamics simulations for these heteroclusters coupled to Gaussian laser fields (pulse width tau=25 fs), elucidating outer ionization dynamics, nanoplasma screening effects (being significant for Iless than or equal to10(17) W cm(-2)), and the attainment of cluster vertical ionization (CVI) (at I=10(17) W cm(-2) for cluster radius R(0)less than or equal to31 A). Nuclear kinematic effects on heterocluster Coulomb explosion are governed by the kinematic parameter eta=q(C)m(A)/q(A)m(C) for (CA(4))(n) clusters (A=H,D), where q(j) and m(j) (j=A,C) are the ionic charges and masses. Nonuniform heterocluster Coulomb explosion (eta>1) manifests an overrun effect of the light ions relative to the heavy ions, exhibiting the expansion of two spatially separated subclusters, with the light ions forming the outer subcluster at the outer edge of the spatial distribution. Important features of the energetics of heterocluster Coulomb explosion originate from energetic triggering effects of the driving of the light ions by the heavy ions (C4+ for I=10(17)-10(18) W cm(-2) and C6+ for I=10(19) W cm(-2)), as well as for kinematic effects. Based on the CVI assumption, scaling laws for the cluster size (radius R-0) dependence of the energetics of uniform Coulomb explosion of heteroclusters (eta=1) were derived, with the size dependence of the average (E-j,E-av) and maximal (E-j,E-M) ion energies being E-j,E-av=aR(0)(2) and E-j,E-M=(5a/3)R-0(2), as well as for the ion energy distributions P(E-j)proportional toE(j)(1/2); E(j)less than or equal toE(j,M). These results for uniform Coulomb explosion serve as benchmark reference data for the assessment of the effects of nonuniform explosion, where the CVI scaling law for the energetics still holds, with deviations of the a coefficient, which increase with increasing eta. Kinematic effects (for eta>1) result in an isotope effect, predicting the enhancement (by 9%-11%) of E-H,E-av for Coulomb explosion of (C4+H4+)(eta) (eta=3) relative to E-D,E-av for Coulomb explosion of (C4+D4+)(eta) (eta=1.5), with the isotope effect being determined by the ratio of the kinematic parameters for the pair of Coulomb exploding clusters. Kinematic effects for nonuniform explosion also result in a narrow isotope dependent energy distribution (of width DeltaE) of the light ions (with DeltaE/E(H,av)similar or equal to0.3 and DeltaE/E(D,av)similar or equal to0.4), with the distribution peaking at the high energy edge, in marked contrast with the uniform explosion case. Features of laser-heterocluster interactions were inferred from the analyses of the intensity dependent boundary radii (R-0)(I) and the corresponding average D+ ion energies (E-D,E-av)(I), which provide a measure for optimization of the cluster size at intensity I for the neutron yield from dd nuclear fusion driven by Coulomb explosion (NFDCE) of these heteroclusters. We infer on the advantage of deuterium containing heteronuclear clusters, e.g., (CD4)(n) in comparison to homonuclear clusters, e.g., (D-2)(n/2), for dd NFDCE, where the highly charged heavy ions (e.g., C4+ or C6+) serve as energetic and kinematic triggers driving the D+ ions to a high (10-200 keV) energy domain. (C) 2004 American Institute of Physics.