A potential energy force field designed for modeling nucleic acids and particularly their complexes with cationic ligands is presented. The force field is a modified version of that developed by Weiner, S.J., Kollman, P.A., Nguyen, D.T. and Case, D.A.,J. Comp. Chem. 7,230-252 (1986) and is based upon the use of a distance dependent dielectric constant, epsilon = 4rij, and partially neutralized phosphates to represent solvent and counterion. Changes from the Weiner et al. force field include additional atom types and modifications to van der Waals, electrostatic, hydrogen bonding and torsional parameters. Molecular modeling test cases of the force field are presented for a number of simple small molecules, as well as uracil and benzene dimerization, thymine-adenine and cytosine-guanine base pair formation, and adenosine/deoxyadenosine pseudorotation. Several DNA and RNA oligomers and DNA/RNA intercalation complexes with ethidium are also modeled with the force field. In all cases, the modeling results compare favorably with available experimental results. Additionally, conformational trends observed experimentally for nucleic acids by NMR and X-ray crystallographic techniques are reproduced. The modeling results for ethidium intercalation indicate a complex in which the favorable interactions are primarily van der Waals contacts, and in which electrostatic interactions are a relatively minor component. We feel the force field is particularly useful for molecular mechanics aided drug design, and an analysis of modeling results with respect to design of drugs which bind selectively to RNA is presented.