THE SPITZER SPACE TELESCOPE OBSERVATION OF SPECTRA MOST LIKELY ATTRIBUTABLE TO DIVERSE AND ABUNDANT POPULATIONS OF POLYCYCLIC AROMATIC HYDROCARBONS (PAHS) IN SPACE HAS LED TO TREMENDOUS INTEREST IN THESE MOLECULES AS TRACERS OF THE PHYSICAL CONDITIONS IN DIFFERENT ASTROPHYSICAL REGIONS. A MAJOR CHALLENGE IN USING PAHS AS MOLECULAR TRACERS IS THE COMPLEXITY OF THE SPECTRAL FEATURES IN THE 3-20 M REGION. THE LARGE NUMBER AND VIBRATIONAL SIMILARITY OF THE PUTATIVE PAHS RESPONSIBLE FOR THESE SPECTRA NECESSITATE DETERMINATION FOR THE MOST ACCURATE BASIS SPECTRA POSSIBLE FOR COMPARISON. IT IS ESSENTIAL THAT THESE SPECTRA BE ESTABLISHED IN ORDER FOR THE REGIONS EXPLORED WITH THE NEWEST GENERATION OF OBSERVATORIES SUCH AS SOFIA AND JWST TO BE UNDERSTOOD. CURRENT STRATEGIES TO DEVELOP THESE SPECTRA FOR INDIVIDUAL PAHS INVOLVE EITHER MATRIX-ISOLATION IR MEASUREMENTS OR QUANTUM CHEMICAL CALCULATIONS OF HARMONIC VIBRATIONAL FREQUENCIES. THESE STRATEGIES HAVE BEEN EMPLOYED TO DEVELOP THE SUCCESSFUL PAH IR SPECTRAL DATABASE AS A REPOSITORY OF BASIS FUNCTIONS USED TO FIT ASTRONOMICALLY OBSERVED SPECTRA BUT THEY ARE LIMITED IN IMPORTANT WAYS. BOTH TECHNIQUES PROVIDE AN ADEQUATE DESCRIPTION OF THE MOLECULES IN THEIR ELECTRONIC VIBRATIONAL AND ROTATIONAL GROUND STATE BUT THESE CONDITIONS DO NOT REPRESENT ENERGETICALLY HOT REGIONS FOR PAHS NEAR STRONG RADIATION FIELDS OF STARS AND ARE NOT DIRECT REPRESENTATIONS OF THE GAS PHASE. SOME NON-NEGLIGIBLE MATRIX EFFECTS ARE KNOWN IN CONDENSED-PHASE STUDIES AND THE INCLUSION OF ANHARMONICITY IN QUANTUM CHEMICAL CALCULATIONS IS ESSENTIAL TO GENERATE PHYSICALLY-RELEVANT RESULTS ESPECIALLY FOR HOT BANDS. WHILE SCALING FACTORS IN EITHER CASE CAN BE USEFUL THEY ARE AGNOSTIC TO THE SYSTEM STUDIED AND ARE NOT ROBUSTLY PREDICTIVE. ONE STRATEGY THAT HAS EMERGED TO CALCULATE THE MOLECULAR VIBRATIONAL STRUCTURE USES VIBRATIONAL PERTURBATION THEORY ALONG WITH A QUARTIC FORCE FIELD (QFF) TO ACCOUNT FOR HIGHER-ORDER DERIVATIVES OF THE POTENTIAL ENERGY SURFACE. QFFS CAN REGULARLY PREDICT THE FUNDAMENTAL VIBRATIONAL FREQUENCIES TO WITHIN 5 CM-1 OF EXPERIMENTALLY MEASURED VALUES. THIS LEVEL OF ACCURACY REPRESENTS A REDUCTION IN DISCREPANCIES BY AN ORDER OF MAGNITUDE COMPARED WITH HARMONIC FREQUENCIES CALCULATED WITH DENSITY FUNCTIONAL THEORY (DFT). THE MAJOR LIMITATION OF THE QFF STRATEGY IS THAT THE LEVEL OF ELECTRONIC-STRUCTURE THEORY REQUIRED TO DEVELOP A PREDICTIVE FORCE FIELD IS PROHIBITIVELY TIME CONSUMING FOR MOLECULAR SYSTEMS LARGER THAN 5 ATOMS. RECENT ADVANCES IN QFF TECHNIQUES UTILIZING INFORMED DFT APPROACHES HAVE PUSHED THE SIZE OF THE SYSTEMS STUDIED UP TO 24 HEAVY ATOMS BUT RELEVANT PAHS CAN HAVE UP TO HUNDREDS OF ATOMS. WE HAVE DEVELOPED ALTERNATIVE ELECTRONIC-STRUCTURE METHODS THAT MAINTAIN THE ACCURACY OF THE COUPLED-CLUSTER CALCULATIONS EXTRAPOLATED TO THE COMPLETE BASIS SET LIMIT WITH RELATIVISTIC AND CORE CORRELATION CORRECTIONS APPLIED: THE CCCR QFF. THESE ALTERNATIVE METHODS ARE BASED ON SIMPLIFICATIONS OF HARTREE FOCK THEORY IN WHICH THE COMPUTATIONALLY INTENSIVE TWO-ELECTRON INTEGRALS ARE APPROXIMATED USING EMPIRICAL PARAMETERS. THESE METHODS REDUCE COMPUTATIONAL TIME TO ORDERS OF MAGNITUDE LESS THAN THE CCCR CALCULATIONS. WE HAVE DERIVED A SET OF OPTIMIZED EMPIRICAL PARAMETERS TO MINIMIZE THE DIFFERENCE BETWEEN THE CCCR QFF ENERGIES FOR THE GEOMETRIES USED TO BUILD THE FORCE FIELD FOR TWO MOLECULAR IONS OF ASTROCHEMICAL SIGNIFICANCE. WE HAVE SHOWN THAT IT IS POSSIBLE TO DERIVE A SET OF EMPIRICAL PARAMETERS THAT WILL PRODUCE RMS ENERGY DIFFERENCES OF LESS THAN 2 CM-1 FOR OUR TEST SYSTEMS. WE ARE PROPOSING TO ADOPT THIS REPARAMETERIZATION STRATEGY AND SOME OF THE LESSONS LEARNED FROM THE INFORMED DFT STUDIES TO CREATE A SEMI-EMPIRICAL METHOD WHOSE TREMENDOUS SPEED WILL ALLOW US TO STUDY THE ROVIBRATIONAL STRUCTURE OF LARGE PAHS WITH UP TO 100S OF CARBON ATOMS.